Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,504,968 Δ1 bp 5.7% coding (1113/1176 nt) ydcO ← BenE family inner membrane putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,504,9680T. 5.7% 231.9 / 6.5 70coding (1113/1176 nt)ydcOBenE family inner membrane putative transporter
Reads supporting (aligned to +/- strand):  ref base T (32/34);  new base . (2/2);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

GTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCAGCGCCTGATACAAACTGCCGCCGATGGTACTTAACAGCGCCAGACCTGC  >  NC_000913/1504823‑1505111
                                                                                                                                                 |                                                                                                                                               
gTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTTAACGTAGCAGGGATCAAAGTACTTCATGGCGTGGTGGATGTCAATTTCGGTTTGTGTCAGCGATTAAATTACACTCCACGTAACAAACGCCTCCGGCAATTAATccc                                                                                                                                              >  3:217369/1‑149 (MQ=255)
gTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATccc                                                                                                                                              >  5:151780/1‑149 (MQ=255)
         ggAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCa                                                                                                                                     <  8:67115/149‑1 (MQ=255)
            aCCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGaa                                                                                                                                  >  8:68698/1‑149 (MQ=255)
            aCCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGaa                                                                                                                                  <  5:88547/149‑1 (MQ=255)
            aCCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGaa                                                                                                                                  <  4:60800/149‑1 (MQ=255)
                aCCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCaa                                                                                                                              >  4:189179/1‑149 (MQ=255)
                   gAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAACCCCAAAACGCAGAACCAATccc                                                                                                                          >  1:358632/1‑149 (MQ=255)
                   gAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAACCCCAAAACGCAGAACCAATccc                                                                                                                          >  1:363960/1‑149 (MQ=255)
                    aaaTGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATccc                                                                                                                          <  4:76215/149‑1 (MQ=255)
                            cAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGc                                                                                                                   <  6:151780/148‑1 (MQ=255)
                             aGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGt                                                                                                                 <  8:126411/149‑1 (MQ=255)
                              ggCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTc                                                                                                                <  2:365690/149‑1 (MQ=255)
                                       cTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACtt                                                                                                       >  4:118628/1‑149 (MQ=255)
                                            cGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGt                                                                                                  <  1:36533/149‑1 (MQ=255)
                                             gTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGt                                                                                                  >  8:340298/1‑148 (MQ=255)
                                               aGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTAc                                                                                               >  7:162639/1‑149 (MQ=255)
                                                 cAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTAc                                                                                               >  5:122400/1‑147 (MQ=255)
                                                  aGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAg                                                                                            <  5:39674/149‑1 (MQ=255)
                                                      aTCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAAt                                                                                        >  8:320378/1‑149 (MQ=255)
                                                         cACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGcc                                                                                     >  1:165931/1‑149 (MQ=255)
                                                         cACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAACCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGcca                                                                                    <  2:363960/149‑1 (MQ=255)
                                                         cACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAACCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGcca                                                                                    <  2:358632/149‑1 (MQ=255)
                                                          aCGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGcca                                                                                    >  2:381511/1‑149 (MQ=255)
                                                             tCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGccacca                                                                                 <  8:35914/149‑1 (MQ=255)
                                                             tCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGccacca                                                                                 <  8:29975/149‑1 (MQ=255)
                                                                ttCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCAccg                                                                              >  5:177168/1‑149 (MQ=255)
                                                                   aGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCAccgccg                                                                           >  6:67054/1‑149 (MQ=255)
                                                                    gggCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACcgccgc                                                                          <  1:381511/149‑1 (MQ=255)
                                                                    gggCGTGGAGGATGTCAATAGCGGTTGCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACcgccgc                                                                          <  4:233987/149‑1 (MQ=255)
                                                                         tGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGtctc                                                                     >  4:198621/1‑149 (MQ=255)
                                                                           gaggaTGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTcgc                                                                   <  7:68698/149‑1 (MQ=255)
                                                                              gaTGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTcgctca                                                                >  8:235665/1‑149 (MQ=255)
                                                                                 gTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTcgctcacgc                                                             >  6:99548/1‑149 (MQ=255)
                                                                                    aaTATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGc                                                                                                                   <  4:340755/92‑1 (MQ=255)
                                                                                    aaTATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGc                                                                                                                   >  4:367621/1‑92 (MQ=255)
                                                                                    aaTATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGc                                                                                                                   <  3:367621/92‑1 (MQ=255)
                                                                                    aaTATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGc                                                                                                                   >  3:340755/1‑92 (MQ=255)
                                                                                    aaTATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTcgctcacgctca                                                          <  6:177168/149‑1 (MQ=255)
                                                                                    aaTATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTcgctcacgctca                                                          <  4:242739/149‑1 (MQ=255)
                                                                                         cGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCAAAAACGCAGAACCAATCCCGACCAGCGTCAAACCACTTGCCGTTACCAGAAATGCCACAACCGCCGCGTCTCGCTCACGCTCATTATg                                                     >  2:45909/1‑149 (MQ=255)
                                                                                          ggTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGc                                                    <  1:40535/149‑1 (MQ=255)
                                                                                           gTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCa                                                   <  7:320378/149‑1 (MQ=255)
                                                                                                      cGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGGTCACGCTCATTATGCAGCGCCTGATAc                                        <  1:175781/149‑1 (MQ=255)
                                                                                                       gATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCAGCGCCTGATAc                                        >  7:5528/1‑148 (MQ=255)
                                                                                                        aTTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCAGCGCCTGATACa                                       >  1:232715/1‑148 (MQ=255)
                                                                                                           aaaTTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCAGCGCCTGATACAAAc                                    <  5:99548/148‑1 (MQ=255)
                                                                                                            aaTTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCAGCGCCTGATACAAACTg                                  >  5:30553/1‑149 (MQ=255)
                                                                                                             aTTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCAGCGCCTGATACAAACTgc                                 <  5:67054/149‑1 (MQ=255)
                                                                                                             aTTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCAGCGCCTGATACAAACTgc                                 >  6:128365/1‑149 (MQ=255)
                                                                                                                     ccacGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCAGCGCCTGATACAAACTGCCGCCGATg                         >  1:158556/1‑149 (MQ=255)
                                                                                                                     ccacGAAACAAAAGCCACCGGCAAGTAATCCCCAAAACGCAGAACCAATCCCGCCCAGCGTCAAACCAATTGCCATTACCAGTAATGCCACCACCGCCGCGTCACGCTAACTCTCATGATCCAGCTCCAGAAACAAACTGCCGCAGATg                         >  7:348461/1‑149 (MQ=255)
                                                                                                                      cacGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCAGCGCCTGATACAAACTGCCGCCGATgg                        >  1:322460/1‑149 (MQ=255)
                                                                                                                       acGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCAGCGCCTGATACAAACTGCCGCCGATgg                        >  1:385078/1‑148 (MQ=255)
                                                                                                                       acGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCAGCGCCTGATACAAACTGCCGCCGATGGt                       <  2:165931/149‑1 (MQ=255)
                                                                                                                        cGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCAGCGCCTGATACAAACTGCCGCCGATGGTa                      <  4:20738/149‑1 (MQ=255)
                                                                                                                            aCAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCAGCGCCTGATACAAACTGCCGCCGATGGTACtt                   <  7:340298/148‑1 (MQ=255)
                                                                                                                             cAAAAGCCTACGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCAGCGCCTGATACAAACTGCCGCCTATGGTACTTaa                 <  3:277383/149‑1 (MQ=255)
                                                                                                                              aaaCGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCAGCGCCTGATACAAACTGCCGCCGATGGTACTTaa                 >  2:224277/1‑148 (MQ=255)
                                                                                                                              aaaCGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCAGCGCCTGATACAAACTGCCGCCGATGGTACTTAAc                <  7:248354/149‑1 (MQ=255)
                                                                                                                              aaaCGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACATGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCAGCGCCTGATACAACATGACGCCGAAGGTACtt                   >  6:323589/1‑146 (MQ=255)
                                                                                                                              aaaCGCCGCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCAGCGCCTTATACAAACTGCCGCCGATGGTACtt                   <  5:323589/146‑1 (MQ=255)
                                                                                                                               aaCGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCAGCGCCTGATACAAACTGCCGCCGATGGTACTTAACa               <  2:232715/149‑1 (MQ=255)
                                                                                                                               aaCGACTCCGGCAATTAATCCCCTAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCTCGTCTCGCTCACGCTCATTAGGCAGCGCCTGATACAAACTGCCGCCGATGGTACTTaaaa               <  4:217369/149‑3 (MQ=255)
                                                                                                                                 cGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACCTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCAGCGCCTGATACAAACTGCCGCCGATGGTACTTAACAgc             >  5:196608/1‑149 (MQ=255)
                                                                                                                                     tCCCGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCAGCGCCTGATACAAACTGCCACCGATGGTACTTAACAgcgc           >  2:35760/1‑147 (MQ=255)
                                                                                                                                       cGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCAGCGCCTGATACAAACTGCAGCCGATGGTACTTAACAGCGCCAGa       >  6:142026/1‑149 (MQ=255)
                                                                                                                                        ggCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCAGCGCCTGATACAAACTGCCGCCGATGGTACTTAACAGCGCCAGAc      <  3:198621/149‑1 (MQ=255)
                                                                                                                                            aTTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCAGCGCCTGATACAAACTGCCGCCGATGGTACTTAACAGCGCCAGACCTGc  <  6:30553/149‑1 (MQ=255)
                                                                                                                                            aTTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCAGCGCCTGATACAAACTGCCGCCGATGGTACTTAACAGCGCCAGACCTGc  <  5:142026/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTTGCCGTTACCAGAAATGCCACCACCGCCGCGTCTCGCTCACGCTCATTATGCAGCGCCTGATACAAACTGCCGCCGATGGTACTTAACAGCGCCAGACCTGC  >  NC_000913/1504823‑1505111

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: