Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,207,576 A→G 5.6% intergenic (‑44/‑571) purH ← / → rrsE IMP cyclohydrolase and phosphoribosylaminoimidazolecarboxamide formyltransferase/16S ribosomal RNA of rrnE operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,207,5760AG 5.6% 1472.5 / 49.6 451intergenic (‑44/‑571)purH/rrsEIMP cyclohydrolase and phosphoribosylaminoimidazolecarboxamide formyltransferase/16S ribosomal RNA of rrnE operon
Reads supporting (aligned to +/- strand):  ref base A (206/212);  new base G (11/14);  total (217/234)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.54e-01

Only 100 of 456 total aligned reads displayed.

ACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTGGGCAAAATACAG  >  NC_000913/4207427‑4207725
                                                                                                                                                     |                                                                                                                                                     
   gcagcTCCACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAAttt                                                                                                                                                    <  1:5695564/150‑1 (MQ=255)
    cagcTCCACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAAtttt                                                                                                                                                   >  2:4670057/1‑150 (MQ=255)
     agcTCCACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGTTTTGACTATTACAGAGAGCGTTAGCTGAAttttt                                                                                                                                                  >  1:6382847/1‑150 (MQ=255)
      gcTCCACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTc                                                                                                                                                 <  1:3887575/150‑1 (MQ=255)
       cTCCACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTcg                                                                                                                                                >  2:3844569/1‑150 (MQ=255)
         ccACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTcgcg                                                                                                                                              >  2:1267674/1‑150 (MQ=255)
         ccACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGAGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTcgcg                                                                                                                                              >  2:2561490/1‑150 (MQ=255)
          cacaCCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGa                                                                                                                                             <  2:3533437/150‑1 (MQ=255)
          cacaCCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGa                                                                                                                                             <  2:3329053/150‑1 (MQ=255)
           acacCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGaa                                                                                                                                            >  1:5623823/1‑150 (MQ=255)
              ccGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGaaaaa                                                                                                                                         <  2:6203750/150‑1 (MQ=255)
                gcgTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACt                                                                                                                                       <  1:3844569/150‑1 (MQ=255)
                gcgTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACt                                                                                                                                       <  1:160411/150‑1 (MQ=255)
                  gTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAGTTTTTCGCGAAAAATTCa                                                                                                                                     <  2:3561982/150‑1 (MQ=255)
                  gTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAGTTTTTCGCGAAAAATTCa                                                                                                                                     <  2:4768751/150‑1 (MQ=255)
                          aGTGCGTGGGCGATTTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCTTGGTAAATCCCCTGGTTTTGACTATTACAGAGTGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGc                                                                                                                             <  1:883611/150‑1 (MQ=255)
                          aGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGc                                                                                                                             <  2:38966/150‑1 (MQ=255)
                           gTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGcc                                                                                                                            <  2:5054143/150‑1 (MQ=255)
                            tGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGccc                                                                                                                           >  2:2409487/1‑150 (MQ=255)
                            tGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGccc                                                                                                                           >  2:3772463/1‑150 (MQ=255)
                               cTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTa                                                                                                                        >  1:469435/1‑150 (MQ=255)
                                   gCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACgg                                                                                                                    >  2:5359178/1‑150 (MQ=255)
                                   gCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACgg                                                                                                                    >  1:1253672/1‑150 (MQ=255)
                                    cGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAAttttt                                                                                                                                                  <  2:5841946/119‑1 (MQ=255)
                                    cGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACggg                                                                                                                   <  2:4831487/150‑1 (MQ=255)
                                    cGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACggg                                                                                                                   <  2:172571/150‑1 (MQ=255)
                                       aTTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCa                                                                                                                <  2:2473347/150‑1 (MQ=255)
                                       aTTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCa                                                                                                                >  2:558850/1‑150 (MQ=255)
                                         tCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATc                                                                                                              >  1:1794850/1‑150 (MQ=255)
                                          cgacgaTACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATcc                                                                                                             <  2:1959608/150‑1 (MQ=255)
                                           gacgaTACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCt                                                                                                            >  2:6370890/1‑150 (MQ=255)
                                            acgaTACCGGCTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTa                                                                                                          <  2:1526719/150‑1 (MQ=255)
                                                 aCCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTAtttt                                                                                                      >  1:708995/1‑150 (MQ=255)
                                                   cGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTc                                                                                                    <  2:2758766/150‑1 (MQ=255)
                                                      cTTTGTCAGAAACCCTGCGCCGCGCGGGGGGGACTGGCCGAGGTTGTTGCTTGGTAAATCCCCTGGTTTTTACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGcc                                                                                                 <  1:2942119/150‑1 (MQ=255)
                                                      cTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTAGCGAAAAACTCAGCTAACGCCCCTAACCGGGCATCCATATTATTCgac                                                                                                 >  1:3702889/1‑148 (MQ=255)
                                                          gTCAGAAACACTGCGCCGCGCGCGGCGGACTGGACGACGTTGTTGCAGGTTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACACAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCa                                                                                             <  1:1123669/150‑1 (MQ=255)
                                                           tCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAGTTTTTCGCGAAAAATTCAGCTa                                                                                                                                 <  2:3188461/113‑1 (MQ=255)
                                                           tCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACgg                                                                                                                    >  1:6607486/1‑126 (MQ=255)
                                                           tCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACgg                                                                                                                    <  2:6607486/126‑1 (MQ=255)
                                                           tCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCTCCCGCAt                                                                                            >  1:4312687/1‑150 (MQ=255)
                                                           tCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCAt                                                                                            >  1:4619025/1‑150 (MQ=255)
                                                           tCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCAt                                                                                            >  1:907180/1‑150 (MQ=255)
                                                            cAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCAtt                                                                                           <  2:4452228/150‑1 (MQ=255)
                                                             aGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTg                                                                                          >  1:1011310/1‑150 (MQ=255)
                                                              gAACCACTGAGCCGCGCGCGGCGGACTGGACGACGTTGTTGCATGTTAAATCCCCTGTATTTGACTATTACAGAGAGCGTTAGCTTAATTTTTCGCTAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGt                                                                                         <  1:2409102/150‑1 (MQ=255)
                                                               aaaCACTGAGCTGAGCGCGGCGGACTGGACGAGGTTGTTTCTTGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTTAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTa                                                                                        <  1:158199/150‑1 (MQ=255)
                                                               aaaCACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTa                                                                                        >  2:1766980/1‑150 (MQ=255)
                                                                      gAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGaaaa                                                                                 <  2:1418696/150‑1 (MQ=255)
                                                                              cgcgGCGGTCTGGACGTCGTTGTTGCCTGTTAATTCCCGTGGTTTTGACTTTTACTGTGTGTGTTTTCTTAATTTTTCGCGAAAAACACAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTgcg                                                                         <  1:6208530/150‑1 (MQ=255)
                                                                                cggcggACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCa                                                                       >  2:5506986/1‑150 (MQ=255)
                                                                                 ggcggACTGGACGTGGTTGTTTCATTGTAAATCCCCTGGTTTTTACTTTTACAGAGAGCGTTAGCTGAATTTTTCGCTAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCaa                                                                      <  1:6163676/150‑1 (MQ=255)
                                                                                        tGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCAAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAAAGTTTGCGCAACGCTcgc                                                               >  1:1403833/1‑150 (MQ=255)
                                                                                         ggACGACGTTTTTGCATGGTAAATCCCCTGGTTTTGACTATTACAGAGAGCTTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTcgcg                                                              <  1:3997067/150‑1 (MQ=255)
                                                                                          gacgacGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGa                                                             >  2:3359749/1‑150 (MQ=255)
                                                                                          gacgacGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGAAACGCTCGCGa                                                             >  2:4201374/1‑150 (MQ=255)
                                                                                                gttgttGCATGGTAAATCCCCTGGATTTGACTATTACAGAGTGCGTTAGCTTAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAAttttt                                                       <  1:4928921/150‑1 (MQ=255)
                                                                                                gttgttGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAACAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAAttttt                                                       <  1:5066889/150‑1 (MQ=255)
                                                                                                    ttgCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTctct                                                   >  2:6559419/1‑150 (MQ=255)
                                                                                                    ttgCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTctct                                                   >  2:6012032/1‑150 (MQ=255)
                                                                                                      gCATTGTAAATCCCCTGGATTTGAATATTACAGAGAGCGTTAGCTTAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCttt                                                 <  1:759044/150‑1 (MQ=255)
                                                                                                      gCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTTAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCttt                                                 <  1:3456031/150‑1 (MQ=255)
                                                                                                       cATGGTAAATCCCCTGGATTTTACTATTACAGAGAGCGTTATCTGAATTTATCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTc                                                <  1:3355946/150‑1 (MQ=255)
                                                                                                       cATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCAGTGTAACGAAAACGTGTGCGCAACGCTCGCGAATTTTTCTCTTTc                                                >  2:3128133/1‑150 (MQ=255)
                                                                                                        aTGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCa                                               <  1:3168379/150‑1 (MQ=255)
                                                                                                        aTGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCa                                               >  1:5337582/1‑150 (MQ=255)
                                                                                                         tGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCa                                                                                                                <  2:2117772/84‑1 (MQ=255)
                                                                                                         tGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCaa                                              >  2:2827287/1‑150 (MQ=255)
                                                                                                         tGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCaa                                              >  2:2032507/1‑150 (MQ=255)
                                                                                                         tGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACACAGCTAACGCCCCTAACGGGGCa                                                                                                                >  1:2117772/1‑84 (MQ=255)
                                                                                                           ttAAATCCCCTGGATTTTACTATTACAGAGAGCGTTTGCTTAATTTTTCGCGAAAAACACAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATg                                            <  1:2807710/149‑1 (MQ=255)
                                                                                                           gTAAATCCACTGGATTTGAATATTACAGAGAGCGTTAGCTTAATTTTTCTCTAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATg                                            <  1:2866555/150‑1 (MQ=255)
                                                                                                                tCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTAAATTTTCGCCCGCCTGGTAACGCAAACGGTGGCGCACCGCTCGCGAATGTTTCTCTTGCAATGGTGag                                       >  1:5889056/1‑149 (MQ=255)
                                                                                                                   ccTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATcac                                    >  2:4956363/1‑150 (MQ=255)
                                                                                                                     tGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCTAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACaa                                  <  1:4496901/150‑1 (MQ=255)
                                                                                                                     tGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATGCTTATTTTTCGCCCGCATTGTAACGAATACGTTTGGGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACaa                                  >  1:1714998/1‑150 (MQ=255)
                                                                                                                       gATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAAtt                                >  2:5100911/1‑150 (MQ=255)
                                                                                                                       gATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTAGCCCGCATTGTAAAGAAAACGTTTGCGCAACGCACGCGAATTTTTCTCTTTCAATGGTGATCACAAtt                                >  1:4540530/1‑150 (MQ=255)
                                                                                                                        aTTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAAttt                               >  2:2990288/1‑150 (MQ=255)
                                                                                                                           tGACTATTACAGAGCGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGa                            <  1:4206377/150‑1 (MQ=255)
                                                                                                                            gACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTAGCCCGCAGTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGAc                           >  1:3115881/1‑150 (MQ=255)
                                                                                                                             aCTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACt                          >  2:4107642/1‑150 (MQ=255)
                                                                                                                             aCTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACt                          >  2:5146918/1‑150 (MQ=255)
                                                                                                                             aCTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCAAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTTACt                          >  1:4872418/1‑150 (MQ=255)
                                                                                                                               tATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCGATGGTGATCACAATTTTGACtgt                        >  2:4603071/1‑150 (MQ=255)
                                                                                                                               tATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCAAACGGGGCATCCTTATTTTTAGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACtgt                        >  1:3249440/1‑150 (MQ=255)
                                                                                                                               tATTACAGAGAGCGTTAGCTGAATTTTTAGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTAGCCCGCATTGTAACGTAAAGGTTTGCGCAACGCTCGCGAATTTTTGTCATTCAATTGTGATGACAATTTTGAgtgt                        >  1:6059282/1‑150 (MQ=255)
                                                                                                                                aTTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACtgtg                       <  2:5658783/150‑1 (MQ=255)
                                                                                                                                  tACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGt                     >  2:57143/1‑150 (MQ=255)
                                                                                                                                   aCAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGtt                    <  1:5417062/150‑1 (MQ=255)
                                                                                                                                     agagagCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTAc                  <  2:2388644/150‑1 (MQ=255)
                                                                                                                                       agagCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCg                >  2:710174/1‑150 (MQ=255)
                                                                                                                                       agagCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCg                <  2:1011325/150‑1 (MQ=255)
                                                                                                                                           cGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTggg            >  2:1058682/1‑150 (MQ=255)
                                                                                                                                               aGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTGGGCaaa        <  2:2277574/150‑1 (MQ=255)
                                                                                                                                                gCTTAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTGGGCaaaa       <  1:582855/150‑1 (MQ=255)
                                                                                                                                                gCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTGGGCaaaa       <  2:4300243/150‑1 (MQ=255)
                                                                                                                                                  tGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTGGGCAAAATa     >  2:5241784/1‑150 (MQ=255)
                                                                                                                                                    aaTTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTGGGCAAAATACa   <  1:3553257/150‑1 (MQ=255)
                                                                                                                                                     aTTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTGGGCAAAATACAg  <  2:1842495/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCACAATTTTGACTGTGGTTACCGTGGGCAAAATACAG  >  NC_000913/4207427‑4207725

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