Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 251,639 A→G 5.5% T248A (ACG→GCG)  dinB → DNA polymerase IV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913251,6390AG 5.5% 275.6 / 3.7 89T248A (ACG→GCG) dinBDNA polymerase IV
Reads supporting (aligned to +/- strand):  ref base A (39/45);  new base G (3/2);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

GACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAGC  >  NC_000913/251491‑251786
                                                                                                                                                    |                                                                                                                                                   
gACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGGGGAACGCa                                                                                                                                                     <  3:269886/149‑1 (MQ=255)
 aCCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCAc                                                                                                                                                    <  5:116181/149‑1 (MQ=255)
 aCCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAACTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTGGCGAAAATCCGTCGGCGTGGAACGCAc                                                                                                                                                    <  7:73342/149‑1 (MQ=255)
     gCGGGGATGTAGAATAGTGTGATCGGGTGCGGCTGCTTGGACGGTTTGGGGTGTTTGGGCGCTGTTTGTGGGGGCGTTGTGGGGGGATGGGCGATCGCGATGGTGGCCGCGCGCGGGGGCGAAAATCCGTCGGCGTGGCGCGCGCGATg                                                                                                                                                >  4:160680/1‑149 (MQ=255)
      cGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAATCGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATg                                                                                                                                                <  8:361088/148‑1 (MQ=255)
             gTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCgg                                                                                                                                            >  4:21616/1‑145 (MQ=255)
              tACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGa                                                                                                                                        >  6:287653/1‑148 (MQ=255)
                cAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTACCAGCGCACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGatat                                                                                                                                     >  2:44816/1‑149 (MQ=255)
                 aaaaGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATAtt                                                                                                                                    >  5:350721/1‑149 (MQ=255)
                    aGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATtcat                                                                                                                                 <  3:5401/149‑1 (MQ=255)
                      tgtgATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATtcatca                                                                                                                               >  8:103394/1‑149 (MQ=255)
                      tgtgATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATtcatca                                                                                                                               <  8:253874/149‑1 (MQ=255)
                        tgATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCAc                                                                                                                              <  5:287653/148‑1 (MQ=255)
                          aTCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTgg                                                                                                                           >  6:244495/1‑149 (MQ=255)
                           tCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGt                                                                                                                          >  1:263826/1‑149 (MQ=255)
                                  aTGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCTAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAAt                                                                                                                    <  5:132513/148‑1 (MQ=255)
                                  aTGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAAt                                                                                                                    <  5:293932/148‑1 (MQ=255)
                                    gctgctTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGc                                                                                                                                              <  3:146264/120‑1 (MQ=255)
                                    gctgctTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGc                                                                                                                                              >  4:146264/1‑120 (MQ=255)
                                     ctgctTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAAtgtg                                                                                                                 >  3:130483/1‑148 (MQ=255)
                                     ctgctTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAAAGCACGATGGCGGAAGATATTCATCACTGGTCTGAAtgtg                                                                                                                 <  3:226856/148‑1 (MQ=255)
                                     ctgctTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGAGATGTTAACAGCGAACTGTTGAGAAAATCCTTAGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAAtgtg                                                                                                                 <  1:65854/148‑1 (MQ=255)
                                      tgctTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAAtgtg                                                                                                                 >  7:165830/1‑147 (MQ=255)
                                      tgctTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAAtgtg                                                                                                                 >  4:168254/1‑147 (MQ=255)
                                       gctTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAg                                                                                                              >  2:35487/1‑149 (MQ=255)
                                          tAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCGCGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGa                                                                                                           >  7:192911/1‑149 (MQ=255)
                                          tAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCGCGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGa                                                                                                           >  7:187187/1‑149 (MQ=255)
                                           aaaCGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGa                                                                                                           <  5:244495/148‑1 (MQ=255)
                                           aaaCGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGa                                                                                                           <  5:110900/148‑1 (MQ=255)
                                              cGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCGCGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTAt                                                                                                       <  8:187187/149‑1 (MQ=255)
                                              cGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCGCGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTAt                                                                                                       <  8:192911/149‑1 (MQ=255)
                                                     gCAAAATTGGCTGCAGTTTGTGGGAGCGGAGTCAGGGGATGGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGGCTGAATGTGAAGCGATTATCGAGCg                                                                                                 <  5:314611/148‑1 (MQ=255)
                                                      cAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAAt                                                                                                                    <  3:84539/128‑1 (MQ=255)
                                                      cAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAAt                                                                                                                    >  4:84539/1‑128 (MQ=255)
                                                      cAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGc                                                                                               >  1:425376/1‑149 (MQ=255)
                                                        aaTTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTg                                                                                             >  7:24357/1‑149 (MQ=255)
                                                               ccGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACTCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCg                                                                                       >  5:175723/1‑148 (MQ=255)
                                                               ccGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGaa                                                                                                                                          >  2:66449/1‑97 (MQ=255)
                                                               ccGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGaa                                                                                                                                          <  1:66449/97‑1 (MQ=255)
                                                               ccGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCg                                                                                       >  2:153226/1‑148 (MQ=255)
                                                               cagCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGAGAACGGTTGCGAAAAACCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTATGAATCTGAAGCGATAATCGAGCGGCTGTATCCg                                                                                       >  2:278089/3‑148 (MQ=255)
                                                                cGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCgg                                                                                      >  3:227178/1‑148 (MQ=255)
                                                                 gCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGaa                                                                                    <  3:21616/149‑1 (MQ=255)
                                                                    ttttGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCt                                                                                              <  2:219716/136‑1 (MQ=255)
                                                                    ttttGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCt                                                                                              >  1:219716/1‑136 (MQ=255)
                                                                    ttttGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACtt                                                                                 >  6:146599/1‑149 (MQ=255)
                                                                     tttGTGGGAGCGTAGTCAGGGGATTGACGAACTCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGAAGGATGGCGGTAGATATTCATCACTGGTGTGAATGTGAAGCGATTATCTAGCGGCTGTATCCGGAACtt                                                                                 <  2:45348/148‑1 (MQ=255)
                                                                         tgGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACg                                                                            <  6:175723/149‑1 (MQ=255)
                                                                         tgGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACg                                                                            <  1:35487/149‑1 (MQ=255)
                                                                                      aGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTg                                                                     >  7:291572/1‑143 (MQ=255)
                                                                                      aGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTg                                                                     <  8:291572/143‑1 (MQ=255)
                                                                                      aGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCaaa                                                                <  1:153226/148‑1 (MQ=255)
                                                                                      aGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAg                                                               <  3:168254/149‑1 (MQ=255)
                                                                                        gggATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGACCGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGCAGCGATTATCGAGCGGCTGTATCCGGAACTTGACCGCCGGCTGGCAAAgg                                                              >  4:306652/1‑148 (MQ=255)
                                                                                          gATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATtc                                                                                                                                   >  7:420243/1‑77 (MQ=255)
                                                                                          gATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATtc                                                                                                                                   <  8:420243/77‑1 (MQ=255)
                                                                                             tGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATtc                                                                                                                                   >  4:251526/1‑74 (MQ=255)
                                                                                             tGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATtc                                                                                                                                   <  3:251526/74‑1 (MQ=255)
                                                                                              gACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAAcc                                                       >  2:399122/1‑149 (MQ=255)
                                                                                               aCGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACg                                                                            >  4:17860/1‑127 (MQ=255)
                                                                                               aCGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACg                                                                            <  3:17860/127‑1 (MQ=255)
                                                                                                cGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCt                                                      >  1:425444/1‑148 (MQ=255)
                                                                                                 gAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGa                                                    <  4:227178/149‑1 (MQ=255)
                                                                                                  aaCGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGAt                                                   >  1:213005/1‑149 (MQ=255)
                                                                                                    cgcgATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGAtt                                                  >  4:124727/1‑148 (MQ=255)
                                                                                                             aaCAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTg                                         <  2:425376/148‑1 (MQ=255)
                                                                                                              aCAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCt                                       >  3:40986/1‑149 (MQ=255)
                                                                                                              aCAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACCTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCt                                       >  5:208792/1‑149 (MQ=255)
                                                                                                               cAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTc                                      <  8:165830/149‑1 (MQ=255)
                                                                                                                aGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGc                                                                                               <  5:350846/91‑1 (MQ=255)
                                                                                                                aGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGc                                                                                               >  6:350846/1‑91 (MQ=255)
                                                                                                                 gCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGc                                    <  3:160680/149‑1 (MQ=255)
                                                                                                                       ggTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACCTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAgggg                              <  6:208792/149‑1 (MQ=255)
                                                                                                                           gCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGa                                            >  3:256243/1‑131 (MQ=255)
                                                                                                                           gCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGa                                            <  4:256243/131‑1 (MQ=255)
                                                                                                                            cGAATATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTTGTCTGATTGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGaa                         <  6:117619/149‑1 (MQ=255)
                                                                                                                             gAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGaaa                        <  7:103394/149‑1 (MQ=255)
                                                                                                                             gAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGaaa                        >  8:323882/1‑149 (MQ=255)
                                                                                                                               aaaTCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAAtt                      <  8:247182/149‑1 (MQ=255)
                                                                                                                                  tCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTaaa                   <  2:213005/149‑1 (MQ=255)
                                                                                                                                  tCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGCTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTaaa                   <  1:278089/149‑1 (MQ=255)
                                                                                                                                    cGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGt                 <  1:399122/149‑1 (MQ=255)
                                                                                                                                     gTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGaaa                        >  5:253491/1‑141 (MQ=255)
                                                                                                                                     gTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGaaa                        <  6:253491/141‑1 (MQ=255)
                                                                                                                                          cGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACtt                                                                                 >  7:407911/1‑79 (MQ=255)
                                                                                                                                          cGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACtt                                                                                 <  8:407911/79‑1 (MQ=255)
                                                                                                                                              gAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAATGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGAtt        <  5:146599/148‑1 (MQ=255)
                                                                                                                                                 cGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGAtttt      >  7:53805/1‑147 (MQ=255)
                                                                                                                                                  gCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAg   <  6:350721/149‑1 (MQ=255)
                                                                                                                                                  gCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAg   <  2:263826/149‑1 (MQ=255)
                                                                                                                                                    aCGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAGc  >  1:319932/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAGC  >  NC_000913/251491‑251786

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