Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,347,168 C→A 5.1% intergenic (+18/‑216) 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,1680CA 5.1% 314.1 / 5.1 111intergenic (+18/‑216)nrdA/nrdBribonucleoside‑diphosphate reductase 1, alpha subunit/ribonucleoside‑diphosphate reductase 1, beta subunit, ferritin‑like protein
Reads supporting (aligned to +/- strand):  ref base C (52/53);  new base A (4/2);  total (56/55)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.73e-01

Only 100 of 115 total aligned reads displayed.

CACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCG  >  NC_000913/2347021‑2347310
                                                                                                                                                   |                                                                                                                                              
cACCGCCTACAAATTCGGGGTCAAAACACTGTTTTATCAGAACCCCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATgcggc                                                                                                                                             <  2:748831/150‑1 (MQ=255)
 aCCGCCTACAAATTAGGGGTCATAACTCTGTATTATCAGAACACCGGTTACGGCGCTTAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGAGGCGCATGTAAGATCTGATATTGAGATGCCGGATgcggc                                                                                                                                             <  2:1673693/150‑1 (MQ=255)
  ccGCCTATAAATTCGGGGTCAAAACACTGTATAATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGTCGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATgcggcg                                                                                                                                            <  2:544690/150‑1 (MQ=255)
   cGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGt                                                                                                                                           >  1:442101/1‑150 (MQ=255)
    cccgactaaTTCGGGGACAAACCCCTGTATTATCAGTACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGGTGCGAACGCGGCGCATGTAAGATCTGATAATGAGATGCCGGATGCGGCGTa                                                                                                                                          <  1:234790/143‑1 (MQ=255)
         cAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCAAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGc                                                                                                                                     >  2:1685414/1‑150 (MQ=255)
          ataTTCGTGGTCAAAACACTGTATTATCAGAACACCCTTGACGGCGCTGAAGACGCACAAGACGATCTGGTGGCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTGAGATCTTATATTGAGATTCCGGATGCGGCGTAAACGcc                                                                                                                                    <  2:1590652/148‑1 (MQ=255)
          aactacGGGGTCAACACACTGTATTAGCACAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTTGTGCCGTAAATCCAGGACGATGGCTGGGACAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACggc                                                                                                                                    <  2:170657/145‑3 (MQ=255)
           aaTTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACTAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCt                                                                                                                                   <  1:940217/150‑1 (MQ=255)
             ttCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTa                                                                                                                                 <  1:770517/150‑1 (MQ=255)
                 gggTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                                                             <  2:1641867/150‑1 (MQ=255)
                   gTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                           >  2:1184400/1‑150 (MQ=255)
                   gTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                           <  2:1609856/150‑1 (MQ=255)
                   gTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGAGTAAACGCCTTATCCGGc                                                                                                                           >  1:517759/1‑150 (MQ=255)
                    ccAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                          <  1:302103/149‑1 (MQ=255)
                     cAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTAATCCGGCCt                                                                                                                         >  1:363769/1‑150 (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                                                                                                                  >  2:876135/1‑150 (MQ=255)
                            cTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTc                                                                                                                  <  1:1700379/150‑1 (MQ=255)
                             tGTATTATGAGAACACCCGTGCCGGCTCTGAAGACGCACAAGACGATCTGGTTCCGTCAATCCAGGACGATGGCTGCGAAAGCGGAGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCg                                                                                                                 <  2:1287952/150‑1 (MQ=39)
                              gTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCgg                                                                                                                <  1:280402/150‑1 (MQ=255)
                              gTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCgg                                                                                                                >  2:205178/1‑150 (MQ=255)
                                attatCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGtt                                                                                                              <  2:1296967/150‑1 (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)
                                   atCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGt                                                                                                           >  2:1766414/1‑150 (MQ=255)
                                     cAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAg                                                                                                         >  1:1336360/1‑150 (MQ=255)
                                       gAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGc                                                                                                       >  2:1652171/1‑150 (MQ=255)
                                       gAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGc                                                                                                       >  1:1517479/1‑150 (MQ=255)
                                        tacacCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGACGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGcc                                                                                                      <  2:1158287/149‑1 (MQ=255)
                                            ccagTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGAGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAGGGCTCGGTTTGTAGGCCTGAt                                                                                                  >  2:1379704/4‑150 (MQ=255)
                                              cGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATaa                                                                                                >  2:1471911/1‑150 (MQ=255)
                                                  aCGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGAcg                                                                                            >  2:1095554/1‑150 (MQ=255)
                                                    ggCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATAGTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGATTGTAGGCCTGATAATAcgcg                                                                                          >  2:294871/1‑150 (MQ=255)
                                                            aGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTc                                                                                  <  2:248578/150‑1 (MQ=255)
                                                             gACGCACAAGACGAGCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCg                                                                                 <  2:366001/150‑1 (MQ=255)
                                                               cGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCa                                                                               >  2:329096/1‑150 (MQ=255)
                                                               cGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCa                                                                               <  2:1314279/150‑1 (MQ=255)
                                                               cGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCa                                                                               >  1:1648188/1‑150 (MQ=255)
                                                                  acaAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCa                                                                            >  1:1601068/1‑150 (MQ=255)
                                                                   caAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCAGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGATTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAg                                                                           >  2:1369472/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)
                                                                       aCGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATGCGACGCTGGCGCGTCTTATCAGGcct                                                                       <  2:718953/150‑3 (MQ=37)
                                                                        cGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATCGGTTTGTAGGCCTGATAAGACGAGCCAGCGTCGCATCAGGCTcc                                                                      >  2:938882/1‑150 (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)
                                                                                   cGTCAACCCAGGACTATGGCTCCGTACGCGGCGCATGTAAGATCTGATATTTAGATGCCGGATGCGGCTTAATCGGCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGAt                                                           <  2:324600/150‑1 (MQ=25)
                                                                                    gTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATg                                                          <  2:1170469/150‑1 (MQ=255)
                                                                                     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                                   <  2:860444/150‑1 (MQ=255)
                                                                                                           aaaGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGcc                                   <  1:867822/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)
                                                                                                                                         ttGCCGGATGCGGCGTAATCGCCTTTTCCGGCATACGGCTCTGTTTTTATGCCTTATTAGACGCGCCATCGTCGCATCAGGCTCCTGGTGCATGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGAc     <  1:1379704/149‑1 (MQ=25)
                                                                                                                                          tGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTg                                                   >  1:1782732/1‑103 (MQ=255)
                                                                                                                                          tGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTg                                                   <  2:1782732/103‑1 (MQ=255)
                                                                                                                                          tGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGAcg    >  2:481011/1‑150 (MQ=255)
                                                                                                                                            ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGAcgcg  >  2:613948/1‑150 (MQ=255)
                                                                                                                                            ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATGTGTAGGCCTGATAAGAcgcg  >  1:837204/1‑150 (MQ=255)
                                                                                                                                            ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTCAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGAcgcg  >  2:918892/1‑150 (MQ=255)
                                                                                                                                            ccGGATAAGGCGTTTACGCCGCATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGAcgcg  >  2:777290/1‑150 (MQ=37)
                                                                                                                                                   |                                                                                                                                              
CACCGCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCG  >  NC_000913/2347021‑2347310

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: