Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 4,183,295 A→C N684T (AAC→ACC)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi4,183,2950AC96.5% 491.8 / 2.6 141N684T (AAC→ACC) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base C (66/70);  total (69/72)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 141 total aligned reads displayed.

TTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATC  >  NC_000913_3_pae_pgi/4183149‑4183440
                                                                                                                                                  |                                                                                                                                                 
          aaaGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAg                                                                                                                                      >  2:1078875/1‑150 (MQ=255)
           aaGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAgg                                                                                                                                     <  2:1119004/150‑1 (MQ=255)
           aaGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAgg                                                                                                                                     <  2:574199/150‑1 (MQ=255)
                  aTCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTcc                                                                                                                              >  4:419232/1‑150 (MQ=255)
                    ccAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGa                                                                                                                            >  4:333536/1‑150 (MQ=255)
                    ccAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATACAATGTCAGGCCGTTccta                                                                                                                            <  4:606300/150‑3 (MQ=255)
                     cAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGAc                                                                                                                           >  4:54102/1‑150 (MQ=255)
                      aGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTcctgga                                                                                                                          >  3:876036/1‑146 (MQ=255)
                       gCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGGATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGActc                                                                                                                         >  2:1585946/1‑150 (MQ=255)
                          tGTTCAGCCGCGACCAGGGTGAATACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGAGCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTgc                                                                                                                      <  2:431665/150‑1 (MQ=255)
                           gTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTgcg                                                                                                                     >  2:1160881/1‑150 (MQ=255)
                              cAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGt                                                                                                                                         >  4:146074/1‑127 (MQ=255)
                              cAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGt                                                                                                                                         <  3:146074/127‑1 (MQ=255)
                              cAGCCGCGACCAGGTTGACTACATGGACGTACCCACCCNGCAGGTNGNATCCGTCGGTGCGTCCCTGAACCCGTACCTGGAACACGACGACGCCAACCGGGCATTGATGGGAGCGACCAGGGAACGTCAGGCCGTCCCGACTCGgcgcgc                                                                                                                  >  3:468677/1‑150 (MQ=255)
                               aGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTgcgcgcc                                                                                                                 >  2:330088/1‑149 (MQ=255)
                                  cgcgACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGAt                                                                                                              <  2:1219513/150‑1 (MQ=255)
                                      aCCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGc                                                                                                          <  3:877192/150‑1 (MQ=255)
                                       ccAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGcc                                                                                                         >  2:804457/1‑150 (MQ=255)
                                          ggTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCt                                                                                                      >  2:1349819/1‑150 (MQ=255)
                                             tGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCNACGTCAGGCCGTTCCGCCTCTGCGCGCTGATAAGCCGCtggt                                                                                                   <  1:675488/150‑1 (MQ=255)
                                             tGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggt                                                                                                   <  2:239204/150‑1 (MQ=255)
                                              gACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggtt                                                                                                  <  1:82123/150‑1 (MQ=255)
                                              gACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggtt                                                                                                  <  1:1627716/150‑1 (MQ=255)
                                                 tACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAgg                                                                                                                                     <  2:743015/112‑1 (MQ=255)
                                                 tACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAgg                                                                                                                                     >  1:743015/1‑112 (MQ=255)
                                                  aCATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTa                                                                                              <  4:642120/150‑1 (MQ=255)
                                                  aCATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTa                                                                                              <  3:333536/150‑1 (MQ=255)
                                                    aTGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCa                                                                                                                                                 <  3:636556/97‑1 (MQ=255)
                                                    aTGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACt                                                                                            >  1:770847/1‑150 (MQ=255)
                                                    aTGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACt                                                                                            >  1:296911/1‑150 (MQ=255)
                                                     tGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTg                                                                                           <  2:350140/150‑1 (MQ=255)
                                                      ggACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTgg                                                                                          <  3:93606/150‑1 (MQ=255)
                                                       gACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGt                                                                                         <  3:63303/150‑1 (MQ=255)
                                                        aCGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTATCCGGTAATACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTTATAAGCCGCTGGTTGGTACTGGTa                                                                                        <  3:659451/150‑1 (MQ=255)
                                                        aCGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCa                                                                                                                                       <  2:584179/103‑1 (MQ=255)
                                                        aCGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTa                                                                                        >  1:684472/1‑150 (MQ=255)
                                                        aCGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTa                                                                                        <  2:31945/150‑1 (MQ=255)
                                                            aTCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGa                                                                                    >  2:1587050/1‑150 (MQ=255)
                                                             tCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGaa                                                                                   <  3:419232/150‑1 (MQ=255)
                                                                aCCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGt                                                                                >  1:1444272/1‑150 (MQ=255)
                                                                aCCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGt                                                                                >  3:550042/1‑150 (MQ=255)
                                                                 cccAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTg                                                                               <  1:1771446/150‑1 (MQ=255)
                                                                   cagcagGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTg                                                                               <  1:141341/148‑1 (MQ=255)
                                                                   cagcagGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTg                                                                               >  2:141341/1‑148 (MQ=255)
                                                                   cagcagGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCt                                                                             >  4:1093769/1‑150 (MQ=255)
                                                                       agGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTg                                                                         <  2:1444272/150‑1 (MQ=255)
                                                                           ggtATCCGTCGGTGCGTCCCTGCTCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATACGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGt                                                                     <  1:1012890/150‑1 (MQ=255)
                                                                           ggtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCAt                                                                                                                                                >  2:676617/1‑75 (MQ=255)
                                                                           ggtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCAt                                                                                                                                                <  1:676617/75‑1 (MQ=255)
                                                                           ggtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGt                                                                     <  4:468677/150‑1 (MQ=255)
                                                                            gtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAGCGTGCTGTTGCCGtt                                                                    <  4:596849/150‑1 (MQ=255)
                                                                            gtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGtt                                                                    <  3:54102/150‑1 (MQ=255)
                                                                            gtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGtt                                                                    >  2:840650/1‑150 (MQ=255)
                                                                              aTCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGa                                                                  <  1:840650/150‑1 (MQ=255)
                                                                               tCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGAc                                                                 >  2:1568069/1‑150 (MQ=255)
                                                                                   tCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCg                                                             <  1:912765/150‑1 (MQ=255)
                                                                                      gTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAgg                                                                                                                                     >  1:132520/1‑75 (MQ=255)
                                                                                      gTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAgg                                                                                                                                     <  2:132520/75‑1 (MQ=255)
                                                                                      gTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGa                                                                                                                            <  2:511314/84‑1 (MQ=255)
                                                                                      gTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGa                                                                                                                            >  1:511314/1‑84 (MQ=255)
                                                                                       tGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggtt                                                                                                  <  1:276727/109‑1 (MQ=255)
                                                                                       tGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggtt                                                                                                  >  2:276727/1‑109 (MQ=255)
                                                                                        gCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGACCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTa                                                        >  4:778537/1‑150 (MQ=255)
                                                                                         cGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTaa                                                       <  1:188049/150‑1 (MQ=255)
                                                                                          gTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTgg                                                                                          <  1:89740/114‑1 (MQ=255)
                                                                                          gTCCCTGATCCCGTACCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTgg                                                                                          >  2:89740/1‑114 (MQ=255)
                                                                                             ccTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGa                                                                                                                            <  2:1405977/77‑1 (MQ=255)
                                                                                             ccTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGa                                                                                                                            >  1:1405977/1‑77 (MQ=255)
                                                                                                          ggtgAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACgt                                      <  3:1007224/147‑1 (MQ=255)
                                                                                                           tGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACgtg                                     >  1:341786/1‑150 (MQ=255)
                                                                                                           tGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACgtg                                     <  3:544815/150‑1 (MQ=255)
                                                                                                             gAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACgtggt                                   <  2:1222116/150‑1 (MQ=255)
                                                                                                             gAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACgtggt                                   >  1:1434496/1‑150 (MQ=255)
                                                                                                              aaCACGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACgtggtg                                  <  1:1479082/150‑1 (MQ=255)
                                                                                                                cacGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACgtggtggt                                <  3:285717/150‑1 (MQ=255)
                                                                                                                cacGATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACgtggtggt                                <  2:1775261/150‑1 (MQ=255)
                                                                                                                   gATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGgtgt                                                         >  3:412236/1‑122 (MQ=255)
                                                                                                                   gATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGgtgt                                                         <  4:412236/122‑1 (MQ=255)
                                                                                                                   gATGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTc                             >  4:846239/1‑150 (MQ=255)
                                                                                                                    aTGACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCg                            >  1:1023787/1‑150 (MQ=255)
                                                                                                                      gACGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGtt                          >  2:149479/1‑150 (MQ=255)
                                                                                                                        cGCCAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCa                        <  1:804457/150‑1 (MQ=255)
                                                                                                                          ccAACCGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGt                      >  2:444442/1‑150 (MQ=255)
                                                                                                                              ccGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggt                                                                                                   <  1:209178/69‑1 (MQ=255)
                                                                                                                              ccGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggt                                                                                                   >  2:209178/1‑69 (MQ=255)
                                                                                                                              ccGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGt                  >  1:1110586/1‑150 (MQ=255)
                                                                                                                               cGTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTAt                                                                                       >  4:86119/1‑80 (MQ=255)
                                                                                                                                gTGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCg                                                                      >  3:81015/1‑96 (MQ=255)
                                                                                                                                 tGCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGa               >  4:791821/1‑150 (MQ=255)
                                                                                                                                  gCATTGATGGGTGCGACCATGCAACGTCGGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGAt              >  2:182430/1‑150 (MQ=255)
                                                                                                                                  gCATTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGAt              <  1:149479/150‑1 (MQ=255)
                                                                                                                                    aTTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGcc                                                                                                         >  3:1046416/1‑57 (MQ=255)
                                                                                                                                    aTTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGtt                                                                          <  2:1350866/88‑1 (MQ=255)
                                                                                                                                    aTTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGtt                                                                          >  1:1350866/1‑88 (MQ=255)
                                                                                                                                    aTTGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGc            >  3:490632/1‑150 (MQ=255)
                                                                                                                                     ttGATGGGTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCt           >  3:298559/1‑150 (MQ=255)
                                                                                                                                            gTGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTa    >  1:452956/1‑150 (MQ=255)
                                                                                                                                             tGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGtt                          <  2:51598/127‑1 (MQ=255)
                                                                                                                                             tGCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGtt                          >  1:51598/1‑127 (MQ=255)
                                                                                                                                              gCGACCATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATc  >  1:306592/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATC  >  NC_000913_3_pae_pgi/4183149‑4183440

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: