Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,347,181 T→G 5.1% intergenic (+31/‑203) nrdA → / → nrdB ribonucleoside‑diphosphate reductase 1, alpha subunit/ribonucleoside‑diphosphate reductase 1, beta subunit, ferritin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,347,1810TG 5.1% 355.3 / 8.2 118intergenic (+31/‑203)nrdA/nrdBribonucleoside‑diphosphate reductase 1, alpha subunit/ribonucleoside‑diphosphate reductase 1, beta subunit, ferritin‑like protein
Reads supporting (aligned to +/- strand):  ref base T (57/55);  new base G (4/2);  total (61/57)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-01

Only 100 of 120 total aligned reads displayed.

AATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAG  >  NC_000913/2347032‑2347331
                                                                                                                                                     |                                                                                                                                                      
aaTTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACTAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCt                                                                                                                                                        <  1:940217/150‑1 (MQ=255)
  ttCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTa                                                                                                                                                      <  1:770517/150‑1 (MQ=255)
        gTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                                <  2:1609856/150‑1 (MQ=255)
        gTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                                >  2:1184400/1‑150 (MQ=255)
        gTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGAGTAAACGCCTTATCCGGc                                                                                                                                                >  1:517759/1‑150 (MQ=255)
         ccAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                               <  1:302103/149‑1 (MQ=255)
          cAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATAAGGCGTTTACGCCGCATCCGGcat                                                                                                                                              <  1:268738/150‑3 (MQ=255)
           aaacCACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                                                                                             <  1:450245/150‑1 (MQ=255)
           aaaaCACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                                                                                             >  1:797366/1‑150 (MQ=255)
                 cTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTc                                                                                                                                       >  1:1278212/1‑150 (MQ=255)
                 cTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTc                                                                                                                                       <  1:1700379/150‑1 (MQ=255)
                 cTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTc                                                                                                                                       >  2:876135/1‑150 (MQ=255)
                  tGTATTATGAGAACACCCGTGCCGGCTCTGAAGACGCACAAGACGATCTGGTTCCGTCAATCCAGGACGATGGCTGCGAAAGCGGAGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCg                                                                                                                                      <  2:1287952/150‑1 (MQ=39)
                   gTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCgg                                                                                                                                     <  1:280402/150‑1 (MQ=255)
                   gTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCgg                                                                                                                                     >  2:205178/1‑150 (MQ=255)
                      ttatCGGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGttt                                                                                                                                  <  2:874954/150‑1 (MQ=255)
                       tatCAGCACACCCGTGCCGGCGCTGAAGACGCACAAGACGAGCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTATGGCTCGGTTTg                                                                                                                                 <  1:1168817/150‑1 (MQ=255)
                       tatCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTCAGATGCCGGAAGCGGCGTAAACGCCTTATCCGGCCTACAGCTCGGTTTg                                                                                                                                 >  2:321456/1‑150 (MQ=255)
                        atCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGt                                                                                                                                >  2:1759526/1‑150 (MQ=255)
                         tCAGAACACCCGTGACGGCGCTGAAGCCGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTa                                                                                                                               <  2:1746760/150‑1 (MQ=255)
                            gAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGc                                                                                                                            >  1:1517479/1‑150 (MQ=255)
                            gAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGc                                                                                                                            >  2:1652171/1‑150 (MQ=255)
                             tacacCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGACGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGcc                                                                                                                           <  2:1158287/149‑1 (MQ=255)
                                 ccagTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGAGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAGGGCTCGGTTTGTAGGCCTGAt                                                                                                                       >  2:1379704/4‑150 (MQ=255)
                                       aCGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGAcg                                                                                                                 >  2:1095554/1‑150 (MQ=255)
                                         ggCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGACTACGGCTCGGTTTGTAGGCCTCATAAGAcgcg                                                                                                               >  2:1394075/1‑150 (MQ=255)
                                         ggCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATAGTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGATTGTAGGCCTGATAATAcgcg                                                                                                               >  2:294871/1‑150 (MQ=255)
                                                 aGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTc                                                                                                       >  2:656682/1‑150 (MQ=255)
                                                 aGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTc                                                                                                       <  2:248578/150‑1 (MQ=255)
                                                  gACGCACAAGACGAGCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCg                                                                                                      <  2:366001/150‑1 (MQ=255)
                                                    cGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCa                                                                                                    <  2:1314279/150‑1 (MQ=255)
                                                    cGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCa                                                                                                    >  1:1648188/1‑150 (MQ=255)
                                                          aGACGCTCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATCTTGCGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGc                                                                                              <  2:1535619/150‑1 (MQ=255)
                                                          aGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGACTGATAAGACGCGCCAGCGACGCATCAGGc                                                                                              >  2:671560/1‑150 (MQ=255)
                                                           gACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATGCGACGCTGGCGCGTCTTATCAGGcc                                                                                             <  1:226972/150‑2 (MQ=37)
                                                           gACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCt                                                                                             >  2:1463690/1‑150 (MQ=255)
                                                             cGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATCGGTTTGTAGGCCTGATAAGACGAGCCAGCGTCGCATCAGGCTcc                                                                                           >  2:938882/1‑150 (MQ=255)
                                                                tCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCggg                                                                                        <  2:1602379/150‑1 (MQ=255)
                                                                tCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCggg                                                                                        <  1:1735356/150‑1 (MQ=255)
                                                                 cTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGt                                                                                       >  1:1540314/1‑150 (MQ=255)
                                                                   ggTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCCGGCTCCGGGTGc                                                                                     >  1:340981/1‑150 (MQ=255)
                                                                   ggTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAAGGCTCGGTTTGTAGGCCGGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGc                                                                                     >  2:266047/1‑150 (MQ=255)
                                                                    gTGCCGGCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGcc                                                                                    <  1:42149/150‑1 (MQ=255)
                                                                     tGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCg                                                                                   <  2:1307519/150‑1 (MQ=255)
                                                                     ggCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCg                                                                                   <  2:769553/149‑1 (MQ=255)
                                                                      gCCGTCAATCCAGGACTATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCgg                                                                                  <  2:813009/150‑1 (MQ=255)
                                                                      gCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCgg                                                                                  >  2:228430/1‑150 (MQ=255)
                                                                        cGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGAt                                                                                >  1:422873/1‑150 (MQ=255)
                                                                        cGTCAACCCAGGACTATGGCTCCGTACGCGGCGCATGTAAGATCTGATATTTAGATGCCGGATGCGGCTTAATCGGCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGAt                                                                                <  2:324600/150‑1 (MQ=25)
                                                                          tCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGc                                                                              <  1:1095554/150‑1 (MQ=255)
                                                                          tCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGc                                                                              <  1:649506/150‑1 (MQ=255)
                                                                           cAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCa                                                                             <  1:804729/150‑1 (MQ=255)
                                                                             aTCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGc                                                                           <  2:1253158/150‑1 (MQ=255)
                                                                             aTCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGGTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGc                                                                           <  1:702747/150‑1 (MQ=255)
                                                                              tCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCg                                                                          >  1:550118/1‑150 (MQ=255)
                                                                               cccggACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGt                                                                         <  2:915975/147‑1 (MQ=255)
                                                                               ccAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTCCCGGATGCaccct                                                                         >  1:934140/1‑146 (MQ=255)
                                                                                cAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTg                                                                        >  2:494801/1‑150 (MQ=255)
                                                                                cAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTg                                                                        >  2:1571328/1‑150 (MQ=255)
                                                                                   gACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAAc                                                                     >  1:1538249/1‑150 (MQ=255)
                                                                                    aCGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACg                                                                    <  1:554236/150‑1 (MQ=255)
                                                                                      gATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGcc                                                                  <  2:394188/150‑1 (MQ=255)
                                                                                       aTGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCt                                                                 >  1:424296/1‑150 (MQ=255)
                                                                                         ggCTGCGAAAGCGGCGCATGTAAGATCTGATATGGCGAGGCCGGCTGCGGCGGAAACGCCTTATCCGGCCTACGGCTCGGGTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTa                                                               <  1:312699/150‑1 (MQ=37)
                                                                                            ggCGAAAGCGGCGCCTGTAAGAGCTGATATTGAGATGCCGGCTGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGCCGGGTGCCGGATGCAGCGTGAACGCCTTATcc                                                            <  2:1126339/149‑1 (MQ=255)
                                                                                              cGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCgg                                                          >  2:1696972/1‑150 (MQ=255)
                                                                                                acaGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATTCGTCGTAAACGCCTAATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGcc                                                        <  1:181416/148‑1 (MQ=255)
                                                                                                aaaGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGcc                                                        <  1:1479975/150‑1 (MQ=255)
                                                                                                   gcggcgCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTAc                                                     <  2:121204/150‑1 (MQ=255)
                                                                                                     ggcgCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACgg                                                   >  2:1748591/1‑150 (MQ=255)
                                                                                                       cgcATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCt                                                 >  1:201000/1‑150 (MQ=255)
                                                                                                         cATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTAATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCg                                               <  2:363769/150‑1 (MQ=255)
                                                                                                          aTGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCgg                                              >  2:34970/1‑150 (MQ=255)
                                                                                                            gTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGAt                                            >  2:141021/1‑150 (MQ=255)
                                                                                                            gTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGACGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCggtt                                            >  1:1156105/1‑148 (MQ=255)
                                                                                                                 aTCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTa                                       >  1:1786674/1‑150 (MQ=255)
                                                                                                                     gATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGcc                                   <  1:1058172/150‑1 (MQ=255)
                                                                                                                       tatTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGTAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTg                                 <  2:797366/150‑1 (MQ=255)
                                                                                                                          tGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATa                              >  2:847711/1‑150 (MQ=255)
                                                                                                                              ttGCCGGATGCGGCGTAATCGCCTTTTCCGGCATACGGCTCTGTTTTTATGCCTTATTAGACGCGCCATCGTCGCATCAGGCTCCTGGTGCATGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGAc                          <  1:1379704/149‑1 (MQ=25)
                                                                                                                               tGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTg                                                                        <  2:1782732/103‑1 (MQ=255)
                                                                                                                                 ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGAcgcg                       >  2:613948/1‑150 (MQ=255)
                                                                                                                                 ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATGTGTAGGCCTGATAAGAcgcg                       >  1:837204/1‑150 (MQ=255)
                                                                                                                                 ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTCAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGAcgcg                       >  2:918892/1‑150 (MQ=255)
                                                                                                                                 ccGGATAAGGCGTTTACGCCGCATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGGTCGGATTTGTAGGCCTGATAAGAcgcg                       >  2:360850/1‑150 (MQ=37)
                                                                                                                                 ccGGATAAGGCGTTTACGCCGCATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGAcgcg                       >  2:777290/1‑150 (MQ=37)
                                                                                                                                 ccGGATAAGGCGTTTACGCCGCATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGAGTTGTAGGCCTGATAAGAcgcg                       >  1:1768073/1‑150 (MQ=37)
                                                                                                                                         ggcgTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTc               >  1:1711220/1‑150 (MQ=255)
                                                                                                                                         ggcgTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGCATGCAGCGAGAACTCCTTATCAGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACCCCCCAGCgaa               >  1:247676/1‑148 (MQ=25)
                                                                                                                                          gcgTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCCTCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCg              <  2:133639/150‑1 (MQ=255)
                                                                                                                                           cgTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGc             >  1:398959/1‑150 (MQ=255)
                                                                                                                                            gTAAACTCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCa            <  2:302684/150‑1 (MQ=255)
                                                                                                                                            gTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCa            <  2:1278212/150‑1 (MQ=255)
                                                                                                                                            gTAAACGCCTTAGCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCa            <  2:485174/150‑1 (MQ=255)
                                                                                                                                             tAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCAt           >  2:1767253/1‑150 (MQ=255)
                                                                                                                                               aaCGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCa         >  2:1454344/1‑150 (MQ=255)
                                                                                                                                                    cTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGcac    >  2:161817/1‑150 (MQ=255)
                                                                                                                                                    cTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGcac    >  2:1341967/1‑150 (MQ=255)
                                                                                                                                                    cTTATCCGGCCTACGGCTCGGTGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGcac    <  1:575469/150‑1 (MQ=255)
                                                                                                                                                     ttATCCGGCCTACGGAGCGGTTTGGAGGCCTGATAAGACGCGCAGCGTCGCATCAGGCTCAGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGACTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAg  <  2:131785/150‑1 (MQ=25)
                                                                                                                                                     |                                                                                                                                                      
AATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAG  >  NC_000913/2347032‑2347331

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