Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,219,632 G→A 6.2% intergenic (‑151/+60) osmF ← / ← bglX putative ABC transporter periplasmic binding protein/beta‑D‑glucoside glucohydrolase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,219,6320GA 6.2% 251.3 / 4.5 80intergenic (‑151/+60)osmF/bglXputative ABC transporter periplasmic binding protein/beta‑D‑glucoside glucohydrolase, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (39/36);  new base A (3/2);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.81e-01

CCTCATTCTTTTACTGTTGTTTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCATATTTC  >  NC_000913/2219490‑2219779
                                                                                                                                              |                                                                                                                                                   
ccTCATTCTTTTACTGTTGTTTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATAAGgcgc                                                                                                                                               >  3:375167/1‑149 (MQ=255)
  tCATTCTTTTACTGTTGTTTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGgcgc                                                                                                                                               >  3:64455/1‑147 (MQ=255)
    aTTCTTTTACTGTTGTTTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGa                                                                                                                                             <  5:398952/147‑1 (MQ=255)
    aTTCTTTTACTGTTGTTTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAgc                                                                                                                                           >  4:272461/1‑149 (MQ=255)
         tttACTGTTGTTTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGcc                                                                                                                                        >  6:57692/1‑147 (MQ=255)
         tttACTGTTGTTTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGcc                                                                                                                                        <  5:57692/147‑1 (MQ=255)
          ttACTGTTGTTTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCtt                                                                                                                                      <  6:406318/148‑1 (MQ=255)
            aCTGTTGTTTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATc                                                                                                                                   >  8:254654/1‑149 (MQ=255)
             cTGTTGTTTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATcc                                                                                                                                  >  5:356109/1‑149 (MQ=255)
              tgttgtTTTCAGCGAAATAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAAACAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATccc                                                                                                                                 >  5:336590/1‑148 (MQ=255)
                ttgttTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATAAGGCGCTCGCGCCGCATCCgg                                                                                                                                <  4:375167/148‑1 (MQ=255)
                ttgttTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGAATCAGTTGCAGAATCAGATAAATTCTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGc                                                                                                                               <  6:113853/149‑1 (MQ=255)
                   ttttCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTa                                                                                                                            <  3:272461/149‑1 (MQ=255)
                   ttttCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTGGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTa                                                                                                                            <  3:105777/149‑1 (MQ=255)
                     ttCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTAc                                                                                                                           <  8:145550/148‑1 (MQ=255)
                        aGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACaaa                                                                                                                        <  4:124013/148‑1 (MQ=255)
                         gCGAATTAAGATAAAAGCATAGTTGATAATGGCGGTGTTAGATTGAGCTTAAGGATTCAGTACCAGAATAAGATAATTACTTTAGATCCATGTTCTGTCATGCGTTTACTGCAGGATGAGGCGGGAGAGACTTATCAGGACTACATAgc                                                                                                                      >  2:342203/1‑148 (MQ=255)
                         gCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAgg                                                                                                                      <  7:254654/149‑1 (MQ=255)
                               tAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCa                                                                                                                 >  5:333448/1‑148 (MQ=255)
                                      aaaGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTcgc                                                                                                         <  6:226260/149‑1 (MQ=255)
                                       aaGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTcgcg                                                                                                        >  3:209946/1‑149 (MQ=255)
                                       aaGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGCATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGTTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTcgcg                                                                                                        >  5:235235/1‑149 (MQ=255)
                                         gCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGcc                                                                                                      >  3:152788/1‑149 (MQ=255)
                                           aTAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCt                                                                                                    <  1:298813/149‑1 (MQ=255)
                                           aTAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCt                                                                                                    <  3:35908/149‑1 (MQ=255)
                                               ttGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACt                                                                                                <  7:15411/149‑1 (MQ=255)
                                                 gATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTa                                                                                               >  7:87048/1‑148 (MQ=255)
                                                 gATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTa                                                                                               <  4:64455/148‑1 (MQ=255)
                                                  aTAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCaaa                                                                                                               <  1:125793/131‑1 (MQ=255)
                                                  aTAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCaaa                                                                                                               >  2:125793/1‑131 (MQ=255)
                                                   tAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAg                                                                                            <  4:209946/149‑1 (MQ=255)
                                                     aTGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGTCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCa                                                                                          <  8:431957/149‑1 (MQ=255)
                                                      tGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCAt                                                                                         >  3:275602/1‑149 (MQ=255)
                                                        gCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTa                                                                                       <  1:355723/149‑1 (MQ=255)
                                                        gCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTa                                                                                       >  7:172614/1‑149 (MQ=255)
                                                             gTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGc                                                                                   >  2:222761/1‑148 (MQ=255)
                                                              ttcgCTTGATCGTAAGTATTCAGTTGCAGACTCAGATAAATACTTAAGAGGCATAGTCGGTCAGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCaa                                                                                 <  4:384671/146‑1 (MQ=38)
                                                                   ttGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCa                                                                             >  8:149652/1‑148 (MQ=255)
                                                                        cGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCg                                                                       <  5:66915/149‑1 (MQ=255)
                                                                              gATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTc                                                                 <  4:275602/149‑1 (MQ=255)
                                                                               aTTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTc                                                                        <  7:55999/141‑1 (MQ=255)
                                                                               aTTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTc                                                                        >  8:55999/1‑141 (MQ=255)
                                                                               aTTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCt                                                                >  2:383227/1‑149 (MQ=255)
                                                                                ttCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGc                                                                      <  8:190237/142‑1 (MQ=255)
                                                                                ttCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGc                                                                      >  7:190237/1‑142 (MQ=255)
                                                                                    gTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATgcg                                                                                                                                                   <  3:141557/61‑1 (MQ=255)
                                                                                    gTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATgcg                                                                                                                                                   >  4:141557/1‑61 (MQ=255)
                                                                                         aGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGc                                                      <  2:418987/149‑1 (MQ=255)
                                                                                         aGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGc                                                      >  3:33715/1‑149 (MQ=255)
                                                                                               aGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAAt                                                 <  8:159691/148‑1 (MQ=255)
                                                                                               aGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAAt                                                 >  2:74472/1‑148 (MQ=255)
                                                                                                  tAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGt                                             >  2:214525/1‑149 (MQ=255)
                                                                                                      tACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCa                                                                             >  1:134950/1‑113 (MQ=255)
                                                                                                      tACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCa                                                                             <  2:134950/113‑1 (MQ=255)
                                                                                                        cTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAgg                                                                                                                      >  6:48209/1‑70 (MQ=255)
                                                                                                        cTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAgg                                                                                                                      <  5:48209/70‑1 (MQ=255)
                                                                                                         ttAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGAt                                       >  7:409853/1‑148 (MQ=255)
                                                                                                              agGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGc                                                                                           >  2:117290/1‑91 (MQ=255)
                                                                                                              agGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGc                                                                                           <  1:117290/91‑1 (MQ=255)
                                                                                                                   tatTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTg                             >  1:323800/1‑148 (MQ=255)
                                                                                                                   tatTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTg                             >  6:169669/1‑148 (MQ=255)
                                                                                                                   tatTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTg                             >  7:246457/1‑148 (MQ=255)
                                                                                                                     ttCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGc                                                                                   <  1:321993/92‑1 (MQ=255)
                                                                                                                     ttCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGc                                                                                   >  2:321993/1‑92 (MQ=255)
                                                                                                                        ggTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCt                                                                <  2:30012/108‑1 (MQ=255)
                                                                                                                        ggTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCt                                                                >  1:30012/1‑108 (MQ=255)
                                                                                                                        ggTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACAt                               <  1:371027/141‑1 (MQ=255)
                                                                                                                        ggTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACAt                               >  2:371027/1‑141 (MQ=255)
                                                                                                                         gTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGc                      >  3:382365/1‑149 (MQ=255)
                                                                                                                         gTCCGGCGTTTCCTGCCGGAGGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGc                      <  6:356109/149‑1 (MQ=255)
                                                                                                                           ccGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCa                    <  1:383227/149‑1 (MQ=255)
                                                                                                                             ggCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCGTTTCTTAACGCGTGCGGAATCAGTGCCGATTAAGACATTGAACTTGCCAgg                  >  2:54320/1‑149 (MQ=255)
                                                                                                                                      tGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGc                      >  6:337705/1‑136 (MQ=255)
                                                                                                                                      tGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGc                      <  5:337705/136‑1 (MQ=255)
                                                                                                                                      tGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGt          >  4:390052/1‑148 (MQ=255)
                                                                                                                                       gCCGGATAAGGCGCTCGCGCCGCATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCACCTCAACCTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCa        >  5:104734/1‑149 (MQ=255)
                                                                                                                                         cGGATAAGGCGCTCGCGCCGCATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTaa                                                             >  7:68629/1‑94 (MQ=37)
                                                                                                                                         cGGATAAGGCGCTCGCGCCGCATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTaa                                                             <  8:68629/94‑1 (MQ=37)
                                                                                                                                          ggATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCATGCTCGGCGTCatat     <  6:235235/149‑1 (MQ=255)
                                                                                                                                          ggATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAAATCAACCTCGCCTTTCTTAACGCGTGCGGAATCAGTTCCGATAAAGACATTGAACTTGACATGCTCGGCGTCatat     >  4:278608/1‑149 (MQ=255)
                                                                                                                                           gATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCATAtt    >  4:352347/1‑149 (MQ=255)
                                                                                                                                             tGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTCCAGCATCTCAAACTCGCCTTTGTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCATATTTc  <  6:336590/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
CCTCATTCTTTTACTGTTGTTTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCATATTTC  >  NC_000913/2219490‑2219779

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: