Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,923,720 T→C 6.9% intergenic (‑40/+24) rsmG ← / ← mnmG 16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein/5‑methylaminomethyl‑2‑thiouridine modification at tRNA U34

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,923,7200TC 6.9% 218.8 / 6.0 72intergenic (‑40/+24)rsmG/mnmG16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein/5‑methylaminomethyl‑2‑thiouridine modification at tRNA U34
Reads supporting (aligned to +/- strand):  ref base T (27/40);  new base C (2/3);  total (29/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

TTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTT  >  NC_000913/3923574‑3923868
                                                                                                                                                  |                                                                                                                                                    
ttGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGc                                                                                                                                                    <  1:42948/149‑1 (MQ=255)
      cATTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCa                                                                                                                                              <  1:580421/149‑1 (MQ=255)
      cATTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCACCCGCTTACCa                                                                                                                                              <  1:273356/149‑1 (MQ=255)
       aTTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAg                                                                                                                                             >  3:338272/1‑149 (MQ=255)
       aTTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAg                                                                                                                                             >  1:38969/1‑149 (MQ=255)
       aTTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTACTTACCAg                                                                                                                                             <  3:540506/149‑1 (MQ=255)
         ttATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCACCCGCTTACCAGGc                                                                                                                                           <  3:198547/149‑1 (MQ=255)
              cagcaTATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCAtttt                                                                                                                                      >  4:414520/1‑149 (MQ=255)
               agcaTATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCAttttt                                                                                                                                     >  4:395323/1‑149 (MQ=255)
               agcaTATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCAttttt                                                                                                                                     <  1:587847/149‑1 (MQ=255)
                 caTATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTaa                                                                                                                                   >  1:304174/1‑149 (MQ=255)
                     ttCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCg                                                                                                                               >  1:440134/1‑149 (MQ=255)
                                   aaGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTAcc                                                                                                                                               <  3:86098/119‑1 (MQ=255)
                                   aaGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTAcc                                                                                                                                               >  4:86098/1‑119 (MQ=255)
                                   aaGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTacgac                                                                                                                 <  2:376819/149‑1 (MQ=255)
                                     gCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACgc                                                                                                               <  3:576324/149‑1 (MQ=255)
                                       tGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACgcag                                                                                                             >  4:13457/1‑149 (MQ=255)
                                         gTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACgcagca                                                                                                           >  3:118965/1‑149 (MQ=255)
                                            tttCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATAc                                                                                                        <  2:476506/149‑1 (MQ=255)
                                              tCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATAccc                                                                                                      <  3:105004/149‑1 (MQ=255)
                                                tGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTg                                                                                                    >  3:95830/1‑149 (MQ=255)
                                                       cGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTc                                                                                             <  2:142992/149‑1 (MQ=255)
                                                        ggTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCa                                                                                            <  2:215654/149‑1 (MQ=255)
                                                        ggTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCa                                                                                            <  1:341975/149‑1 (MQ=255)
                                                        ggTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCa                                                                                            <  1:212653/149‑1 (MQ=255)
                                                          tAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTAAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGc                                                                                          <  3:586722/149‑1 (MQ=255)
                                                          tAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATAACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCATTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGc                                                                                          >  3:424795/1‑149 (MQ=255)
                                                            aGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCca                                                                                        >  4:88822/1‑149 (MQ=255)
                                                                     ccTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACcagcag                                                                               >  2:227314/1‑149 (MQ=255)
                                                                      cTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGa                                                                              >  4:251868/1‑149 (MQ=255)
                                                                      cTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGa                                                                              >  3:88411/1‑149 (MQ=255)
                                                                      cTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGa                                                                              >  1:547698/1‑149 (MQ=255)
                                                                       tGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGaa                                                                             >  1:489690/1‑149 (MQ=255)
                                                                        gCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAAt                                                                            <  4:88411/149‑1 (MQ=255)
                                                                          gTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACAGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATgg                                                                          <  2:256867/149‑1 (MQ=255)
                                                                             tttCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCa                                                                                            >  2:60795/1‑128 (MQ=255)
                                                                             tttCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCa                                                                                            <  1:60795/128‑1 (MQ=255)
                                                                              ttCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGAt                                                                      <  3:470690/149‑1 (MQ=255)
                                                                              ttCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGAt                                                                      <  1:121016/149‑1 (MQ=255)
                                                                              ttCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGAt                                                                      <  2:229748/149‑1 (MQ=255)
                                                                              ttCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGAt                                                                      <  2:12031/149‑1 (MQ=255)
                                                                              ttCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGAt                                                                      <  2:468891/149‑1 (MQ=255)
                                                                              ttCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGAt                                                                      <  2:563188/149‑1 (MQ=255)
                                                                              ttCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAACGGAGAt                                                                      <  1:283784/149‑1 (MQ=255)
                                                                                cagcagTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATgg                                                                    <  1:75476/149‑1 (MQ=255)
                                                                                    agTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGCGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGc                                                                >  4:237877/1‑149 (MQ=255)
                                                                                          gagaGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGt                                                          >  3:269758/1‑149 (MQ=255)
                                                                                          gagaGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGt                                                          <  3:475833/149‑1 (MQ=255)
                                                                                            gagTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGa                                                        <  2:547698/149‑1 (MQ=255)
                                                                                              gTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAAt                                                                                                                                  <  3:182047/73‑1 (MQ=255)
                                                                                              gTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACg                                                      >  4:47497/1‑149 (MQ=255)
                                                                                              gTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACg                                                      >  1:118154/1‑149 (MQ=255)
                                                                                              gTTTGTTGAGCACGGTGATTACCTGTTCTTGACGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAAt                                                                                                                                  >  4:182047/1‑73 (MQ=255)
                                                                                                     gAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGaaa                                               <  4:118965/149‑1 (MQ=255)
                                                                                                      aGCACGGTGATCACCTGTTGTTGATGCGCTGCCTGGTAAGCGCGTGCTCACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCCTAGCCTGTTTTTTCAGCGACTCTACCAGAATGGAGATGGTGGCAGGCGTGGTGCCAGAAAt                                              <  2:586992/149‑1 (MQ=255)
                                                                                                                   ccTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTgg                                 <  4:95830/149‑1 (MQ=255)
                                                                                                                    cTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGc                                <  1:425112/149‑1 (MQ=255)
                                                                                                                      gTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGgag                                                                        <  2:455052/107‑1 (MQ=255)
                                                                                                                      gTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGgag                                                                        >  1:455052/1‑107 (MQ=255)
                                                                                                                      gTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCg                              >  3:326591/1‑149 (MQ=255)
                                                                                                                        tCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGAt                            <  2:304174/149‑1 (MQ=255)
                                                                                                                         cTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATa                           <  3:379134/149‑1 (MQ=255)
                                                                                                                          ttGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATag                          >  2:25444/1‑149 (MQ=255)
                                                                                                                                 aatGCCTGGTAAGCACCCGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTgg                   >  3:109275/3‑149 (MQ=255)
                                                                                                                                    gCCTGGTAAGCACCCGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGttt                <  4:109275/149‑1 (MQ=255)
                                                                                                                                      cTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTtgt              >  3:200903/1‑149 (MQ=255)
                                                                                                                                        ggTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATa                           <  4:343087/134‑1 (MQ=255)
                                                                                                                                        ggTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATa                           >  3:343087/1‑134 (MQ=255)
                                                                                                                                         gTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGAt           <  3:397489/149‑1 (MQ=255)
                                                                                                                                          tAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATc          <  1:581764/149‑1 (MQ=255)
                                                                                                                                          tAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATc          <  4:119969/149‑1 (MQ=255)
                                                                                                                                          tAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATc          <  1:25444/149‑1 (MQ=255)
                                                                                                                                                  tGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGtt  >  3:87794/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTGGTTTTTCTGGTGATCGGTAAGCGAAATACCTGCGTCTTTCAGCAGTAAGGAGAGTTTGTTGAGCACGGTGATTACCTGTTCTTGATGCGTTGCCTGGTAAGCGGGTGCTTACCAGGCATTTTTAATGCGTTATGCGCTACGACGCAGCATACCCTGTTTTTTCAGCCACACCAGCAGAATGGAGATGGCCGCAGGCGTGACGCCAGAAATACGCGAAGCTTGGCCGATAGAGGCTGGTTTGTGATCGTTAAGTT  >  NC_000913/3923574‑3923868

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: