Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,276,886 C→A 7.9% intergenic (+33/‑116) garD → / → prlF D‑galactarate dehydrogenase/antitoxin of the SohA(PrlF)‑YhaV toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,276,8860CA 7.9% 506.7 / 23.5 149intergenic (+33/‑116)garD/prlFD‑galactarate dehydrogenase/antitoxin of the SohA(PrlF)‑YhaV toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base C (71/66);  new base A (6/6);  total (77/72)
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.63e-01

Only 100 of 149 total aligned reads displayed.

GTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTAC  >  NC_000913/3276745‑3277032
                                                                                                                                              |                                                                                                                                                  
gttgttCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑ATTAGTGc                                                                                                                                             >  2:36571/1‑149 (MQ=255)
 ttgttCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGcc                                                                                                                                            <  2:852854/149‑1 (MQ=255)
   gttCCACTTTATTCTCGACGTCGCCAGCGGGAAGCAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTGTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCAt                                                                                                                                          >  3:928024/1‑149 (MQ=255)
   gttCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCAt                                                                                                                                          <  1:586003/149‑1 (MQ=255)
      ccACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑ACTAAGGCCATGtt                                                                                                                                       <  3:71275/149‑1 (MQ=255)
       cACTTTATTCTCGACGTCGTCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTGTTTTGGAACGGCCCTCGCAAAACATGGC‑ACTAAGCCCATGttt                                                                                                                                      <  4:1033426/149‑1 (MQ=255)
       cACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAAGCCGGCACCGGTGACCTGATTTCTTTTCGAACGTCCCTCGCAAAACATGGCACTAAG‑GCCATGttt                                                                                                                                      <  1:819327/149‑1 (MQ=255)
        aCTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGtttt                                                                                                                                     >  3:133958/1‑149 (MQ=255)
        aCTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGtttt                                                                                                                                     <  2:631830/149‑1 (MQ=255)
          tttATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTa                                                                                                                                   >  4:408904/1‑149 (MQ=255)
          tttATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTa                                                                                                                                   >  4:423304/1‑149 (MQ=255)
          tttATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTa                                                                                                                                   <  4:928024/149‑1 (MQ=255)
          tttATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTa                                                                                                                                   <  1:422887/149‑1 (MQ=255)
          tttATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑ACTAAGGCCATGttttgc                                                                                                                                   <  4:759798/149‑3 (MQ=255)
                ctcGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGttt                                                                                                                             <  3:731411/149‑1 (MQ=255)
                   gACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTaaa                                                                                                                          >  4:342792/1‑149 (MQ=255)
                   gACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTaaa                                                                                                                          <  3:983208/149‑1 (MQ=255)
                   gACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTaaa                                                                                                                          <  1:346714/149‑1 (MQ=255)
                      gTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGcc                                                                                                                       <  2:16289/149‑1 (MQ=255)
                      gTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGcc                                                                                                                       <  3:262059/149‑1 (MQ=255)
                          ccAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCa                                                                                                                   <  2:492703/149‑1 (MQ=255)
                             gCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGt                                                                                                                <  1:884516/149‑1 (MQ=255)
                             gCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGt                                                                                                                <  1:900544/149‑1 (MQ=255)
                                 gaagaagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCAt                                                                                                            >  3:805547/1‑149 (MQ=255)
                                  aagaagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGccc                                                                                                                      >  1:610006/1‑138 (MQ=255)
                                  aagaagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCAtt                                                                                                           <  2:748227/149‑1 (MQ=255)
                                    gaagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaa                                                                                                         <  3:543195/149‑1 (MQ=255)
                                     aagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaat                                                                                                        >  3:199015/1‑149 (MQ=255)
                                     aagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaat                                                                                                        <  1:253715/149‑1 (MQ=255)
                                     aagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaat                                                                                                        <  4:960457/149‑1 (MQ=255)
                                     aagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaat                                                                                                        >  1:469020/1‑149 (MQ=255)
                                     aagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaat                                                                                                        <  2:357215/149‑1 (MQ=255)
                                     aagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCCTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCGATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaat                                                                                                        >  2:70898/1‑149 (MQ=255)
                                       gaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaataa                                                                                                      >  4:1038454/1‑149 (MQ=255)
                                       gaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaataa                                                                                                      >  4:746039/1‑149 (MQ=255)
                                       gaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaataa                                                                                                      >  2:152990/1‑149 (MQ=255)
                                          aaCCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgc                                                                                                   <  4:766832/149‑1 (MQ=255)
                                            ccTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcat                                                                                                 <  1:45666/149‑1 (MQ=255)
                                            ccTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcat                                                                                                 >  1:320951/1‑149 (MQ=255)
                                             cTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatt                                                                                                >  3:642061/1‑149 (MQ=255)
                                              ttCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcaaat                                                                                               <  3:759003/149‑4 (MQ=255)
                                                ctcGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcaaatgc                                                                                             <  4:963713/149‑6 (MQ=255)
                                                     aTCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTAtt                                                                                                                                 <  2:71030/108‑1 (MQ=255)
                                                     aTCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTAc                                                                                        >  2:319778/1‑149 (MQ=255)
                                                     aTCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTAc                                                                                        >  1:11782/1‑149 (MQ=255)
                                                         aTGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTa                                                                                    <  4:624058/149‑1 (MQ=255)
                                                          tGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCAAATGCATTattaatg                                                                                   <  4:668767/149‑7 (MQ=255)
                                                            ggggCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGtt                                                                                 <  3:47474/149‑1 (MQ=255)
                                                               gCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCaa                                                                              >  1:674228/1‑149 (MQ=255)
                                                                cTGCATAACCAGATGGCGGTGTTTAACCCGGCACCGGTGACATGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAg                                                                             <  4:464568/149‑1 (MQ=255)
                                                                 tGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGc                                                                            >  2:515627/1‑149 (MQ=255)
                                                                        ccAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAAt                                                                     >  1:103570/1‑149 (MQ=255)
                                                                          aGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGt                                                                                                                               >  3:374306/1‑89 (MQ=255)
                                                                          aGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGt                                                                                                                               <  4:374306/89‑1 (MQ=255)
                                                                          aGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAAttt                                                                   <  3:423456/149‑1 (MQ=255)
                                                                           gCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGccc                                                                                                                      <  1:278929/97‑1 (MQ=255)
                                                                           gCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGccc                                                                                                                      >  2:278929/1‑97 (MQ=255)
                                                                            cTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGa                                                                 <  4:1007460/149‑1 (MQ=255)
                                                                            cTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGa                                                                 <  1:515627/149‑1 (MQ=255)
                                                                            cTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGa                                                                 <  1:701240/149‑1 (MQ=255)
                                                                              ggCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGc                                                               >  1:800290/1‑149 (MQ=255)
                                                                                 ggTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcattt                                                                                               >  2:76819/1‑114 (MQ=255)
                                                                                 ggTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcattt                                                                                               <  1:76819/114‑1 (MQ=255)
                                                                                  gtgtTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGaa                                                           >  4:218087/1‑149 (MQ=255)
                                                                                  gtgtTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGaa                                                           >  1:516247/1‑149 (MQ=255)
                                                                                       tAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAg                                                      <  2:396225/149‑1 (MQ=255)
                                                                                        aaCCCGGGAGCGGTGACTTGGTTACGTTGGGAACGTCCCTCGCAAAACATGGT‑CTTAGTGCCATGTTTTTATTGTTTAAAGACCCCACGTCCATTTATAATGCATTTGCATTACCTTAGTTCAATCTTATAATTTGAGCAGAAAACAgg                                                     >  2:309256/1‑149 (MQ=255)
                                                                                        aaCCCGGCACCGGTGACCTGAGTTCTTTTGGGACGTCCCTCGCAAGACATAGC‑CTTAGCGCCATGTTTTTATTGTTTAAAGGCCCCACGTCCATTGATAATGCAATTGCCTTACCTTAGTGCAAGCGTATCATTTGAGCAGAAAACAgg                                                     >  2:11523/1‑149 (MQ=255)
                                                                                          cccGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGTCt                                                   >  1:726753/1‑149 (MQ=255)
                                                                                           ccGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTg                                                  >  2:18090/1‑149 (MQ=255)
                                                                                                     tGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACa                                        <  2:800290/149‑1 (MQ=255)
                                                                                                      gACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAg                                       <  2:674228/149‑1 (MQ=255)
                                                                                                         cTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTaaaa                                             >  4:272675/1‑140 (MQ=255)
                                                                                                         cTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGa                                    >  2:649088/1‑149 (MQ=255)
                                                                                                         cTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGa                                    <  1:722920/149‑1 (MQ=255)
                                                                                                          tGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGaa                                   <  3:631075/149‑1 (MQ=255)
                                                                                                           gATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAAt                                  <  2:111173/149‑1 (MQ=255)
                                                                                                            aTTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAGACAGGCTGTAAAAgg                                           <  2:870090/139‑1 (MQ=255)
                                                                                                            aTTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATc                                 <  2:492562/149‑1 (MQ=255)
                                                                                                             tttCTTTTGGAACGTTCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCa                                >  4:774068/1‑149 (MQ=255)
                                                                                                                 ttttGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGcc                            >  1:374496/1‑149 (MQ=255)
                                                                                                                 ttttGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGcc                            >  2:373448/1‑149 (MQ=255)
                                                                                                                 ttttGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGcc                            >  3:310694/1‑149 (MQ=255)
                                                                                                                  tttGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGccc                           <  3:423304/149‑1 (MQ=255)
                                                                                                                    tGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGc                         <  4:915251/149‑1 (MQ=255)
                                                                                                                     ggAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCt                        >  1:674507/1‑149 (MQ=255)
                                                                                                                        aCGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGACAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAAt                     <  4:864652/149‑1 (MQ=255)
                                                                                                                               tCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAAt                                  <  3:78062/129‑1 (MQ=255)
                                                                                                                               tCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAAt                                  >  4:78062/1‑129 (MQ=255)
                                                                                                                               aataaaaaCATGGC‑ACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATgctcgct              >  3:1032029/5‑149 (MQ=255)
                                                                                                                                 taaaaaCATGGC‑ACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGctct            >  1:836792/3‑149 (MQ=255)
                                                                                                                                   aaaaCATGGC‑ACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCa          >  4:842347/1‑149 (MQ=255)
                                                                                                                                    aaaCATGGC‑ACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGcc                            >  4:472664/1‑130 (MQ=255)
                                                                                                                                    aaaCATGGC‑ACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGcc                            <  3:472664/130‑1 (MQ=255)
                                                                                                                                     aaCATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACg        >  3:544100/1‑149 (MQ=255)
                                                                                                                                     aaCATGGC‑ACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACg        >  4:55731/1‑149 (MQ=255)
                                                                                                                                      aCATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGc       <  1:1014099/149‑1 (MQ=255)
                                                                                                                                          ggC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTa   <  1:70898/149‑1 (MQ=255)
                                                                                                                                          ggC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTa   <  2:674507/149‑1 (MQ=255)
                                                                                                                                           gC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTAc  >  2:530143/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTAC  >  NC_000913/3276745‑3277032

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