Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,038,390 A→G 6.4% L217L (TTA→TTG appC → cytochrome bd‑II ubiquinol oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,038,3900AG 6.4% 285.4 / 7.4 93L217L (TTA→TTGappCcytochrome bd‑II ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base A (37/50);  new base G (2/4);  total (39/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

ATGACCAGCTTCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACC  >  NC_000913/1038241‑1038536
                                                                                                                                                     |                                                                                                                                                  
gtGACCAGCTTCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTa                                                                                                                                                    <  1:125318/149‑1 (MQ=255)
    ccAGCTTCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCt                                                                                                                                                <  1:938101/150‑1 (MQ=255)
      aGCTTCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGTAATGCGACATTGCGCTcc                                                                                                                                              <  2:13382/150‑2 (MQ=255)
      aGCTTCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGTAATGCGACATTGCGCTcc                                                                                                                                              <  2:40358/150‑2 (MQ=255)
        cTTCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGtt                                                                                                                                            <  2:456242/150‑1 (MQ=255)
         ttCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGttt                                                                                                                                           <  2:224810/150‑1 (MQ=255)
           cAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGc                                                                                                                                         <  2:1342050/150‑1 (MQ=255)
           cAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGc                                                                                                                                         <  2:1355005/150‑1 (MQ=255)
            aGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGcc                                                                                                                                        >  2:868755/1‑150 (MQ=255)
             gCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCa                                                                                                                                       >  1:1246597/1‑150 (MQ=255)
              cGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCAt                                                                                                                                      >  2:821202/1‑150 (MQ=255)
               gAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATc                                                                                                                                     <  1:1225928/150‑1 (MQ=255)
                aGCTGGTCTTTAATCCTGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCTGTGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGTAGCGCAATGTCGCATTACGCTCGTTTGCCATCg                                                                                                                                    >  2:1640823/1‑150 (MQ=255)
                aGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCg                                                                                                                                    <  2:792103/150‑1 (MQ=255)
                 gCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCgg                                                                                                                                   >  2:695181/1‑150 (MQ=255)
                    ggTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTc                                                                                                                                <  1:1295918/150‑1 (MQ=255)
                    ggTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTc                                                                                                                                <  1:684368/150‑1 (MQ=255)
                    ggTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTc                                                                                                                                <  2:1593704/150‑1 (MQ=255)
                    ggTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTc                                                                                                                                <  2:1020467/150‑1 (MQ=255)
                    ggTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTc                                                                                                                                >  2:1457188/1‑150 (MQ=255)
                      tCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCg                                                                                                                              <  2:1641520/150‑1 (MQ=255)
                       cTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGt                                                                                                                             >  2:1206065/1‑150 (MQ=255)
                       cTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGt                                                                                                                             >  1:159521/1‑150 (MQ=255)
                        tttAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTc                                                                                                                            >  2:1252704/1‑150 (MQ=255)
                            aTCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTTTGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCg                                                                                                                        >  1:993782/1‑150 (MQ=255)
                             tCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCgg                                                                                                                       >  2:1395384/1‑150 (MQ=255)
                              ccGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGt                                                                                                                      <  1:113160/150‑1 (MQ=255)
                                ggTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTg                                                                                                                                          <  1:655769/128‑1 (MQ=255)
                                ggTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTg                                                                                                                                          >  2:655769/1‑128 (MQ=255)
                                  tCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTActc                                                                                                                  <  1:665888/150‑1 (MQ=255)
                                  tCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTActc                                                                                                                  <  1:868755/150‑1 (MQ=255)
                                   cAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTActct                                                                                                                 <  1:572015/150‑1 (MQ=255)
                                     gCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTgg                                                                                                               >  1:792758/1‑150 (MQ=255)
                                         ggTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCGTCGGTTCCGTCTTCGGTACTCTGGCGAt                                                                                                           >  2:599365/1‑150 (MQ=255)
                                          gTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGAtt                                                                                                          >  2:1368548/1‑150 (MQ=255)
                                           tGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTa                                                                                                         >  2:422114/1‑150 (MQ=255)
                                           tGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTa                                                                                                         >  1:1078548/1‑150 (MQ=255)
                                              aaTTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCg                                                                                                      <  1:1589361/150‑1 (MQ=255)
                                                 ttGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTa                                                                                                   <  1:1206065/150‑1 (MQ=255)
                                                   gtgCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTAcc                                                                                                 <  2:1084196/150‑1 (MQ=255)
                                                    tgCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTTTGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTAccc                                                                                                <  2:993782/150‑1 (MQ=255)
                                                    tgCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTAccc                                                                                                >  1:40318/1‑150 (MQ=255)
                                                     gCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCt                                                                                               >  1:993220/1‑150 (MQ=255)
                                                        cacCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCa                                                                                            <  2:1310314/150‑1 (MQ=255)
                                                          ccGTAATGGCGGGCTACGTTACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTTCAAc                                                                                          <  2:157392/150‑1 (MQ=255)
                                                          ccGTAATGGCGGGCTACGTGACCGGGGCAATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAAc                                                                                          >  2:748492/1‑150 (MQ=255)
                                                              aaTGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCgg                                                                                      <  2:1364385/150‑1 (MQ=255)
                                                               aTGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCAATTATCGGTACCCTGCAACTCGga                                                                                     <  1:528046/150‑1 (MQ=255)
                                                                  gCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGAc                                                                                  <  1:695181/150‑1 (MQ=255)
                                                                   cGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACa                                                                                 >  1:427810/1‑150 (MQ=255)
                                                                        tACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCt                                                                            <  1:266732/150‑1 (MQ=255)
                                                                            ggACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAATTCGGAGACAGTTCTGCGt                                                                        <  2:1554973/149‑1 (MQ=255)
                                                                             gACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTa                                                                       >  2:1074083/1‑150 (MQ=255)
                                                                             gACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTa                                                                       <  1:315832/150‑1 (MQ=255)
                                                                              cccGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTAt                                                                      <  1:203160/149‑1 (MQ=255)
                                                                              aCCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTAt                                                                      >  1:311882/1‑150 (MQ=255)
                                                                               ccttggCCATGTTTATTATGGCGATCAGCGCCTTGTATTTACTTCGCGTACGGTAGCGCAATGTCTCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATg                                                                     <  1:634537/146‑1 (MQ=255)
                                                                               ccGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATg                                                                     >  1:763351/1‑150 (MQ=255)
                                                                                cGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGa                                                                    <  2:159521/150‑1 (MQ=255)
                                                                                 ggggCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGaa                                                                   >  1:686314/1‑150 (MQ=255)
                                                                                  gggCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAg                                                                  <  1:1046817/150‑1 (MQ=255)
                                                                                   ggCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGt                                                                 <  1:1608060/150‑1 (MQ=255)
                                                                                       aTGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTcgcg                                                             <  1:1240839/150‑1 (MQ=255)
                                                                                        tGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCgcgc                                                            >  1:1125935/1‑150 (MQ=255)
                                                                                           ttattaTGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAg                                                         >  2:1184293/1‑150 (MQ=255)
                                                                                           ttattaTGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGt                                                        <  2:152599/150‑1 (MQ=255)
                                                                                            tattatGGCGATCAGCGCCTGGTATTTACTGAGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGt                                                        >  2:430849/1‑150 (MQ=255)
                                                                                              ttatGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTAc                                                      >  1:982495/1‑150 (MQ=255)
                                                                                                 tGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAAc                                                   <  1:1074083/150‑1 (MQ=255)
                                                                                                 tGGCGATCAGAGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAAc                                                   >  2:278926/1‑150 (MQ=255)
                                                                                                   gCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCg                                                 >  2:1405325/1‑150 (MQ=255)
                                                                                                   gCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCg                                                 >  1:656888/1‑150 (MQ=255)
                                                                                                    cGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCgg                                                >  2:1046762/1‑150 (MQ=255)
                                                                                                              cTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTggc                                      >  2:286592/1‑150 (MQ=255)
                                                                                                                  tATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGgcggcg                                  <  2:656888/150‑1 (MQ=255)
                                                                                                                       aCTGAGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGa                             <  1:430849/150‑1 (MQ=255)
                                                                                                                         tGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAg                           >  2:857334/1‑150 (MQ=255)
                                                                                                                         tGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAg                           <  1:794219/150‑1 (MQ=255)
                                                                                                                               gACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGt                     <  2:1311719/150‑1 (MQ=255)
                                                                                                                                aCGGTAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTg                    <  2:263222/150‑1 (MQ=255)
                                                                                                                                 cGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTgg                   <  2:993220/150‑1 (MQ=255)
                                                                                                                                 aacgAGCGTAATGCGACATTGCGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATg                               <  2:1334874/135‑1 (MQ=255)
                                                                                                                                 aacgAGCGTAATGCGACATTGCGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATg                               >  1:1334874/4‑138 (MQ=255)
                                                                                                                                   ggAGCGCAATGTCGCATTACGCTCGTTTGCCGTCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCa                 <  1:599365/150‑1 (MQ=255)
                                                                                                                                   cgAGCGTAATGCGACATTGCGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCa                 >  1:376302/2‑150 (MQ=255)
                                                                                                                                      gcgcAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACGGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAAc              <  1:977270/150‑1 (MQ=255)
                                                                                                                                        gcAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTAGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACgg            >  1:786324/1‑150 (MQ=255)
                                                                                                                                         aaaTGTCGCATTACGCTCGTTTTCCATCGTTTCCGTCTTCGGTACTCTTGCTATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGTAGTCGCGCAAGTACAACCGGTAAATCTGGCGGAGATGGAAGGGGAGTGGCAAACGGa           <  2:1246597/149‑1 (MQ=255)
                                                                                                                                            tGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAacc        <  2:792758/150‑1 (MQ=255)
                                                                                                                                            tGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAacc        <  2:1125935/150‑1 (MQ=255)
                                                                                                                                            tGTCGCATTACGCTCGTTTGCCATCGGTTCAGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAacc        <  1:1469199/150‑1 (MQ=255)
                                                                                                                                              tCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAacctg      <  1:748492/150‑1 (MQ=255)
                                                                                                                                              tCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAacctg      >  1:525461/1‑150 (MQ=255)
                                                                                                                                                  aTTGCGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAacctgcacc  <  2:376302/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ATGACCAGCTTCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAGTACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACC  >  NC_000913/1038241‑1038536

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