Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,069,829 G→A 6.8% V331M (GTG→ATG) ‡ aroP‑I → aromatic amino acid transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,069,8290GA 6.8% 215.0 / 4.0 73V331M (GTG→ATG) ‡aroP‑Iaromatic amino acid transport protein
Reads supporting (aligned to +/- strand):  ref base G (38/30);  new base A (2/3);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.13e-01

GGCGACCATCCTCAATATCGTGGTGCTGACGGCTGCGCTGTCGGTGTACAACAGCTGTGTGTACGCCAATAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGTTGGCAGTGTCCGCGTTGGTGATCAACTGGGCGAGCATCAGCATTACCCA  >  NC_002947/1069687‑1069977
                                                                                                                                               |                                                                                                                                                     
ggCGACCATCCTCAATATCGTGGTGCTGACGGCTGCGCTGTCGGTGTACAACAGCTGTGTGTACGCCAATAGCCGTATGCTGTATGGCCTGGCAAGCCAGGGCGATGCGCCGCTTCAGTTGCTGAAGGTGAGACGTAGCGGCGTGCCGc                                                                                                                                                 >  6:288735/1‑149 (MQ=255)
 gCGACCATCCTCAATATCGTGGTGCTGACGGCTGCGCTGTCGGTGTACAACAGCTGTGTGTACGCCAATAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGc                                                                                                                                                 <  2:73318/148‑1 (MQ=255)
   gACCATCCTCAATATCGTGGTGCTGACGGCTGCGCTGTCGGTGTACAACAGCTGTGTGTACGCCAATAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGa                                                                                                                                              >  2:332219/1‑149 (MQ=255)
   gACCATCCTCAATATCGTGGTGCTGACGGCTGCGCTGTCGGTGTACAACAGCTGTGTGTACGCCAATAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCACGCCGCtac                                                                                                                                              <  2:279441/149‑3 (MQ=255)
   gACCATCCTCAATATCGTGGTGCTGACGGCTGCGCTGTCGGTGTACAACAGCTGTGTGTACGCCAATAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATACGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCACGCCGCtac                                                                                                                                              <  4:122424/149‑3 (MQ=255)
   gACCATCCTCAATATCGTGGTGCTGACGGCTGCGCTGTCGGTGTACAACAGCTGTGTGTACGCCAATAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATACGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCACGCCGCtac                                                                                                                                              <  4:127341/149‑3 (MQ=255)
       aTCCTCAATATCGTGGTGCTGACGGCTGCGCTGTCGGTGTACAACAGCTGTGTGTTCGCCAATAGCAGTATGCTGTTTTGCATGGCCAGAAAGGGAGATGAGACGCGTCAGGTTGGTAAGGTGAGCAGTAGCGGCGTGCCGGTGACcgc                                                                                                                                          >  6:252132/1‑149 (MQ=255)
             aaTATCGTGGTGCTGACGGCTGCGCTGTCGGTGTACAACAGCTGTGTGTACGCCAATAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGc                                                                                                                                    >  4:100876/1‑149 (MQ=255)
              atatCGTGGTGCTGACGGCTGCGCTGTCGGTGTACAACAGCTGTGTGTACGCCAATAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGcc                                                                                                                                   >  4:68622/1‑149 (MQ=255)
                 tCGTGGTGCTGACGGCTGCGCTGTCGGTGTACAACAGCGGTGTGTACGCCAATAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCgt                                                                                                                                 <  2:79321/148‑1 (MQ=255)
                   gtggtgCTGACGGCTGCGCTGTCGGTGTACAACAGCTGTGTGTACGCCAATAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCgtgt                                                                                                                               <  3:68622/148‑1 (MQ=255)
                      gtgCTGACGGCTGCGCTGTCGGTGTACAACAGCTGTGTGTACGCCAATAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGc                                                                                                                           >  2:243395/1‑149 (MQ=255)
                                 tGCGCTGTCGGTGTACAACAGCTGTGTGTACGCCAATAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCg                                                                                                               >  1:6433/1‑149 (MQ=255)
                                 tGCGCTGTCGGTGTACAACAGCTGTGTGTACGCCAATAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCg                                                                                                               >  1:11958/1‑149 (MQ=255)
                                 tGCGCTGTCGGTGTACAACAGCTGTGTGTACGCCAATAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCg                                                                                                               >  6:25740/1‑149 (MQ=255)
                                   cgcTGTCGGTGTACAACAGCTGTGTGTACGCCAATAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGc                                                                                                             <  4:107646/149‑1 (MQ=255)
                                     cTGTCGGTGTACAACAGCTGTGTGTACGCCAATAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCt                                                                                                           <  6:250097/149‑1 (MQ=255)
                                          ggTGTACAACAGCTGTGTGTACGCCAATAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTCCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTgtgc                                                                                                       <  1:243395/148‑1 (MQ=255)
                                            tgtACAACAGCTGTGTGTACGCCAATAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTgtgcgtg                                                                                                    >  5:145895/1‑149 (MQ=255)
                                              tACAACAGCTGTGTGTACGCCAATAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTgtgcgtgc                                                                                                   <  6:145895/148‑1 (MQ=255)
                                              tACAACAGCTGTGTGTACGCCAATAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTgtgcgtgc                                                                                                   <  1:187375/148‑1 (MQ=255)
                                              tACAACAGCTGTGTGTACGCCAATAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCt                                                                                                  >  2:254791/1‑149 (MQ=255)
                                                 aacaGCTGTGTGTACGCCAATAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGAt                                                                                               >  6:77946/1‑149 (MQ=255)
                                                        gtgtgtACGCCAATAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTAc                                                                                        >  3:282066/1‑149 (MQ=255)
                                                              aCGCCAATAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATg                                                                                  >  4:97784/1‑149 (MQ=255)
                                                                gCCAATAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGcc                                                                                >  4:75930/1‑149 (MQ=255)
                                                                   aaTATCCGTCTGCTGTTTGGCCGTGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCAGCTGACCGAGCTGGCCGTGTCGGCCTTTGCGACCGGc                                                                                                             <  2:207406/117‑1 (MQ=255)
                                                                   aaTAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCGGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGc                                                                                                             >  1:207406/1‑117 (MQ=255)
                                                                     tAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCcgggcg                                                                           <  1:332219/149‑1 (MQ=255)
                                                                             tGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCttt                                                                   >  6:288109/1‑149 (MQ=255)
                                                                             tGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCtt                                                                    <  3:97784/148‑1 (MQ=255)
                                                                             tGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCtt                                                                    <  8:74740/148‑1 (MQ=255)
                                                                              gCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTg                                                                  >  6:158805/1‑149 (MQ=255)
                                                                               cggTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCGGCAGGTGAGCCGTAGCGGCGGGCCGCTGACCGCGCTGGCCGTGGCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTg                                                                  <  2:11958/146‑1 (MQ=255)
                                                                               cTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTg                                                                  <  2:6433/148‑1 (MQ=255)
                                                                                  tttGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGtt                                                               <  5:77946/148‑1 (MQ=255)
                                                                                      gcctggccAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCt                                                           >  5:69594/1‑148 (MQ=255)
                                                                                       cctggccAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGa                                                         >  7:227436/1‑149 (MQ=255)
                                                                                        ctggccAGCCAGGGCGATGCGCCGCGTCAGTTGGTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATGAATTACCTGATGCCGGGCGGGGCCTTGGGTgtggtggt                                                        >  4:112542/1‑147 (MQ=255)
                                                                                         tggccAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATg                                                       >  1:319766/1‑149 (MQ=255)
                                                                                         tggccAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATg                                                       >  1:28203/1‑149 (MQ=255)
                                                                                          ggccAGCCAGGGCGATGCGCCGCGTCAGTTGGTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCGTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATgg                                                      <  5:288735/149‑1 (MQ=255)
                                                                                           gccagccaGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATgg                                                      <  3:100876/148‑1 (MQ=255)
                                                                                           gccagccaGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATg                                                       <  1:254791/147‑1 (MQ=255)
                                                                                              agccagGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGtt                                                  >  1:227568/1‑149 (MQ=255)
                                                                                              agccagGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGtt                                                  <  6:69594/149‑1 (MQ=255)
                                                                                                ccagGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCttt                                                                                                                       <  1:182749/79‑1 (MQ=255)
                                                                                                ccagGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCttt                                                                                                                       >  2:182749/1‑79 (MQ=255)
                                                                                                 cagGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGTTgg                                                >  3:242386/1‑148 (MQ=255)
                                                                                                   gggCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGa                                                         <  7:287174/137‑1 (MQ=255)
                                                                                                   gggCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGa                                                         >  8:287174/1‑137 (MQ=255)
                                                                                                    ggCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGTTGGCAgt                                            <  2:28203/149‑1 (MQ=255)
                                                                                                       gATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGc                                                                                                                           >  4:246792/1‑68 (MQ=255)
                                                                                                       gATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGc                                                                                                                           <  3:246792/68‑1 (MQ=255)
                                                                                                            gccgcGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGTTGGCAGTGTCcgc                                      <  2:298640/147‑1 (MQ=255)
                                                                                                            gccgcGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGTTGGCAGTGTCcgc                                      >  1:298640/1‑147 (MQ=255)
                                                                                                            gccgcGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGTTGGCAGTGTCCGCGt                                    >  3:14032/1‑149 (MQ=255)
                                                                                                               gcgTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGTTGGCAGTGTCCGCGTTgg                                 >  5:179537/1‑149 (MQ=255)
                                                                                                                cgTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGTTGGCAGTGTCCGCGTTGGt                                <  5:158805/149‑1 (MQ=255)
                                                                                                                  tCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGTTGGCAGTGTCCGCGTTGGTGa                              >  8:173628/1‑149 (MQ=255)
                                                                                                                    aGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGTTGGCAGTGTCCGCGTTGGTGAt                             >  3:105531/1‑148 (MQ=255)
                                                                                                                    aGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGTTGGCAGTGTCCGCGTTGGTGATc                            <  7:173628/149‑1 (MQ=255)
                                                                                                                         gCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGTTGGCAGTGTCCGCGTTGGTGATCAAc                         <  4:282066/148‑1 (MQ=255)
                                                                                                                         gCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGTTGGCAGTGTCCGCGTTGGTGATCAACt                        >  7:1015/1‑149 (MQ=255)
                                                                                                                           tGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCAGCCCGGCCTGTGCGTGCTGATCAATTACATGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGTTGGCAGTGTCCGCGTTGGTGATCAACTg                       >  1:266974/1‑149 (MQ=255)
                                                                                                                             aaGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGTTGGCAGTGTCCGCGTTGGTGATCAACTggg                     <  2:319766/148‑1 (MQ=255)
                                                                                                                                  gAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGTTGGCAGTGTCCGCGTTGGTGATCAACTGGGCg                   <  1:53354/145‑1 (MQ=255)
                                                                                                                                  gAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGTTGGCAGTGTCCGCGTTGGTGATCAACTGGGCg                   >  2:53354/1‑145 (MQ=255)
                                                                                                                                    gCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGTTGGCAGTGTCCGCGTTGGTGATCAACTGGGCGAGCAt              >  5:121261/1‑148 (MQ=255)
                                                                                                                                      gtcaGCGGCACGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGTTGGCAGTGTCCGCGTTGGTGATCAACTGGGCGAGCATCAg           >  5:164459/4‑149 (MQ=255)
                                                                                                                                      gtcaGCGGCACGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGTTGGCAGTGTCCGCGTTGGTGATCAACTGGGCGAGCATCAg           >  2:188449/4‑149 (MQ=255)
                                                                                                                                      cGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGTTGGCAGTGTCCGCGTTGGTGATCAACTGGGCGAGCATCAg           >  3:142831/1‑149 (MQ=255)
                                                                                                                                      cGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGTTGGCAGTGTCCGCGTTGGTGATCAACTGGGCGAGCATCAg           <  4:142831/149‑1 (MQ=255)
                                                                                                                                           cggagTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATTCCGGGAGGTGCATTTGGTTTGCTGATGGCGTTGGCAGTGTCCGCGTTGGTGATCAACTGGGCGAGCATCAGCATTa      <  5:189003/149‑1 (MQ=255)
                                                                                                                                              cgTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGTTGGCAGTGTCCGCGTTGGTGATCAACTGGGCGAGCATCAGCATTAccc   <  5:288109/149‑1 (MQ=255)
                                                                                                                                               gTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGTTGGCAGTGTCCGCGTTGGTGATCAACTGGGCGAGCATCAGCATTACCCa  <  6:47010/149‑1 (MQ=255)
                                                                                                                                               gTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGTTGGCAGTGTCCGCGTTGGTGATCAACTGGGCGAGCATCAGCATTACCCa  >  1:6965/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
GGCGACCATCCTCAATATCGTGGTGCTGACGGCTGCGCTGTCGGTGTACAACAGCTGTGTGTACGCCAATAGCCGTATGCTGTTTGGCCTGGCCAGCCAGGGCGATGCGCCGCGTCAGTTGCTGAAGGTGAGCCGTAGCGGCGTGCCGCTGACCGCGCTGGCCGTGTCGGCCTTTGCGACCGGCCTGTGCGTGCTGATCAATTACCTGATGCCGGGCGGGGCCTTTGGTTTGCTGATGGCGTTGGCAGTGTCCGCGTTGGTGATCAACTGGGCGAGCATCAGCATTACCCA  >  NC_002947/1069687‑1069977

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: