Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,417,468 C→T 7.2% intergenic (‑82/+20) 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,4680CT 7.2% 259.9 / 7.1 83intergenic (‑82/+20)recE/racCRac prophage; exonuclease VIII, 5' to 3' specific dsDNA exonuclease/Rac prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (40/37);  new base T (2/4);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

GGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTCCGATATGCACCAGGCACATTTTTTTGACACCTTCATCAAGTGTTTTTA  >  NC_000913/1417320‑1417612
                                                                                                                                                    |                                                                                                                                                
ggACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTgcgc                                                                                                                                                  >  4:52257/1‑149 (MQ=255)
ggACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTgcgc                                                                                                                                                  <  2:259440/149‑1 (MQ=255)
        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                                                                                                                         >  2:166415/1‑149 (MQ=255)
                         ttttttCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGAt                                                                                                                         >  3:72307/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                                                                                                      >  2:256302/1‑149 (MQ=255)
                                            aGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCAtt                                                                                                      >  8:95191/1‑149 (MQ=255)
                                            aGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCAtt                                                                                                      <  1:166415/149‑1 (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                                                                                                                                    >  7:11029/1‑101 (MQ=255)
                                                              cTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGTATGGCGCAGGCTAt                                                                                                                                    <  8:11029/101‑1 (MQ=255)
                                                               tCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCa                                                                                                                                          >  1:186254/1‑94 (MQ=255)
                                                               tCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCa                                                                                                                                          <  2:186254/94‑1 (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                                                                             <  1:41934/149‑1 (MQ=255)
                                                                     ttttttAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCgg                                                                             >  4:217871/1‑149 (MQ=255)
                                                                           aaCCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCtttt                                                                        <  7:86794/148‑1 (MQ=255)
                                                                            aaaaaagccgaatttctggcaactctaattgcaatgaaaataaaaggtgaaatctgaaagaaATAGCCTGCGTATGGCGCAGGCTATTTCTTTGAGATTTCACCTATTTATTTCATTGCAATTAGAGTAGCCAGAAATTCGGCtttttt                                                                      <  8:143136‑M1/87‑1 (MQ=255)
                                                                              cTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCtttttttt                                                                    >  2:30772/1‑149 (MQ=255)
                                                                              cTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCtttttttt                                                                    <  2:102782/149‑1 (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                                                  <  4:38620/148‑1 (MQ=255)
                                                                                                 cGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTc                                                  <  7:105276/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                                                                                                <  3:58829/70‑1 (MQ=255)
                                                                                                                                 acacTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAAtt                                                                                                >  4:58829/1‑70 (MQ=255)
                                                                                                                                   acTGTTCATAGCCTGCGCCATACGCGGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTCGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTCCGATATGCACCAGGCACATTTTTTTGACACCTtca               >  3:182244/1‑149 (MQ=255)
                                                                                                                                   acTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTCCGATATGCa                                         <  2:261078/123‑1 (MQ=255)
                                                                                                                                   acTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTCCGATATGCa                                         >  1:261078/1‑123 (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)
                                                                                                                                                    |                                                                                                                                                
GGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGCAGGCTATTTCTTTCAGATTTCACCTTTTAATTTCATTGCAATTAGAGTTGCCAGAAATTCGGCTTTTTTTTCTGCGGGCAGATTCTTTCCGATATGCACCAGGCACATTTTTTTGACACCTTCATCAAGTGTTTTTA  >  NC_000913/1417320‑1417612

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: