Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,913,684 A→T 5.8% intergenic (‑33/+15) relA ← / ← rlmD GDP/GTP pyrophosphokinase/23S rRNA m(5)U1939 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,913,6840AT 5.8% 324.5 / 6.2 101intergenic (‑33/+15)relA/rlmDGDP/GTP pyrophosphokinase/23S rRNA m(5)U1939 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (46/49);  new base T (3/3);  total (49/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 105 total aligned reads displayed.

CCAGGTTTCGGCTAAGCACTCACACGACTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCAGGGTT  >  NC_000913/2913535‑2913833
                                                                                                                                                     |                                                                                                                                                     
ccAGGTTTCGGCTAAGCACTCACACGACTTCTGGCTGGTAATACCCAGACTTGCGATACATTTTTCCGGATCAAATTCACCAGACTTATTGATATGTGCACTTCTTACAGAAACCATCGTCATCTCCTTTAGGGACCAGACCTGCCGaaa                                                                                                                                                       >  2:136801/1‑150 (MQ=255)
ccAGGTTTCGGCTAAGCACTCACACGACTTCTGGCTGGTAATACCCAGACTTGCCATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGaaa                                                                                                                                                       >  2:1457504/1‑150 (MQ=255)
 cAGGTTTCGGCTAAGCACTCACACGACTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAAt                                                                                                                                                      >  1:742642/1‑150 (MQ=255)
  tgggttCGGCTAAGCACTCACACGACTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGGACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATc                                                                                                                                                     <  1:371547/146‑1 (MQ=255)
   ggTTTCGGCTAAGCACTCACACGACTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCg                                                                                                                                                    >  1:254341/1‑150 (MQ=255)
     tttCGGCTAAGCACTCACACGACTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGc                                                                                                                                                  <  1:32213/150‑1 (MQ=255)
      ttCGGCTAAGCACTCACACGACTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCa                                                                                                                                                 <  2:1889910/150‑1 (MQ=255)
         ggCTAAGCACTCACACGACTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACAAGACCTGCCGAAATCGGCAAAt                                                                                                                                              >  2:745184/1‑150 (MQ=255)
           cTAAGCACTCACACGACTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGATTTCGGCaggtct                                                                                                                                            <  1:133363/150‑6 (MQ=255)
           cTAAGCACTCACACGACTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCg                                                                                                                                            <  1:822284/150‑1 (MQ=255)
            tAAGCACTCACACGACTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGATTTCGGCaggtctg                                                                                                                                           <  1:1996532/150‑7 (MQ=255)
             aaGCACTCACACGACTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGATTTCGGCaggtctgg                                                                                                                                          <  1:1653557/150‑8 (MQ=255)
               gCACTCACACGACTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAAc                                                                                                                                        >  1:1714590/1‑150 (MQ=255)
                 aCTCACACGACTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTa                                                                                                                                      >  1:1162124/1‑150 (MQ=255)
                     acacGACTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTa                                                                                                                                  >  2:628381/1‑150 (MQ=255)
                       acGACTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAGATCGCAACTATTTAAc                                                                                                                                <  2:227445/150‑1 (MQ=255)
                         gACTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAAcgc                                                                                                                              <  1:1125649/150‑1 (MQ=255)
                          aCTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAAcgcg                                                                                                                             >  1:1724178/1‑150 (MQ=255)
                          aCTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAAcgcg                                                                                                                             <  2:1910745/150‑1 (MQ=255)
                           cTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAAcgcgc                                                                                                                            <  1:1672979/150‑1 (MQ=255)
                           cTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAAcgcgc                                                                                                                            <  2:145111/150‑1 (MQ=255)
                           cTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCAGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAAcgcgc                                                                                                                            >  1:418388/1‑150 (MQ=255)
                            ttCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAAcgcgcg                                                                                                                           <  2:1183610/150‑1 (MQ=255)
                              ctggctggTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCgac                                                                                                                         >  2:1862006/1‑149 (MQ=255)
                                ggctggTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACGATTTAACGCGCGAGaa                                                                                                                       >  2:985263/1‑150 (MQ=255)
                                     gTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTa                                                                                                                  <  1:400022/150‑1 (MQ=255)
                                       aaTACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTATCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTAcc                                                                                                                >  2:383833/1‑150 (MQ=255)
                                       aaTACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTAcc                                                                                                                <  1:1416880/150‑1 (MQ=255)
                                          aCCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATc                                                                                                             >  1:1161078/1‑150 (MQ=255)
                                           cccAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCg                                                                                                            >  1:321913/1‑150 (MQ=255)
                                                  ttGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAg                                                                                                     <  1:628381/150‑1 (MQ=255)
                                                  ttGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAg                                                                                                     >  2:234738/1‑150 (MQ=255)
                                                     ccaTCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGacg                                                                                                  <  1:1457504/148‑3 (MQ=255)
                                                       aTCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTc                                                                                                <  2:1122913/150‑1 (MQ=255)
                                                          cATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCg                                                                                             >  2:1764796/1‑150 (MQ=255)
                                                            tttttCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCgtg                                                                                           <  2:342890/150‑1 (MQ=255)
                                                            tttttCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCgtg                                                                                           <  1:1053538/150‑1 (MQ=255)
                                                                 ccGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTgg                                                                                      >  2:1350802/1‑150 (MQ=255)
                                                                  cGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTggg                                                                                     <  1:881209/150‑1 (MQ=255)
                                                                     aTCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAAc                                                                                  >  2:17775/1‑150 (MQ=255)
                                                                     aTCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAAc                                                                                  <  2:831065/150‑1 (MQ=255)
                                                                     aTCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAAc                                                                                  >  1:1105877/1‑150 (MQ=255)
                                                                      tCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTGCAGATGTCCCGTGTGTGGGAACa                                                                                 <  2:418388/150‑1 (MQ=255)
                                                                      tCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACa                                                                                 <  2:742642/150‑1 (MQ=255)
                                                                      tCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACa                                                                                 <  2:137710/150‑1 (MQ=255)
                                                                       cAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACat                                                                                <  2:1883479/150‑1 (MQ=255)
                                                                         aaTTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACatat                                                                              >  2:1984911/1‑150 (MQ=255)
                                                                         aaTTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACatat                                                                              >  1:988974/1‑150 (MQ=255)
                                                                             cACCAGTCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAg                                                                          >  2:642910/1‑150 (MQ=255)
                                                                               ccAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCa                                                                        <  2:1162124/150‑1 (MQ=255)
                                                                                  gCCTTATTGATATGTGCACTTCTTACCGTAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCg                                                                     <  1:136801/150‑1 (MQ=255)
                                                                                    cTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGcc                                                                   <  1:1613882/150‑1 (MQ=255)
                                                                                    cTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGcc                                                                   <  2:254341/150‑1 (MQ=255)
                                                                                     ttATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCa                                                                  >  2:209210/1‑150 (MQ=255)
                                                                                      tATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAg                                                                 >  1:360135/1‑150 (MQ=255)
                                                                                        ttGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTc                                                               <  2:1093036/150‑1 (MQ=255)
                                                                                        ttGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCAACGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTc                                                               >  1:248709/1‑150 (MQ=255)
                                                                                         tGATATGTGCACTTCTTACCGCAACTATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTcg                                                              >  1:1998229/1‑150 (MQ=255)
                                                                                         tGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTcg                                                              <  1:987628/150‑1 (MQ=255)
                                                                                         tGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTcg                                                              <  1:1752943/150‑1 (MQ=255)
                                                                                          gATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTcgc                                                             >  2:1675923/1‑150 (MQ=255)
                                                                                            tatGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCgcgc                                                           <  1:1388989/150‑1 (MQ=255)
                                                                                            tatGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCATGTGTGGGAACATATCCAGCATCGCCAGTCgcgc                                                           >  1:1114450/1‑150 (MQ=255)
                                                                                              tgtgCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCaa                                                         >  1:1509965/1‑150 (MQ=255)
                                                                                                tgCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATg                                                       <  1:216954/150‑1 (MQ=255)
                                                                                                 gCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATgg                                                      <  1:673259/150‑1 (MQ=255)
                                                                                                   aCTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGta                                                    <  1:1764796/150‑1 (MQ=255)
                                                                                                    cttcttACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGtat                                                   >  1:1147727/1‑150 (MQ=255)
                                                                                                     ttcttACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGtata                                                  <  1:1500846/150‑1 (MQ=255)
                                                                                                        ttACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATcc                                               >  2:196974/1‑150 (MQ=255)
                                                                                                          aCCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTg                                             >  1:360889/1‑150 (MQ=255)
                                                                                                           ccGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGc                                            >  2:1839374/1‑150 (MQ=255)
                                                                                                             gCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCtt                                          <  2:63160/150‑1 (MQ=255)
                                                                                                             gCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCtt                                          >  2:244811/1‑150 (MQ=255)
                                                                                                                 ccATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTtaa                                      <  1:1350802/150‑1 (MQ=255)
                                                                                                                  cATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTtaat                                     <  2:1471256/150‑1 (MQ=255)
                                                                                                                  cATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTtaat                                     <  1:182256/150‑1 (MQ=255)
                                                                                                                          tctcCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTc                             <  2:691985/150‑1 (MQ=255)
                                                                                                                          tctcCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTc                             >  1:403294/1‑150 (MQ=255)
                                                                                                                           ctcCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCg                            >  2:1899193/1‑150 (MQ=255)
                                                                                                                             ccTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCt                          <  1:1862006/150‑1 (MQ=255)
                                                                                                                             ccTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCt                          <  2:1714590/150‑1 (MQ=255)
                                                                                                                              cTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCtt                                          >  1:1728630/1‑133 (MQ=255)
                                                                                                                              cTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCtt                                          <  2:1728630/133‑1 (MQ=255)
                                                                                                                                   gggACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCg                    >  1:1433843/1‑150 (MQ=255)
                                                                                                                                   gggACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCg                    >  2:822361/1‑150 (MQ=255)
                                                                                                                                     gACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAAc                                  >  2:748615/1‑134 (MQ=255)
                                                                                                                                     gACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAAc                                  <  1:748615/134‑1 (MQ=255)
                                                                                                                                     gACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAg                  <  1:1990093/150‑1 (MQ=255)
                                                                                                                                     gACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAg                  >  2:1426027/1‑150 (MQ=255)
                                                                                                                                        gcgatttGCCGATTTCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCGCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCa               >  1:815209/7‑150 (MQ=255)
                                                                                                                                        gcgatttGCCGATTTCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGccc               >  1:1393805/7‑149 (MQ=255)
                                                                                                                                         cgatttGCCGATTTCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAg              >  2:1036333/6‑150 (MQ=255)
                                                                                                                                            ccTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGt           >  1:707390/1‑150 (MQ=255)
                                                                                                                                              tGCCGAAATCGGCAAATCGCAACTATTTAACGCGCTAGAAAAGTACCATCGATTCCAGATGTCCCGTGTCTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTg         <  2:1621982/150‑1 (MQ=255)
                                                                                                                                              tGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTg         >  1:899400/1‑150 (MQ=255)
                                                                                                                                               gCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGc        <  2:1509965/150‑1 (MQ=255)
                                                                                                                                                ccGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCa       <  2:411380/150‑1 (MQ=255)
                                                                                                                                                   aaaTCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCAggg    >  2:536868/1‑150 (MQ=255)
                                                                                                                                                     aTCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCAACGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCAGGGtt  <  2:248709/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCAGGTTTCGGCTAAGCACTCACACGACTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCAGGGTT  >  NC_000913/2913535‑2913833

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