Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,085,863 A→G 6.7% I188V (ATC→GTC)  amaD → D‑lysine oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,085,8630AG 6.7% 220.2 / 4.8 74I188V (ATC→GTC) amaDD‑lysine oxidase
Reads supporting (aligned to +/- strand):  ref base A (34/35);  new base G (2/3);  total (36/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 = 5.68e-01

GATCGAAGTGTTCCGCGACCCCGCGCTGTTCGAGCAGGCAAAAACCGATGCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCG  >  NC_002947/4085715‑4086011
                                                                                                                                                    |                                                                                                                                                    
gATCGAAGTGTTCCGCGACCCCGCGCTGTTCGAGCAGGCAAAAACCGATGCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCa                                                                                                                                                      <  4:362397/149‑1 (MQ=255)
 aTCGAAGTGTTCCGCGACCCCGCGCTGTTCGAGCAGGCAAAAACCGATGCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCAt                                                                                                                                                     <  7:158444/149‑1 (MQ=255)
  tCGAAGTGTTCCGCGACCCCGCGCTGTTCGAGCAGGCAAAAACCGATGCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATc                                                                                                                                                    >  3:302409/1‑149 (MQ=255)
  tCGAAGGGTTCCGCGACCCCGCGCTGTTCGAGCAGGCAAAAACCGATGCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATc                                                                                                                                                    <  2:292244/149‑1 (MQ=255)
     aaGTGTTCCGCGACCCCGCGCTGTTCGAGCAGGCAAAAACCGATGCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGATGCCACCGACCACGGTGGCGTCCag                                                                                                                                                 <  5:219699/149‑2 (MQ=255)
      aGTGTTCCGCGACCCCGCGCTGTTCGAGCAGGCAAAAACCGATGCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGATGCCACCGACCACGGTGGCGTCCag                                                                                                                                                 <  6:113703/148‑2 (MQ=255)
       gtgtTCCGCGACCCCGCGCTGTTCGAGCAGGCAAAAACCGATGCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTg                                                                                                                                               <  6:118357/149‑1 (MQ=255)
            ccGCGACCCCGCGCTGTTCGAGCAGGCAAAAACCGATGCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTc                                                                                                                                           >  4:120556/1‑148 (MQ=255)
             cgcgACCCCGCGCTGTTCGAGCAGGCAAAAACCGATGCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCg                                                                                                                                          <  4:302409/148‑1 (MQ=255)
             cgcgACCCCGCGCTGTTCGAGCAGGCAAAAACCGATGCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGa                                                                                                                                         >  7:145672/1‑149 (MQ=255)
             cgcgACCCCGCGCTGTTCGAGCAGGCAAAAACCGATGCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGa                                                                                                                                         <  1:179002/149‑1 (MQ=255)
              gcgACCCCGCGCTGTTCGAGCAGGCAAAAACCGATGCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGAc                                                                                                                                        >  6:59621/1‑149 (MQ=255)
              gcgACCCCGCGCTGTTCGAGCAGGCAAAAACCGATGCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGAc                                                                                                                                        <  8:145672/149‑1 (MQ=255)
               cgACCCCGCGCTGTTCGAGCAGGCAAAAACCGATGCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGAcc                                                                                                                                       >  2:366498/1‑149 (MQ=255)
                         cTGTTCGAGCAGGCAAAAACCGATGCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGt                                                                                                                             <  7:278971/149‑1 (MQ=255)
                         cTGTTCGAGCAGGCAAAAACCGATGCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGt                                                                                                                             <  8:33550/149‑1 (MQ=255)
                               gAGCAGGCAAAAACCGATGCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGaacaa                                                                                                                       >  4:344811/1‑149 (MQ=255)
                                    ggCAAAAACCGATGCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTg                                                                                                                                               >  6:268610/1‑120 (MQ=255)
                                    ggCAAAAACCGATGCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTg                                                                                                                                               <  5:268610/120‑1 (MQ=255)
                                     gCAAAAACCGATGCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCgg                                                                                                                 >  7:179409/1‑149 (MQ=255)
                                       aaaaaCCGATGCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCATTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCTGTGGCATCCACTGGCTCGACCCGAAAACCCTGAACAACACCGgcg                                                                                                               >  8:163889/1‑149 (MQ=255)
                                          aaCCGATGCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGgcgcg                                                                                                             >  2:206786/1‑148 (MQ=255)
                                           aCCGATGCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCt                                                                                                           <  8:179409/149‑1 (MQ=255)
                                             cGATTCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGACCATCAGCTGGCCGCCACCGTGGTCGGTGGCATCCACTGGCTCGCCCCGAAAACCGTGAACAACCCAGCCGAGCTCa                                                                                                         >  4:31382/1‑149 (MQ=255)
                                             cGATGCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCa                                                                                                         >  4:357396/1‑149 (MQ=255)
                                                   cAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAAccc                                                                                                                    >  1:268455/1‑132 (MQ=255)
                                                   cAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAAccc                                                                                                                    <  2:268455/132‑1 (MQ=255)
                                                   cAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTg                                                                                                   <  3:120556/149‑1 (MQ=255)
                                                   cAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTg                                                                                                   <  6:41738/149‑1 (MQ=255)
                                                    aaGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTgg                                                                                                  <  3:344811/149‑1 (MQ=255)
                                                             aGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATgcg                                                                                           <  1:206786/147‑1 (MQ=255)
                                                                cGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCgtt                                                                                      >  8:227251/1‑149 (MQ=255)
                                                                  cTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGc                                                                                                                                             >  5:336405/1‑92 (MQ=255)
                                                                  cTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGc                                                                                                                                             <  6:336405/92‑1 (MQ=255)
                                                                     cGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTc                                                                                  <  7:163889/148‑1 (MQ=255)
                                                                       gCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTg                                                                               <  4:131281/149‑1 (MQ=255)
                                                                        cTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTg                                                                               >  5:176863/1‑148 (MQ=255)
                                                                        cTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTg                                                                               >  2:160740/1‑148 (MQ=255)
                                                                         ttGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCgttgtt                                                                                   >  2:1098/1‑143 (MQ=255)
                                                                         ttGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCgttgtt                                                                                   <  1:1098/143‑1 (MQ=255)
                                                                           gcgcTTCGAAATCCTTGAATGCGGTCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACAGTGGTCGGTGGCATCCACTGGCTCGACCCGATAACAGTGAACAACCCCGGCGCGCTCACCCGTGTCTATGCGGCGTTGTTCCTGc                                                                              <  3:31382/146‑1 (MQ=255)
                                                                                  gAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGcggcgg                                                                    >  1:195739/1‑149 (MQ=255)
                                                                                    aaTCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCggg                                                                   >  1:182974/1‑148 (MQ=255)
                                                                                    aaTCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCggg                                                                   <  2:182974/148‑1 (MQ=255)
                                                                                         ttGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTc                                                             >  3:280959/1‑149 (MQ=255)
                                                                                           gAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGt                                                           >  4:361913/1‑149 (MQ=255)
                                                                                                 gggCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATgg                                                     <  6:176863/149‑1 (MQ=255)
                                                                                                     aGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGAt                                                 <  2:195739/149‑1 (MQ=255)
                                                                                                       ttGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATgc                                               <  3:357396/149‑1 (MQ=255)
                                                                                                       ttGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATgc                                               >  6:137206/1‑149 (MQ=255)
                                                                                                           aGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTg                                                                                                   <  7:251589/93‑1 (MQ=255)
                                                                                                           aGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTg                                                                                                   >  8:251589/1‑93 (MQ=255)
                                                                                                           aGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATgcgcg                                            >  6:197012/1‑148 (MQ=255)
                                                                                                            ggCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCa                                          <  5:59621/149‑1 (MQ=255)
                                                                                                             gCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTgg                                                                                                  >  7:355784/1‑92 (MQ=255)
                                                                                                             gCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTgg                                                                                                  <  8:355784/92‑1 (MQ=255)
                                                                                                             gCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCa                                                        <  3:302656/134‑1 (MQ=255)
                                                                                                             gCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCa                                                        >  4:302656/1‑134 (MQ=255)
                                                                                                             gCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAg                                         >  2:58601/1‑149 (MQ=255)
                                                                                                               aCGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGcc                                       >  6:298540/1‑149 (MQ=255)
                                                                                                                 gTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATgcgg                                                                                          <  3:227765/96‑1 (MQ=255)
                                                                                                                 gTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATgcgg                                                                                          >  4:227765/1‑96 (MQ=255)
                                                                                                                  tGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTgc                                    <  4:280959/149‑1 (MQ=255)
                                                                                                                          cagtGGATGCCACCGACCACGGTGGCGTCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGAGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAggc                             >  5:198972/4‑148 (MQ=255)
                                                                                                                           gCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTg                                                                                                   <  2:5327/77‑1 (MQ=255)
                                                                                                                           gCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTg                                                                                                   >  1:5327/1‑77 (MQ=255)
                                                                                                                           agtGGATGCCACCGACCACGGTGGCGTCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAggcca                           >  3:270119/3‑149 (MQ=255)
                                                                                                                               gACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGTAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAggccaatgg                       <  7:227251/149‑1 (MQ=255)
                                                                                                                               gACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAggccaatg                        <  5:197012/148‑1 (MQ=255)
                                                                                                                                         tGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCgg             >  1:202606/1‑149 (MQ=255)
                                                                                                                                          ggTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCgttgt                                                                                    <  2:966/77‑1 (MQ=255)
                                                                                                                                          ggTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCgttgt                                                                                    >  1:966/1‑77 (MQ=255)
                                                                                                                                            tCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTa          <  1:58601/149‑1 (MQ=255)
                                                                                                                                             cGGTGGCGTCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTag         <  2:166890/149‑1 (MQ=255)
                                                                                                                                                  gCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGc    >  7:15379/1‑149 (MQ=255)
                                                                                                                                                    aTCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCgc                                                  <  6:254342/101‑2 (MQ=255)
                                                                                                                                                    aTCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCgc                                                  >  5:254342/1‑100 (MQ=255)
                                                                                                                                                    aTCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCTGGTAGAGAGCCg  >  2:79662/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GATCGAAGTGTTCCGCGACCCCGCGCTGTTCGAGCAGGCAAAAACCGATGCCAAGGGCCTGAGCCGCTACGGCTTGCGCTTCGAAATCCTTGAATGCGGGCAGTTGCAGGCACGTGAACATCAGCTGGACGCCACCGTGGTCGGTGGCATCCACTGGCTCGACCCGAAAACCGTGAACAACCCCGGCGCGCTCACCCGTGGCTATGCGGCGTTGTTCCTGCAACGCGGCGGGCAGTTCGTGCATGGCGATGCGCGCAGCCTGCGCCAGGCCAATGGCCAATGGCGGGTAGAGAGCCG  >  NC_002947/4085715‑4086011

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: ‑