Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,219,633 C→A 7.0% intergenic (‑152/+59) 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,6330CA 7.0% 261.0 / 6.3 80intergenic (‑152/+59)osmF/bglXputative ABC transporter periplasmic binding protein/beta‑D‑glucoside glucohydrolase, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (38/36);  new base A (4/2);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-01

CCTCATTCTTTTACTGTTGTTTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCATATTTCAT  >  NC_000913/2219490‑2219781
                                                                                                                                               |                                                                                                                                                    
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                                                                                         >  7:172614/1‑149 (MQ=255)
                                                        gCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTa                                                                                         <  1:355723/149‑1 (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                                                                          >  8:55999/1‑141 (MQ=255)
                                                                               aTTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTc                                                                          <  7:55999/141‑1 (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                                                                                                                        <  5:48209/70‑1 (MQ=255)
                                                                                                        cTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAgg                                                                                                                        >  6:48209/1‑70 (MQ=255)
                                                                                                         ttAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGAt                                         >  7:409853/1‑148 (MQ=255)
                                                                                                              agGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGc                                                                                             >  2:117290/1‑91 (MQ=255)
                                                                                                              agGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGc                                                                                             <  1:117290/91‑1 (MQ=255)
                                                                                                                   tatTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTg                               >  7:246457/1‑148 (MQ=255)
                                                                                                                   tatTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTg                               >  6:169669/1‑148 (MQ=255)
                                                                                                                   tatTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTg                               >  1:323800/1‑148 (MQ=255)
                                                                                                                     ttCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGc                                                                                     >  2:321993/1‑92 (MQ=255)
                                                                                                                     ttCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGc                                                                                     <  1:321993/92‑1 (MQ=255)
                                                                                                                        ggTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCt                                                                  >  1:30012/1‑108 (MQ=255)
                                                                                                                        ggTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCt                                                                  <  2:30012/108‑1 (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                        <  5:337705/136‑1 (MQ=255)
                                                                                                                                      tGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGc                        >  6:337705/1‑136 (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)
                                                                                                                                               cggcgCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTATCGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCATATTTCAt  <  6:333448/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CCTCATTCTTTTACTGTTGTTTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCATATTTCAT  >  NC_000913/2219490‑2219781

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