Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,918,916 T→A 7.5% K314N (AAA→AAT) ‡ atpA ← F1 sector of membrane‑bound ATP synthase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,918,9160TA 7.5% 448.8 / 20.3 161K314N (AAA→AAT) ‡atpAF1 sector of membrane‑bound ATP synthase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (72/77);  new base A (7/5);  total (79/82)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01

Only 100 of 164 total aligned reads displayed.

TACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAGGGAGATCTGACGGTAAGCAACAGCCTGTTT  >  NC_000913/3918768‑3919065
                                                                                                                                                    |                                                                                                                                                     
tACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCtt                                                                                                                                                      <  1:1642908/150‑1 (MQ=255)
  ccGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGGAACGCACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGCACCGGTTTTCCCTTTCACTTCCCCTTTg                                                                                                                                                    >  2:1291817/1‑150 (MQ=255)
   cGGCGTTGAACAGGTTGGTTTCCAGGAAGTTCTGACCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTgg                                                                                                                                                   <  2:346387/150‑1 (MQ=255)
     gCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATTGGCAGTGCGGTCAGAGAACCGGTTATCCCTTTCACTTCACCTTTGGTg                                                                                                                                                 >  1:1269525/1‑150 (MQ=255)
     gCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTg                                                                                                                                                 <  2:463593/150‑1 (MQ=255)
        ttGAACAGGTTGTTTTCCAGTAAGAGATGACCATCGGTAATGGATATTACGTTGGTCGTAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGCCAGAGCGGGCAGAGAACCGGTTTTCCCTTTCACTTCACCAAAGGTGAAg                                                                                                                                              <  2:295738/150‑1 (MQ=255)
         tGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAgg                                                                                                                                             <  1:994649/150‑1 (MQ=255)
             cAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTc                                                                                                                                         >  2:1492758/1‑150 (MQ=255)
             cAGGGTGGTTTCCAGGAAGATCTGCCCCTCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTc                                                                                                                                         <  1:363213/150‑1 (MQ=255)
               ggttggttTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCaa                                                                                                                                       >  2:588485/1‑150 (MQ=255)
                 ttggttTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACg                                                                                                                                     <  1:1600132/150‑1 (MQ=255)
                 ttggttTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACg                                                                                                                                     >  2:637015/1‑150 (MQ=255)
                    gttTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTAt                                                                                                                                  >  2:200679/1‑150 (MQ=255)
                    gttTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCAAAGGTGAAgtgaaa                                                                                                                                         >  2:1510614/1‑138 (MQ=255)
                       tCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGCACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCg                                                                                                                               >  1:817432/1‑150 (MQ=255)
                       tCCAGGAAGATCTGACCATCGGTAATGCAGATTACGTTGGTCGTAACTAACGCCGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCg                                                                                                                               <  1:1391354/150‑1 (MQ=255)
                       tCCAGGAAGATCTGACCAGCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGGGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCg                                                                                                                               <  1:1075587/150‑1 (MQ=255)
                        ccAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCgg                                                                                                                              >  1:542809/1‑150 (MQ=255)
                         cAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGc                                                                                                                             <  2:1529646/150‑1 (MQ=255)
                         cAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGc                                                                                                                             >  1:647552/1‑150 (MQ=255)
                           ggAAGAGCTGCCCATCGGTAATGGAGATTACGTGGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGGGCGGGCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGt                                                                                                                           <  2:97633/150‑1 (MQ=255)
                            gAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGtt                                                                                                                          <  1:637015/150‑1 (MQ=255)
                               gATCTGACCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCAGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACAGGTTTTCCCTTTCACTTGACGTTTGGTGAAGGCTTCAACGTATTCGGCGTTAac                                                                                                                       >  1:1258984/1‑150 (MQ=255)
                               gATCTGACCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGGCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAac                                                                                                                       <  2:720198/150‑1 (MQ=255)
                                aTCTGACCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAaca                                                                                                                      <  1:1294264/150‑1 (MQ=255)
                                   tGACCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGt                                                                                                                   <  1:1292016/150‑1 (MQ=255)
                                     aCCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTgc                                                                                                                 <  1:1535697/150‑1 (MQ=255)
                                      ccATCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTgcc                                                                                                                >  1:1060522/1‑149 (MQ=255)
                                      ccATCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTgca                                                                                                                <  2:1197538/150‑1 (MQ=255)
                                           ggTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTAAACACGTGCAGCACg                                                                                                           >  1:577284/1‑150 (MQ=255)
                                              aaTGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTc                                                                                                        >  1:569146/1‑150 (MQ=255)
                                               aTGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTcc                                                                                                       <  1:324727/150‑1 (MQ=255)
                                               aTGGAGATTACGTTGGTCGGAACGAACGCAGAAACGGCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTcc                                                                                                       <  1:1283824/150‑1 (MQ=255)
                                                 ggAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCcag                                                                                                     <  1:159863/150‑1 (MQ=255)
                                                  gagaTTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCcagc                                                                                                    <  2:1181608/150‑1 (MQ=255)
                                                  ctgaTTACGTTGGTCGGAACGAACGCAGAAAGGTCACCCGGATGAGTTTCGATAATCGGCAGTGCGGGCATAGATCCGGTTTTCCCTTTCACTTCACCTTTGTTGAATGCTTAAACGTATTGGGCGTTTACACGTGCAGCACGCTCcagc                                                                                                    <  1:914442/148‑1 (MQ=255)
                                                    gaTTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTACCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCACcagcaa                                                                                                  >  2:914442/1‑149 (MQ=255)
                                                         cGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACgag                                                                                             <  1:527120/150‑1 (MQ=255)
                                                         cGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACgag                                                                                             <  1:1363674/150‑1 (MQ=255)
                                                          gTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCca                                                                                                      >  2:426906/1‑140 (MQ=255)
                                                          gTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCca                                                                                                      <  1:426906/140‑1 (MQ=255)
                                                             ggTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTg                                                                                         >  1:192628/1‑150 (MQ=255)
                                                             ggTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTg                                                                                         <  2:1145284/150‑1 (MQ=255)
                                                               tCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGa                                                                                       >  1:1192226/1‑150 (MQ=255)
                                                               gcGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGa                                                                                       <  1:200679/149‑1 (MQ=255)
                                                                gggAACGAACGCAGAAACGTCACCCGACTGTGTGTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACATCACCTTTGGGGAAGGCTTCAACGTATTCGGCGTTAACACGTGCATCACGCTCCAGCAGACGAGAGTGGAg                                                                                      <  1:310544/149‑1 (MQ=255)
                                                                   aacgaacgCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCTGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTa                                                                                   <  2:1540275/150‑1 (MQ=255)
                                                                     cgaacgCAGAAACGTCCCCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGa                                                                                 >  2:780411/1‑150 (MQ=255)
                                                                        acgCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGaaaa                                                                              >  1:1175172/1‑150 (MQ=255)
                                                                         cgCAGAAACGGCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAAc                                                                             <  2:647552/150‑1 (MQ=255)
                                                                          gccgAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACg                                                                            <  1:479609/147‑1 (MQ=255)
                                                                           cAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGt                                                                           <  2:30941/150‑1 (MQ=255)
                                                                            aGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTc                                                                          <  2:482954/150‑1 (MQ=255)
                                                                             gAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCg                                                                         <  1:958191/150‑1 (MQ=255)
                                                                                aCGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGccc                                                                      <  2:66702/150‑1 (MQ=255)
                                                                                aCGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGccc                                                                      <  2:723088/150‑1 (MQ=255)
                                                                                    cACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGa                                                                  >  1:897379/1‑150 (MQ=255)
                                                                                     aCCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGaa                                                                 <  2:938513/150‑1 (MQ=255)
                                                                                      cccGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAAt                                                                >  1:799566/1‑150 (MQ=255)
                                                                                      cccGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAAt                                                                <  1:1693232/150‑1 (MQ=255)
                                                                                      cccGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCATCAACGTATTCGGCGTTAACACGTGCAGCAAGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAAt                                                                >  1:901376/1‑150 (MQ=255)
                                                                                        cGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGc                                                              >  2:1355699/1‑150 (MQ=255)
                                                                                          ccTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTTGTGAAGGCTTCAACGTATTCGGCTTTAACACTTTCAGCACTCTCCAGCAGACGAGAGTGGATGTAGAAAACGTCGCCCGGGAATGCtt                                                            <  1:1067202/150‑1 (MQ=255)
                                                                                           cTGAGGTTCGATAATCGGCAGGGCGGTCAGCGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTc                                                           <  1:917442/150‑1 (MQ=255)
                                                                                            tGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCATCAACGTATTCGGCGTTTGCACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTcc                                                          >  1:1511689/1‑149 (MQ=255)
                                                                                              aGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGTAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACg                                                        >  1:215421/1‑150 (MQ=255)
                                                                                               gTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGt                                                       >  2:356606/1‑150 (MQ=255)
                                                                                                tttCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTc                                                      <  1:1444578/150‑1 (MQ=255)
                                                                                                      tAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCg                                                <  2:1655411/150‑1 (MQ=255)
                                                                                                           ggCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGgacga                                           >  1:1520912/1‑150 (MQ=255)
                                                                                                                       agagAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACTTCTCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAgg                               <  1:86973/150‑1 (MQ=255)
                                                                                                                       agagAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAgg                               >  1:114549/1‑150 (MQ=255)
                                                                                                                        gagaACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAggg                              <  1:588485/150‑1 (MQ=255)
                                                                                                                        gagaACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAggg                              >  2:187759/1‑150 (MQ=255)
                                                                                                                           aaCCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAGGgaga                           <  2:1258984/150‑1 (MQ=255)
                                                                                                                             ccGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAGGGAGATc                         <  1:1492758/150‑1 (MQ=255)
                                                                                                                             ccGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAGGGAGATc                         >  2:568458/1‑150 (MQ=255)
                                                                                                                               ggTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAGGGAGATCTg                       <  2:1060571/150‑1 (MQ=255)
                                                                                                                                    tCCCTTTCACTTCACCTTTGGTTAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGATCAGCAGGGAGATCTGACGGt                  <  1:1026106/150‑1 (MQ=255)
                                                                                                                                    tCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAGGGAGATCTGACGGt                  >  2:1435011/1‑150 (MQ=255)
                                                                                                                                     cccTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAGGGAGATCTGACGGTa                 >  2:1207794/1‑150 (MQ=255)
                                                                                                                                     cccTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAGGGAGATCTGACGGTa                 >  1:756533/1‑150 (MQ=255)
                                                                                                                                      ccTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAGGGAGATCTGACGGTaa                >  1:702405/1‑150 (MQ=255)
                                                                                                                                       cTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAGGGAGATCTGACGGTAAg               <  2:897379/150‑1 (MQ=255)
                                                                                                                                        attcacatcacCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGAAGGGAGATCTGACGGTAAGc              <  2:1055785/149‑1 (MQ=255)
                                                                                                                                         aagccttcacCAAAGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAGGGAGATCTGACGGTAAGCa             >  2:1252562/5‑150 (MQ=255)
                                                                                                                                         aagccttcacCAAAGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAGGGAGATCTGACGGTAAGCa             <  1:1180630/146‑1 (MQ=255)
                                                                                                                                         aagccttcacCAAAGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAGGGAGATCTGACGGTAAGCa             >  2:589266/5‑150 (MQ=255)
                                                                                                                                           gccttcacCAAAGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAGGGAGATCTGACGGTAAGCAAc           >  1:1580358/3‑150 (MQ=255)
                                                                                                                                            ccttcacCAAAGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAGGGAGATCTGACGGTAAGCAAc           <  1:846666/148‑1 (MQ=255)
                                                                                                                                            ccttcacCAAAGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAGGGAGATCTGACGGTAAGCAAc           >  2:846666/2‑149 (MQ=255)
                                                                                                                                            acttcacCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAGGGAGATCTGACGGTAAGCAACa          >  1:1705417/1‑150 (MQ=255)
                                                                                                                                             cttcacCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAGGGAGATCTGACGGTAAGCACCAg         >  2:846240/1‑150 (MQ=255)
                                                                                                                                               tcacCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAGGGAGATCTGACGGTAAGCAACAGcc       >  1:76215/1‑150 (MQ=255)
                                                                                                                                               tcacCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAGGGAGATCTGACGGTAAGCAACAGcc       <  2:1387600/150‑1 (MQ=255)
                                                                                                                                               tcacCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAGGGAGATCTGACGGTAAGCAACAGcc       >  2:367234/1‑150 (MQ=255)
                                                                                                                                               tcacCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAGGGAGATCTGACGGTAAGCAACAGcc       <  2:1598282/150‑1 (MQ=255)
                                                                                                                                                cacCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAGGGAGATCTGACGGTAAGCAACAGCCt      >  2:1098244/1‑150 (MQ=255)
                                                                                                                                                  ccTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAGGGAGATCTGACGGTAAGCAACAGCCTGt    >  1:1498538/1‑150 (MQ=255)
                                                                                                                                                    tttGGTGAAGGCTTCAACGTTTTCGGCGTTAACACGTGCAGCACGCTCCAGCAGTCGAGAGTGGAGGTAGAAAACGTCGGCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAGGGAGATCTGACGGTAAGCAACAGCCTGttt  <  1:333399/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGCAGAAACGTCACCCGCCTGAGTTTCGATAATCGGCAGTGCGGTCAGAGAACCGGTTTTCCCTTTCACTTCACCTTTGGTGAAGGCTTCAACGTATTCGGCGTTAACACGTGCAGCACGCTCCAGCAGACGAGAGTGGAGGTAGAAAACGTCGCCCGGGAATGCTTCACGTCCTGGCGGACGACGGAGCAGCAGGGAGATCTGACGGTAAGCAACAGCCTGTTT  >  NC_000913/3918768‑3919065

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