Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,347,175 A→T 5.8% intergenic (+25/‑209) nrdA → / → nrdB ribonucleoside‑diphosphate reductase 1 subunit alpha/ribonucleoside‑diphosphate reductase 1 subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,347,1750AT 5.8% 411.0 / 9.5 127intergenic (+25/‑209)nrdA/nrdBribonucleoside‑diphosphate reductase 1 subunit alpha/ribonucleoside‑diphosphate reductase 1 subunit beta
Reads supporting (aligned to +/- strand):  ref base A (59/59);  new base T (4/3);  total (64/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 128 total aligned reads displayed.

TACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATC  >  NC_000913/2347028‑2347323
                                                                                                                                                   |                                                                                                                                                    
tACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAAc                                                                                                                                                    >  1:483464/1‑150 (MQ=255)
 aCAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACg                                                                                                                                                   <  2:577321/150‑1 (MQ=255)
        cGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCTGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATc                                                                                                                                            >  2:1195750/1‑150 (MQ=255)
           ggTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                         <  1:1048396/150‑1 (MQ=255)
           ggTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                         <  2:1164561/150‑1 (MQ=255)
             tCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATAAGGCGTTTACGCCGCATCCGGca                                                                                                                                       <  1:87373/150‑2 (MQ=255)
               aaaaCACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                                                                                     >  2:277717/1‑150 (MQ=255)
                aaaCACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                                                                                    <  1:1051293/150‑1 (MQ=255)
                aaaCACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                                                                                    >  2:631037/1‑150 (MQ=255)
                aaaCACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                                                                                    >  1:460579/1‑150 (MQ=255)
                aaaCACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATAAGGCGTTTACGCCGCATCCGGcatct                                                                                                                                    <  1:1374947/150‑5 (MQ=255)
                    acTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCt                                                                                                                                >  2:903578/1‑150 (MQ=255)
                     cTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTc                                                                                                                               >  2:82828/1‑150 (MQ=255)
                       gTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCgg                                                                                                                             >  1:1285141/1‑150 (MQ=255)
                       gTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCgg                                                                                                                             >  1:781891/1‑150 (MQ=255)
                        tattatCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGt                                                                                                                            <  2:232782/150‑1 (MQ=255)
                         attatCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGtt                                                                                                                           <  1:113517/150‑1 (MQ=255)
                           tatCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATAAGGCGTTTACGCCGCATCCGGcatc                                                                                                                                     <  1:760748/138‑4 (MQ=38)
                           tatCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATAAGGCGTTTACGCCGCATCCGGcatc                                                                                                                                     >  2:760748/1‑135 (MQ=38)
                            atCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGcc                                                                                                                                                 >  2:1513986/1‑125 (MQ=255)
                            atCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGcc                                                                                                                                                 <  1:1513986/125‑1 (MQ=255)
                            atCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                                                                                      >  2:139724/1‑136 (MQ=255)
                            atCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                                                                                      <  1:139724/136‑1 (MQ=255)
                             tCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTa                                                                                                                       >  1:477065/1‑150 (MQ=255)
                                 aaCACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGcc                                                                                                                   <  1:1376842/150‑1 (MQ=255)
                                  acacCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAGGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCt                                                                                                                  <  1:575632/150‑1 (MQ=255)
                                   cacCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTg                                                                                                                 <  1:403453/150‑1 (MQ=255)
                                   cacCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTg                                                                                                                 <  1:971665/150‑1 (MQ=255)
                                   cacCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTg                                                                                                                 <  2:460579/150‑1 (MQ=255)
                                      ccGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATa                                                                                                              <  1:957949/150‑1 (MQ=255)
                                      ccGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATa                                                                                                              <  2:1233404/150‑1 (MQ=255)
                                      ccGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATa                                                                                                              <  1:1163512/150‑1 (MQ=255)
                                        gTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAg                                                                                                            <  2:1182725/150‑1 (MQ=255)
                                          gACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGAc                                                                                                          <  1:78948/150‑1 (MQ=255)
                                              gcgcTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACgcgc                                                                                                      >  2:1206938/1‑150 (MQ=255)
                                              gcgcTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATAAGGCGTTCACGCTGCATCCGGcacc                                                                                                                                     >  2:1223207/1‑116 (MQ=37)
                                              gcgcTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATAAGGCGTTCACGCTGCATCCGGcacc                                                                                                                                     <  1:1223207/119‑4 (MQ=37)
                                                gcTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCa                                                                                                    >  2:8683/1‑150 (MQ=255)
                                                gcTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGATTCTAGGCCTGATAAGACGCGCCa                                                                                                    >  2:627212/1‑150 (MQ=255)
                                                     aGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTc                                                                                               <  2:1223591/150‑1 (MQ=255)
                                                      gACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCg                                                                                              <  1:1315071/150‑1 (MQ=255)
                                                       aCGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGc                                                                                             <  1:1194115/150‑1 (MQ=255)
                                                         gCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCAt                                                                                           >  1:1153825/1‑150 (MQ=255)
                                                          cacaAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATc                                                                                          >  2:112244/1‑150 (MQ=255)
                                                          cacaAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATc                                                                                          >  2:1059589/1‑150 (MQ=255)
                                                          cacaAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATc                                                                                          >  1:1237041/1‑150 (MQ=255)
                                                          cacaAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATc                                                                                          >  2:712664/1‑150 (MQ=255)
                                                           acaAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCa                                                                                         <  2:1285141/150‑1 (MQ=255)
                                                            caAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAg                                                                                        >  2:453402/1‑150 (MQ=255)
                                                             aaGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAgg                                                                                       <  1:1468575/150‑1 (MQ=255)
                                                                 cGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTcc                                                                                   >  2:73517/1‑150 (MQ=255)
                                                                         tGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCg                                                                           <  2:70641/150‑1 (MQ=255)
                                                                           ccGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGa                                                                         <  1:852768/150‑1 (MQ=255)
                                                                              tCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGc                                                                      <  2:1294102/150‑1 (MQ=255)
                                                                               cAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCa                                                                     >  2:1327660/1‑150 (MQ=255)
                                                                               cAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCa                                                                     >  2:359613/1‑150 (MQ=255)
                                                                                aaTCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAg                                                                    <  1:607836/150‑1 (MQ=255)
                                                                                   ccAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGt                                                                 <  2:477065/150‑1 (MQ=255)
                                                                                    cAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTg                                                                >  2:769504/1‑150 (MQ=255)
                                                                                    cAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTg                                                                >  2:1433181/1‑150 (MQ=255)
                                                                                    cAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTg                                                                >  2:638922/1‑150 (MQ=255)
                                                                                    cAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTg                                                                <  2:937288/150‑1 (MQ=255)
                                                                                    cAGGACGATGGCAGCGAAAGAGGCGAATTTAAGATCTGATATTGACATGCCGGATTCGGCGAAAACGCCTTATCCGGCCTACGGCTCGGTATGTACGCCAGATAAGACGAGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTg                                                                >  2:945569/1‑150 (MQ=25)
                                                                                     aGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCCGGCACCCGGAGCCTGATGcgacgctg                                                               <  2:1487101/150‑8 (MQ=25)
                                                                                     aGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGa                                                               <  1:112244/150‑1 (MQ=255)
                                                                                      ggACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGaa                                                              >  1:1071009/1‑150 (MQ=255)
                                                                                          gATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGcc                                                          <  1:359613/150‑1 (MQ=255)
                                                                                            tGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCtt                                                        >  1:770589/1‑150 (MQ=255)
                                                                                            tGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCtt                                                        <  1:418812/150‑1 (MQ=255)
                                                                                              gCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTaa                                                      >  1:618362/1‑149 (MQ=255)
                                                                                                     aaGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCt                                               >  2:375932/1‑150 (MQ=255)
                                                                                                      aGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTa                                              <  1:277717/150‑1 (MQ=255)
                                                                                                           cgcATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCt                                         >  1:359932/1‑150 (MQ=255)
                                                                                                             cATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCg                                       <  1:8683/150‑1 (MQ=255)
                                                                                                               tGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGa                                     >  2:36226/1‑150 (MQ=255)
                                                                                                                gTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCt                                         >  1:53593/1‑145 (MQ=255)
                                                                                                                gTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCt                                         <  2:53593/145‑1 (MQ=255)
                                                                                                                gTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGAt                                    >  1:568291/1‑150 (MQ=255)
                                                                                                                   aGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTg                                 >  2:553061/1‑150 (MQ=255)
                                                                                                                    gATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGt                                >  1:1227101/1‑150 (MQ=255)
                                                                                                                    gATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGt                                <  2:681665/150‑1 (MQ=255)
                                                                                                                     aTCTGATATTGAGATGCCGGATGAGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTa                               >  1:75815/1‑150 (MQ=255)
                                                                                                                        tGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGc                            >  1:58117/1‑150 (MQ=255)
                                                                                                                           tatTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTg                         >  1:371797/1‑150 (MQ=255)
                                                                                                                             ttGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGAt                       <  2:329806/150‑1 (MQ=255)
                                                                                                                             gccgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGAt                       >  1:76046/8‑150 (MQ=25)
                                                                                                                              tGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATa                      <  1:942603/150‑1 (MQ=255)
                                                                                                                              tGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATa                      <  1:769504/150‑1 (MQ=255)
                                                                                                                               gagaTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCGGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATaa                     >  2:121958/1‑150 (MQ=255)
                                                                                                                               gagaTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACCCCTTATCCCGCCTACGGCTCGGATTTGTAGGCCTGATaa                     >  1:941547/1‑150 (MQ=255)
                                                                                                                                   tGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGAcg                 >  1:244768/1‑150 (MQ=255)
                                                                                                                                    gCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGAcgc                >  1:603788/1‑150 (MQ=255)
                                                                                                                                       ggATGGGGCGTAAACGCCTTATCCGGCCTACGGCTCGGGTTGTAGGCCGGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGGGCCGGATGCAGCGTGAACTCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCGGATAAGACGCGcc             <  1:552421/150‑1 (MQ=37)
                                                                                                                                         aTGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAg           <  2:483464/150‑1 (MQ=255)
                                                                                                                                           gcggcgTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCg         <  2:1052933/150‑1 (MQ=255)
                                                                                                                                             ggcgTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTc       >  1:1414196/1‑150 (MQ=255)
                                                                                                                                               ggTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGc     <  1:389584/149‑1 (MQ=255)
                                                                                                                                                gTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGAGTGAACGCCTTATACGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCgaa    >  1:650927/1‑148 (MQ=255)
                                                                                                                                                 tAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCAt   >  1:350517/1‑150 (MQ=255)
                                                                                                                                                  aaaCGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATc  <  2:170247/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATC  >  NC_000913/2347028‑2347323

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: