Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,417,471 T→G 6.9% intergenic (‑85/+17) recE ← / ← racC Rac prophage; exonuclease VIII, 5' to 3' specific dsDNA exonuclease/Rac prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,417,4710TG 6.9% 252.1 / 6.4 84intergenic (‑85/+17)recE/racCRac prophage; exonuclease VIII, 5' to 3' specific dsDNA exonuclease/Rac prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (41/37);  new base G (2/4);  total (43/41)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTCCGATATGCACCAGGCACATTTTTTTGACACCTTCATCAAGTGTTTTTACGTTG  >  NC_000913/1417328‑1417617
                                                                                                                                               |                                                                                                                                                  
gcaggcgcACCGGATGAGTTTTTCGCGTTCCGTACCAGGAAGAGTGGTCTTGGGCGCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGAGCTTGCCGGATACACACTGGTCATAGCCTGCGCCATACGCa                                                                                                                                               <  4:70534/142‑1 (MQ=255)
  aGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTGTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAACTTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAg                                                                                                                                              <  4:291191/148‑1 (MQ=255)
   ggTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGc                                                                                                                                            >  7:168700/1‑149 (MQ=255)
    gTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGc                                                                                                                                            >  8:1163/1‑148 (MQ=255)
    gTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTGTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTAGATTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGc                                                                                                                                            >  4:102268/1‑148 (MQ=255)
     acacttgatgaaggtgtcaaaaaaatgtgcctggtgcatatcggaaagaatctgcccgcagaaaaaaaagccgaatttctggcaactctaattgcaatgaaattaaaaggtgaaatctgaaagaaATAGCCTGCGCCATACGCAGGCTa                                                                                                                                          >  4:156302‑M1/126‑149 (MQ=255)
        aCCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTAtt                                                                                                                                        <  1:6109/148‑1 (MQ=255)
             atgatTTTTTCGATTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTAtttctttc                                                                                                                                  <  2:310177/149‑1 (MQ=255)
              tgatTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCa                                                                                                                                 <  8:168700/149‑1 (MQ=255)
                atTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGa                                                                                                                               <  7:1163/149‑1 (MQ=255)
                 ttttttCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGAt                                                                                                                              >  3:72307/1‑149 (MQ=255)
                 ttttttCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGAt                                                                                                                              >  2:166415/1‑149 (MQ=255)
                     ttCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCa                                                                                                                          >  2:41934/1‑149 (MQ=255)
                        gCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCt                                                                                                                       >  3:38620/1‑149 (MQ=255)
                              cGTAACAGGAAGAGGGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAAt                                                                                                                 <  5:282498/149‑1 (MQ=255)
                               gTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAAtt                                                                                                                >  8:105276/1‑149 (MQ=255)
                               gTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGAGCCATACGCAGGCTATTTCTTTAAGATTTCACCTTTTAAtt                                                                                                                >  3:65349/1‑149 (MQ=255)
                                tAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAAttt                                                                                                               <  7:303064/149‑1 (MQ=255)
                                 aaCAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAAttt                                                                                                               >  5:128998/1‑148 (MQ=255)
                                    atgAAGAGTGGTTTTGTGCGCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCAtt                                                                                                           <  2:37029/147‑1 (MQ=255)
                                    aGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCAtt                                                                                                           >  8:95191/1‑149 (MQ=255)
                                    aGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCAtt                                                                                                           <  1:166415/149‑1 (MQ=255)
                                    aGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCAtt                                                                                                           >  2:256302/1‑149 (MQ=255)
                                    aGGAAGAGTGGTTTTGTGCTCATTGTTATTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTCCGCCATACGCAGGCTATATCGTTCACATTTCACCTTTTAATTTCAtt                                                                                                           >  3:101297/1‑149 (MQ=255)
                                    aGGAAGAGTGGGTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCAtt                                                                                                           <  6:58186/149‑1 (MQ=255)
                                     ggAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTg                                                                                                          >  4:14824/1‑149 (MQ=255)
                                      gAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGTCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTg                                                                                                          >  6:67178/1‑148 (MQ=255)
                                      gAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTg                                                                                                          >  3:185312/1‑148 (MQ=255)
                                           gTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTa                                                                                                    <  7:285079/149‑1 (MQ=255)
                                           gTGGTTTAGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGCTACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTa                                                                                                    <  4:41000/149‑1 (MQ=255)
                                              gTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTaga                                                                                                  <  3:14824/148‑1 (MQ=255)
                                                    tgCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCa                                                                                           <  4:185312/149‑1 (MQ=255)
                                                      cTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGTATGGCGCAGGCTAt                                                                                                                                         <  8:11029/101‑1 (MQ=255)
                                                      cTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGTATGGCGCAGGCTAt                                                                                                                                         >  7:11029/1‑101 (MQ=255)
                                                       tCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCa                                                                                                                                               <  2:186254/94‑1 (MQ=255)
                                                       tCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCa                                                                                                                                               >  1:186254/1‑94 (MQ=255)
                                                         aTTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGaaa                                                                                       <  3:79595/148‑1 (MQ=255)
                                                          ttGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGCTCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAAtt                                                                                     >  2:288501/1‑149 (MQ=255)
                                                            ccgcagaaaaaaaagccgaatttctggcaactctaattgcaatgaaattaaaaggtgaaatctgaaagaaATAGCCTGCGTATGGCGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCAAGAAATTCg                                                                                   <  6:186857‑M1/79‑1 (MQ=255)
                                                             ttttttAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCgg                                                                                  >  4:217871/1‑149 (MQ=255)
                                                             ttttttAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCgg                                                                                  <  1:41934/149‑1 (MQ=255)
                                                                   aaCCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCtttt                                                                             <  7:86794/148‑1 (MQ=255)
                                                                    aaaaaagccgaatttctggcaactctaattgcaatgaaaataaaaggtgaaatctgaaagaaATAGCCTGCGTATGGCGCAGGCTATTTCTTTGAGATTTCACCTATTTATTTCATTGCAATTAGAGTAGCCAGAAATTCGGCtttttt                                                                           <  8:143136‑M1/87‑1 (MQ=255)
                                                                      cTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCtttttttt                                                                         <  2:102782/149‑1 (MQ=255)
                                                                      cTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCtttttttt                                                                         >  2:30772/1‑149 (MQ=255)
                                                                      cTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCttttttt                                                                          >  4:90144/1‑148 (MQ=255)
                                                                         aaCTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTc                                                                        <  3:48152/147‑1 (MQ=255)
                                                                         aaCTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTg                                                                      >  8:152409/1‑149 (MQ=255)
                                                                          aCTCAGATTAAAATTCGTTCTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACTCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGc                                                                     <  2:86118/149‑1 (MQ=255)
                                                                           cTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCg                                                                    >  8:256241/1‑149 (MQ=255)
                                                                            cgaatttctggcaactctaattgcaatgaaattaaaaggtgaaatctgaaagaaATAGCCTGCGTATGGCGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCAAGATATTCGGCTTTTTTTTCTGCgg                                                                   <  6:156151‑M1/95‑1 (MQ=255)
                                                                              aGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCgg                                                                   <  3:52257/147‑1 (MQ=255)
                                                                                   aaaaTTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGAt                                                             <  5:219021/148‑1 (MQ=255)
                                                                                         cGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTc                                                       <  7:105276/148‑1 (MQ=255)
                                                                                         cGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTc                                                       <  4:38620/148‑1 (MQ=255)
                                                                                           ttttGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCa                                                                                           <  6:289707/110‑1 (MQ=255)
                                                                                           ttttGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTCCGa                                                    <  8:270214/149‑1 (MQ=255)
                                                                                             ttGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTCCGata                                                  >  4:212564/1‑149 (MQ=255)
                                                                                                ttCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCg                                                                                   >  3:246440/1‑113 (MQ=255)
                                                                                                ttCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCg                                                                                   <  4:246440/113‑1 (MQ=255)
                                                                                                ttCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTAAGATTTCACCTTTTAATTTCATTGCAATTATAGTTGCCACAAATTAGTCTTTTTTTTCTGCGGGCAGATTCTTTCCGATATg                                                >  5:44/1‑148 (MQ=255)
                                                                                                   aGTGAATGATCTTGCCGGAAACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTGTTTTTTCTGCGGGCAGATTCTTTCCGATATGCAc                                             >  6:91122/1‑148 (MQ=255)
                                                                                                    gTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTCCGATATGCAc                                             >  7:22906/1‑147 (MQ=255)
                                                                                                    gTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTCCGATATGCACCa                                           <  5:186294/149‑1 (MQ=255)
                                                                                                     tGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTc                                                       <  2:27266/136‑1 (MQ=255)
                                                                                                     tGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTc                                                       >  1:27266/1‑136 (MQ=255)
                                                                                                       aaTGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTCCGATATGCACCAgg                                         >  3:113263/1‑148 (MQ=255)
                                                                                                       aaTGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTCCGATATGCACCAGGc                                        <  6:19468/149‑1 (MQ=255)
                                                                                                       aaTGATCTTGCCGGAAACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTATTTCAGATTTCACATTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCCGCTTTTTTTTCTGCGGGCAGATTCTTTCCGATATGCACCAgg                                         >  5:101735/1‑148 (MQ=255)
                                                                                                           ctCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTATGCGGGCAGATTCTTTCCGATATGCACCAGGCACAt                                    >  1:19510/2‑149 (MQ=255)
                                                                                                                 ggtgaaatctgaaagaaATAGCCTGCGTATGGCGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTCCGATATGCACCAGGCACAttttttt                              >  7:143136‑M1/18‑149 (MQ=255)
                                                                                                                  cGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTCCGATATGCACCAGGCACAttttttt                              >  8:305048/1‑148 (MQ=255)
                                                                                                                   ggATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTCCGATATGCACCAGGCACATTTTTTTGa                            <  7:152409/149‑1 (MQ=255)
                                                                                                                         acacTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAAtt                                                                                                     >  4:58829/1‑70 (MQ=255)
                                                                                                                         acacTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAAtt                                                                                                     <  3:58829/70‑1 (MQ=255)
                                                                                                                           acTGTTCATAGCCTGCGCCATACGCGGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTCGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTCCGATATGCACCAGGCACATTTTTTTGACACCTtca                    >  3:182244/1‑149 (MQ=255)
                                                                                                                           acTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTCCGATATGCa                                              >  1:261078/1‑123 (MQ=255)
                                                                                                                           acTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTCCGATATGCa                                              <  2:261078/123‑1 (MQ=255)
                                                                                                                           acTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTCCGATATGCACCAGGCACATTTTTTTGACACCTtca                    <  1:30772/149‑1 (MQ=255)
                                                                                                                                 cATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTCCGATATGCACCAGGCGCATTTTTTTGACACCTTCATCAAgt              >  7:22834/1‑149 (MQ=255)
                                                                                                                                  ATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTCCGATATGCACCAGGCACATTTTTTTGACACCTTCATCAAgt              >  3:217871‑M1/1‑148 (MQ=255)
                                                                                                                                  ATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCCTTCCGATATGCACCAGGCACATTTTTTTGACACCTTCATCAAgt              >  3:90144‑M1/1‑148 (MQ=255)
                                                                                                                                        TGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTCCGATATGCACCAGGCACATTTTTTTGACACCTTCATCAAGTGTTTTTa       <  8:136076‑M1/149‑1 (MQ=255)
                                                                                                                                             cATACGCAGGCTATTTATTTCAGATTTCACGTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTGTTCTGCGGGCAGATTCTTTCAGATATGCACCAGGCACAGTTATTTGACACCTTCATCAAGTGTTTTTACGTTg  >  5:227972/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
TCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTCCGATATGCACCAGGCACATTTTTTTGACACCTTCATCAAGTGTTTTTACGTTG  >  NC_000913/1417328‑1417617

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