Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,218,893 C→T 6.0% G197S (GGC→AGC)  osmF ← putative ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,218,8930CT 6.0% 215.4 / 3.4 67G197S (GGC→AGC) osmFputative ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (28/35);  new base T (2/2);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.25e-01

GCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATCGGCGCGTTCGATAAACTCTGCCGAGGCTGCCAGTTTGAAGGTGCCGCCCTCTTGCAGATAACGACTCAGGTCGGCA  >  NC_000913/2218751‑2219041
                                                                                                                                              |                                                                                                                                                    
gcgTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCAccagc                                                                                                                                                >  4:55335/1‑149 (MQ=255)
gcgTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCAccagc                                                                                                                                                >  2:226938/1‑149 (MQ=255)
     gATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTg                                                                                                                                           >  2:149996/1‑149 (MQ=255)
      aTAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCTACCAGCCAGTGa                                                                                                                                          >  1:416735/1‑149 (MQ=255)
      aTAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCTACCAGCCAGTGa                                                                                                                                          >  1:405976/1‑149 (MQ=255)
      aTAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGa                                                                                                                                          >  4:532374/1‑149 (MQ=255)
         ggTTGCACACCTTGCGGTTCGCTTAAGGTTTGCAGCCCCAGCGCCGCTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAg                                                                                                                                       <  1:247898/149‑1 (MQ=255)
          gTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGc                                                                                                                                      >  1:94631/1‑149 (MQ=255)
          gTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGc                                                                                                                                      >  1:454710/1‑149 (MQ=255)
           ttGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCa                                                                                                                                     <  3:305416/149‑1 (MQ=255)
            tGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCaa                                                                                                                                    <  1:234417/149‑1 (MQ=255)
             gCACACCTTGCGGATCGCTTAAGGGTTGCAGCCCCAGCGCCGCTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAAc                                                                                                                                   <  1:12190/149‑1 (MQ=255)
               acacCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTg                                                                                                                                 <  3:219276/149‑1 (MQ=255)
                cacCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGa                                                                                                                                >  3:416656/1‑149 (MQ=255)
                 acCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGAt                                                                                                                               >  1:210741/1‑149 (MQ=255)
                   cTTGCGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATcc                                                                                                                             >  3:98007/1‑149 (MQ=255)
                       cGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGAc                                                                                                                         <  1:503355/149‑1 (MQ=255)
                             gCTTAAGGTTTGCAGCCCCAGCGCCGCTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCt                                                                                                                   <  1:431150/149‑1 (MQ=255)
                              cTTAAGGTTTGCAGCCCCAGCGCCGCTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCtt                                                                                                                  >  3:174105/1‑149 (MQ=255)
                                 aaGGTTTGCAGCCCCAGCGCCGCTACCGGGCCGTCAGTGCTGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTaaa                                                                                                               >  4:117072/1‑149 (MQ=255)
                                    gTTTGCAGCCCCAGCGCCGCTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGcc                                                                                                            >  1:41329/1‑149 (MQ=255)
                                         cAGCCCCAGCGCCGCTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAgg                                                                                                       >  2:182733/1‑149 (MQ=255)
                                         cAGCCCCAGCGCCGCTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAgg                                                                                                       <  3:55335/149‑1 (MQ=255)
                                              ccAGCGCCGCTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCtttt                                                                                                  >  3:27036/1‑149 (MQ=255)
                                              ccAGCGCCGCTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCtttt                                                                                                  <  1:266355/149‑1 (MQ=255)
                                                   gccgcTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCaaa                                                                                             <  3:12606/149‑1 (MQ=255)
                                                        tACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCgg                                                                                        <  1:149132/149‑1 (MQ=255)
                                                        tACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCgg                                                                                        >  3:553148/1‑149 (MQ=255)
                                                        tACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCgg                                                                                        <  1:341081/149‑1 (MQ=255)
                                                         aCCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCggg                                                                                       >  2:185967/1‑149 (MQ=255)
                                                         aCCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCggg                                                                                       <  4:323616/149‑1 (MQ=255)
                                                         aCCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCggg                                                                                       <  4:90820/149‑1 (MQ=255)
                                                         aCCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCggg                                                                                       <  2:560569/149‑1 (MQ=255)
                                                          ccGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGt                                                                                      <  4:98007/149‑1 (MQ=255)
                                                           cGGGCCGTCAGTGCCGTAAGCCATTGCAGTATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCTGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTa                                                                                     >  2:280808/1‑149 (MQ=255)
                                                           cGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCTGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTa                                                                                     >  4:248249/1‑149 (MQ=255)
                                                            gggCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTaa                                                                                    >  4:565786/1‑149 (MQ=255)
                                                            gggCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTaa                                                                                    >  4:62293/1‑149 (MQ=255)
                                                             ggCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCTACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAAc                                                                                   <  2:405976/149‑1 (MQ=255)
                                                             ggCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCTACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAAc                                                                                   <  2:416735/149‑1 (MQ=255)
                                                             ggCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAAc                                                                                   <  2:210741/149‑1 (MQ=255)
                                                                            aaGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATCGGCGCGTTCg                                                                    >  1:45420/1‑149 (MQ=255)
                                                                               ccATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATCGGCGCGTTCGATa                                                                 <  1:525698/149‑1 (MQ=255)
                                                                                 aTTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATCGGCGCGTTCGATaaa                                                               <  1:182733/149‑1 (MQ=255)
                                                                                  ttGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATc                                                                              >  1:156763/1‑133 (MQ=255)
                                                                                  ttGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATc                                                                              <  2:156763/133‑1 (MQ=255)
                                                                                  ttGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATCGGCGCGTTCGATAAAc                                                              <  2:8523/149‑1 (MQ=255)
                                                                                     cagcaTTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATCGGCGCGTTCGATAAActct                                                           <  1:226938/149‑1 (MQ=255)
                                                                                     cagcaTTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATCGGCGCGTTCGATAAActct                                                           <  2:163618/149‑1 (MQ=255)
                                                                                       gcaTTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATCGGCGCGTTCGATAAAc                                                              <  2:409681/144‑1 (MQ=255)
                                                                                       gcaTTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATCGGCGCGTTCGATAAAc                                                              <  4:170098/144‑1 (MQ=255)
                                                                                       gcaTTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATCGGCGCGTTCGATAAAc                                                              >  1:409681/1‑144 (MQ=255)
                                                                                       gcaTTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATCGGCGCGTTCGATAAAc                                                              >  3:170098/1‑144 (MQ=255)
                                                                                       gcaTTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATCGGCGCGTTCGATAAACTCTGc                                                         <  3:117072/149‑1 (MQ=255)
                                                                                          ttAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATCGGCGCGTTCGATAAACTCTGCCGa                                                      <  3:62293/149‑1 (MQ=255)
                                                                                                  agagGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATCGGCGCGTTCGATAAACTCTGCCGAGGCTGCCa                                              <  4:174105/149‑1 (MQ=255)
                                                                                                      gTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATCGGCGCGTTCGATAAACTCTGCCGAGGCTGCCAGttt                                          >  3:565603/1‑149 (MQ=255)
                                                                                                        ttGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATCGGCGCGTTCGATAAACTCTGCCGAGGCTGCCAGTTTGa                                        <  1:359308/149‑1 (MQ=255)
                                                                                                         tGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGATTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATCGGCGCGTTCGATAAACTCTGCCGAGGCTGCCAGTTTGaa                                       <  2:364568/149‑1 (MQ=255)
                                                                                                            tGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATCGGCGCGTTCGATAAACTCTGCCGAGGCTGCCAGTTTGAAGGt                                    >  3:221487/1‑149 (MQ=255)
                                                                                                                      cTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATCGGCGCGTTCGATAAACTCTGCCGAGGCTGCCAGTTTGAAGGTGCCGCCctct                          <  3:8447/149‑1 (MQ=255)
                                                                                                                          gATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAACGCCATAGGCTTTTTCAACCGCGGGTAACGCATCGGCGCGTTCGATAAACTCTGCCGAGGCTGCCAGTTTGAAGGTGCCGCCCTCTTGCa                      <  4:78266/149‑1 (MQ=255)
                                                                                                                              gTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATCGGCGTGTTCGATAAACTCTGCCGAGGCTGCCAGTTTGAAGGTGCCGCCCTCTTGCAGATa                  <  1:185967/149‑1 (MQ=255)
                                                                                                                               tCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATCGGCGCGTTCGATAAACTCTGCCGAGGCTGCCAGTTTGAAGGTGCCGCCCTCTTGCAGATaa                 <  2:359426/149‑1 (MQ=255)
                                                                                                                                  ccgccgTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATCGGCGCGTTCGATAAACTCTGCCGAGGCTGCCAGTTTGAAGGTGCCGCCCTCTTGCAGATAACGa              <  2:45420/149‑1 (MQ=255)
                                                                                                                                          tCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATCGGCGCGTTCGATAAACTCTGCCGAGGCTGCCAGTTTGAAGGTGCCGCCCTCTTGCAGATAACGACTCAGGTc      >  4:187821/1‑149 (MQ=255)
                                                                                                                                            gCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATCGGCGCGTTCGATAAACTCTGCCGAGGCTGCCAGTTTGAAGGTGCCGCCCTCTTGCAGATAACGACTCAGGTCgg    >  2:261406/1‑149 (MQ=255)
                                                                                                                                              caccaGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATCGGCGCGTTCGATAAACTCTGCCGAGGCTGCCAGTTTGAAGGTGCCGCCCTCTTGCAGATAACGACTCAGGTCGGCa  >  1:231021/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCCAGCGCCGCTACCGGGCCGTCAGTGCCGTAAGCCATTGCAGCATTAACGCCAGAGGTTTGCTGGGCAGCGGCTTTGATCGTCACCGCCGTGTCGCCACCAGCCAGTGACAGCAACTGATCCTGACCGAGCTTAAAGCCATAGGCTTTTTCAAACGCGGGTAACGCATCGGCGCGTTCGATAAACTCTGCCGAGGCTGCCAGTTTGAAGGTGCCGCCCTCTTGCAGATAACGACTCAGGTCGGCA  >  NC_000913/2218751‑2219041

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