Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,397,961 T→G 5.7% S140R (AGC→CGC)  yhdE ← dTTP/UTP pyrophosphatase; m(5)UTP/m(5)CTP/pseudo‑UTP pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,397,9610TG 5.7% 360.4 / 9.8 124S140R (AGC→CGC) yhdEdTTP/UTP pyrophosphatase; m(5)UTP/m(5)CTP/pseudo‑UTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base T (53/64);  new base G (3/4);  total (56/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.13e-01

Only 100 of 124 total aligned reads displayed.

CGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGACCCGATAATTTGCGCAACA  >  NC_000913/3397813‑3398104
                                                                                                                                                    |                                                                                                                                               
cGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAGTACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTg                                                                                                                                                >  2:598721/1‑150 (MQ=255)
cGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCCCCTGCTTTATCTAACGGTTCATCGCTg                                                                                                                                                >  2:1094086/1‑150 (MQ=255)
cGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTg                                                                                                                                                >  1:932810/1‑150 (MQ=255)
cGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCCGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTg                                                                                                                                                >  1:214033/1‑150 (MQ=255)
  cAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGc                                                                                                                                              <  1:892233/150‑1 (MQ=255)
   aGTGCGTTTAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTTTCTTCCTGACAAAACAGCCACCCAGCCCCTGTATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCg                                                                                                                                             <  2:606876/150‑1 (MQ=255)
     tGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGAc                                                                                                                                           <  2:351015/150‑1 (MQ=255)
        gTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATa                                                                                                                                        >  1:606124/1‑150 (MQ=255)
            aaaTTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGccc                                                                                                                                    <  1:1355910/150‑1 (MQ=255)
               ttACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCcgcg                                                                                                                                 <  2:772421/150‑1 (MQ=255)
                 aCTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGAGTACGGCGTGATTGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGAt                                                                                                                               <  1:1573318/150‑1 (MQ=255)
                 aCTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGAt                                                                                                                               <  2:1424619/150‑1 (MQ=255)
                 aCTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGAt                                                                                                                               <  2:1556956/150‑1 (MQ=255)
                   ttAATAATTCATACCTTACACCCACCTGAAAGGCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACTCAACAGCCTCCCACCCCCGGAATACCGTCTGCACCTGCTTTATCTAACGGTGCATCGCTGGCTATGTCGCCAGCGATga                                                                                                                             <  1:243653/150‑2 (MQ=255)
                    taataaTTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTc                                                                                                                            >  2:1615300/1‑150 (MQ=255)
                     aataatTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCGGGCTATGTCGCCAGCGATgaac                                                                                                                           <  2:232861/150‑4 (MQ=255)
                     aataatTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCGGGCTATGTCGCCAGCGATgaac                                                                                                                           <  1:512681/150‑4 (MQ=255)
                       taatTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCGGGCTATGTCGCCAGCGATgaaccg                                                                                                                         <  1:46813/150‑6 (MQ=255)
                        aatTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCg                                                                                                                        >  2:939179/1‑150 (MQ=255)
                           tCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCt                                                                                                                     <  1:1327567/150‑1 (MQ=255)
                            cATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTg                                                                                                                    >  2:535843/1‑150 (MQ=255)
                             aTACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACATCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGt                                                                                                                   <  1:1522445/150‑1 (MQ=255)
                             aTACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGt                                                                                                                   <  2:719821/150‑1 (MQ=255)
                               aCGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTa                                                                                                                 <  1:939179/150‑1 (MQ=255)
                                 gTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGCCCTAGCCCGCGATGTCTTCGTCTGATAAc                                                                                                               >  2:1322069/1‑150 (MQ=255)
                                     cAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCGCTGAATACCGTATGCGCCTGCTTTATCTAACGGTTCATCGCTGGCGACGTAGCCCGCGGTGGCTGCGTCTGTGAACGTTc                                                                                                           >  2:843308/1‑150 (MQ=255)
                                     cAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTc                                                                                                           <  1:205864/150‑1 (MQ=255)
                                        ccAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTTTGCCCCTGCTTTATCTAACGGTTCATCGCTGGCGACCTAGCCCGGGATGTATTCGTCTGTTTACGTTCTGa                                                                                                        >  2:1135181/1‑150 (MQ=255)
                                        ccAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGa                                                                                                        >  2:1263292/1‑150 (MQ=255)
                                        ccAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGa                                                                                                        >  1:1710608/1‑150 (MQ=255)
                                          aGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCCCCCAGCCCCTGAATACCGTATGCACCAGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGaaa                                                                                                      <  2:873245/150‑1 (MQ=255)
                                           gCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAg                                                                                                     <  2:357301/150‑1 (MQ=255)
                                           gCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAg                                                                                                     <  2:233869/150‑1 (MQ=255)
                                            cGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCGCCCAGCAGCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTTACGTTCTGAAAGt                                                                                                    >  1:1294070/1‑150 (MQ=255)
                                            cGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGt                                                                                                    <  2:908313/150‑1 (MQ=255)
                                            cGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGt                                                                                                    >  1:752937/1‑150 (MQ=255)
                                                aaGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTcaca                                                                                                >  1:351575/1‑150 (MQ=255)
                                                  gCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATc                                                                                              <  1:1110319/150‑1 (MQ=255)
                                                   ccGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCg                                                                                             >  2:1386496/1‑150 (MQ=255)
                                                      tcTACGGCTTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTg                                                                                          <  1:569475/149‑1 (MQ=255)
                                                      aCTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTg                                                                                          <  1:444685/150‑1 (MQ=255)
                                                      aCTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTg                                                                                          <  2:1585148/150‑1 (MQ=255)
                                                         aCGGCGTGATAGCTGCCATTTATCTTACTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGTAGACATAGCCCTCGATGTCATCGTCTGTTAACGTTCTGAAAGTCACATCGGTGacc                                                                                       >  1:1576340/1‑150 (MQ=255)
                                                          cGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGTCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGacca                                                                                      >  1:1418958/1‑150 (MQ=255)
                                                           ggCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGaccac                                                                                     <  1:1715333/150‑1 (MQ=255)
                                                               tGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATCCCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAgg                                                                                 <  1:1526793/150‑1 (MQ=255)
                                                                     cTGCCATTTATCTTCCTGACAAACCAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGCCCACCAGGCATTCg                                                                           >  2:208129/1‑150 (MQ=255)
                                                                     cTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCg                                                                           <  2:459748/150‑1 (MQ=255)
                                                                        ccATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGGCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCgaaa                                                                        <  1:1140843/150‑3 (MQ=255)
                                                                         cATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCTGTGACCACCAGGCAATCGAGaa                                                                       <  1:283269/150‑1 (MQ=255)
                                                                          aTTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAAt                                                                      <  2:432840/150‑1 (MQ=255)
                                                                           tttATCTTCCTGACAAAACAGCAACCCAGCCACTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAAtg                                                                     <  2:1602610/150‑1 (MQ=255)
                                                                             tATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCCCCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAAtgtg                                                                   >  2:1497231/1‑150 (MQ=255)
                                                                             tATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAAtgtg                                                                   >  2:513371/1‑150 (MQ=255)
                                                                             tATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAAtgtg                                                                   <  1:278089/150‑1 (MQ=255)
                                                                                     tGACAAAACAGCCACCCAGCCCCTGAATACAGTATGAACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGACTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCAACACGCAAACGAGAATGTGCTGGCTGt                                                           >  1:541666/1‑150 (MQ=255)
                                                                                      gACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTc                                                          >  2:28873/1‑150 (MQ=255)
                                                                                      gACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTc                                                          >  2:1255199/1‑150 (MQ=255)
                                                                                       aCAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCg                                                         >  1:934451/1‑150 (MQ=255)
                                                                                        cAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGTCATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCTTTAGTGAATTTGCTGGCTGTCgg                                                        >  2:1008122/1‑150 (MQ=255)
                                                                                        cAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCgg                                                        >  2:1669814/1‑150 (MQ=255)
                                                                                               gCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGc                                                 <  1:856022/150‑1 (MQ=255)
                                                                                                ccACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGcc                                                >  1:351167/1‑150 (MQ=255)
                                                                                                ccACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGcc                                                <  1:1496055/150‑1 (MQ=255)
                                                                                                  aCCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCAc                                              <  1:535843/150‑1 (MQ=255)
                                                                                                   cccAGCCCCTGAGTACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCAct                                             <  1:598721/150‑1 (MQ=255)
                                                                                                       gCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCActgctg                                         <  2:1108587/150‑1 (MQ=255)
                                                                                                         cccTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGtc                                       >  1:350497/1‑150 (MQ=255)
                                                                                                                 aCCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTg                               >  1:192481/1‑150 (MQ=255)
                                                                                                                 aCCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTg                               <  2:577361/150‑1 (MQ=255)
                                                                                                                 aCCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGACTGTCGGCCAACGCCACTGCTGTCATCACCTg                               <  1:1637436/150‑1 (MQ=255)
                                                                                                                  ccGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGa                              >  2:46297/1‑150 (MQ=255)
                                                                                                                  ccGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGa                              >  1:1159460/1‑150 (MQ=255)
                                                                                                                    gTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATg                            <  2:924873/150‑1 (MQ=255)
                                                                                                                    gTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATg                            <  1:938588/150‑1 (MQ=255)
                                                                                                                    gTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATg                            <  1:415899/150‑1 (MQ=255)
                                                                                                                    gTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATg                            <  2:932810/150‑1 (MQ=255)
                                                                                                                      aTGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATggg                          <  2:466354/150‑1 (MQ=255)
                                                                                                                       tGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGt                         <  1:1322069/150‑1 (MQ=255)
                                                                                                                       tGCACCTGCGTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGt                         <  1:1266389/150‑1 (MQ=255)
                                                                                                                        gCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTc                        <  2:707116/150‑1 (MQ=255)
                                                                                                                        gCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTc                        <  2:752937/150‑1 (MQ=255)
                                                                                                                        gCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTc                        <  1:332429/150‑1 (MQ=255)
                                                                                                                         cACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCGTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGCCCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCt                       >  1:70683/1‑150 (MQ=255)
                                                                                                                           ccTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACAACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGa                     >  1:1136760/1‑150 (MQ=255)
                                                                                                                            cTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGat                    <  1:1386496/150‑2 (MQ=255)
                                                                                                                             tGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAATGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCAAACGCCACTGCTGTCATCACCTGATGGGTCTGAcc                   <  1:1263292/150‑1 (MQ=255)
                                                                                                                              gCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGAccc                  >  1:1483567/1‑150 (MQ=255)
                                                                                                                              gCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGAccc                  >  1:785939/1‑150 (MQ=255)
                                                                                                                                 ttATCTAACGGTTCGTCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGACCCGaa               >  1:641862/1‑149 (MQ=255)
                                                                                                                                 ttATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGACCCGAt               >  1:1181055/1‑150 (MQ=255)
                                                                                                                                     cTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGACCCGATAAtt           >  2:167182/1‑150 (MQ=255)
                                                                                                                                       aaCGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGACCCGATAATTTg         >  2:1490761/1‑150 (MQ=255)
                                                                                                                                       aaCGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGACCCGATAATTTg         <  2:2031/150‑1 (MQ=255)
                                                                                                                                        cgaagacATCGCGGGCTATGTCGCCAGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGACCCGATAATTTgc        >  1:664413/7‑150 (MQ=255)
                                                                                                                                          ggTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGGCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGACCCGATAATTTgcgc      <  1:1438527/150‑1 (MQ=255)
                                                                                                                                            ttCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGACCCGATAATTTGCGaaa    >  2:1489226/1‑150 (MQ=255)
                                                                                                                                              cATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCCCCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCCCCTGATGGGTCTGACCCGATAATTTGCGCAACa  >  2:1084840/1‑150 (MQ=255)
                                                                                                                                              cATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCATTCGAGAATGTGCTGGCTGTCGGCCAACGACACTGCTGTCATCACCTGATGGGTCTGCCCCGATAATTTGCGCAACa  >  1:1040026/1‑150 (MQ=255)
                                                                                                                                              cATCGCGGGCTATGTCGCCAGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGACCCGATAATTTGCGCAACa  >  2:754121/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CGCAGTGCGTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGACCCGATAATTTGCGCAACA  >  NC_000913/3397813‑3398104

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: