Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,512,601 G→C 5.1% A289P (GCG→CCG)  cycA → D‑alanine/beta‑alanine/D‑serine/glycine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,512,6010GC 5.1% 1266.1 / 39.3 372A289P (GCG→CCG) cycAD‑alanine/beta‑alanine/D‑serine/glycine permease
Reads supporting (aligned to +/- strand):  ref base G (139/213);  new base C (9/10);  total (148/224)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 397 total aligned reads displayed.

CCCTGCCGCGGGCGATCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTGTCGAGCCGTTCGGT  >  NC_002947/5512452‑5512749
                                                                                                                                                     |                                                                                                                                                    
 ccTGCCGCGGGCGATCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCcggc                                                                                                                                                     >  2:8562027/1‑150 (MQ=255)
   tGCCGCGGGCGATCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCcggcgg                                                                                                                                                   <  2:7202562/150‑1 (MQ=255)
     ccGCGGGCGATCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCcggcggcg                                                                                                                                                 <  1:8092057/150‑1 (MQ=255)
         gggCGATCAACTCGATCCCGATCAGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCaa                                                                                                                                             <  2:5259490/150‑1 (MQ=255)
             gATCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGcat                                                                                                                                         <  1:4799454/150‑1 (MQ=255)
             gATCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGCCGCCGGCAAGcct                                                                                                                                         <  2:2280793/150‑3 (MQ=255)
              aTCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGcatc                                                                                                                                        <  2:5188132/150‑1 (MQ=255)
              aTCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGcatc                                                                                                                                        <  2:2693643/150‑1 (MQ=255)
               tCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCATGCTTGCCGCCGCCGGCAAGcctgc                                                                                                                                       <  2:654259/150‑5 (MQ=255)
               tCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGcatca                                                                                                                                       <  2:1418627/150‑1 (MQ=255)
               tCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGcatca                                                                                                                                       >  1:8315054/1‑150 (MQ=255)
               tCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGcatca                                                                                                                                       >  1:2086717/1‑150 (MQ=255)
                 aaCTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCatcatc                                                                                                                                     >  1:509927/1‑150 (MQ=255)
                   cTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCaa                                                                                                                                   >  1:5416065/1‑150 (MQ=255)
                    tCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAAc                                                                                                                                  >  1:7284239/1‑150 (MQ=255)
                    tCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAAc                                                                                                                                  >  1:3238372/1‑150 (MQ=255)
                    tCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAAc                                                                                                                                  >  1:7089044/1‑150 (MQ=255)
                     cgatcccgatccGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACt                                                                                                                                 >  2:3071946/1‑150 (MQ=255)
                      gatcccgatccGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACtt                                                                                                                                >  2:7482480/1‑150 (MQ=255)
                          ccgatccGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGTGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCgtg                                                                                                                            >  2:8281654/1‑150 (MQ=255)
                          ccgatccGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGTGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCgtg                                                                                                                            >  2:2525123/1‑150 (MQ=255)
                          ccgatccGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCgtg                                                                                                                            <  1:8769060/150‑1 (MQ=255)
                          ccgatccGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCgtg                                                                                                                            <  2:1341037/150‑1 (MQ=255)
                          ccgatccGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCgtg                                                                                                                            <  2:686877/150‑1 (MQ=255)
                            gatccGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCgtggt                                                                                                                          <  2:7569134/150‑1 (MQ=255)
                            gatccGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCgtggtg                                                                                                                         <  1:2337362/150‑1 (MQ=255)
                                 gTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGa                                                                                                                     <  2:2586164/150‑1 (MQ=255)
                                  tATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGAc                                                                                                                    <  2:4818391/150‑1 (MQ=255)
                                  tATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGAc                                                                                                                    <  2:3624535/150‑1 (MQ=255)
                                     catcGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTc                                                                                                                 <  2:6328477/150‑1 (MQ=255)
                                     catcGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTACTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTc                                                                                                                 <  2:5495713/150‑1 (MQ=255)
                                      atcGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTcg                                                                                                                >  1:4307326/1‑150 (MQ=255)
                                           gtTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCggcggc                                                                                                           <  2:1285407/150‑1 (MQ=255)
                                               tATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTtcg                                                                                                       <  1:5169619/150‑1 (MQ=255)
                                               tATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTtcg                                                                                                       >  1:4649413/1‑150 (MQ=255)
                                                aTGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTtcgt                                                                                                      <  1:6043566/150‑1 (MQ=255)
                                                 tgtgCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTtcgtc                                                                                                     >  2:5595425/1‑150 (MQ=255)
                                                 tgtgCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTtcgtc                                                                                                     <  1:4427378/150‑1 (MQ=255)
                                                  gtgCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTtcgtcg                                                                                                    <  1:1980899/150‑1 (MQ=255)
                                                  gtgCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTtcgtcg                                                                                                    <  1:7464992/150‑1 (MQ=255)
                                                      tGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCa                                                                                                >  2:284243/1‑150 (MQ=255)
                                                      tGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCa                                                                                                >  2:734399/1‑150 (MQ=255)
                                                           cTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAgc                                                                                           <  1:3070610/150‑1 (MQ=255)
                                                              aTCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAgcggc                                                                                        <  2:625376/150‑1 (MQ=255)
                                                              aTCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAgcggc                                                                                        <  1:5633946/150‑1 (MQ=255)
                                                                cgcgATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCgt                                                                                      <  1:7886879/150‑1 (MQ=255)
                                                                 gcgATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCgtg                                                                                     <  2:7405505/150‑1 (MQ=255)
                                                                  cgATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCgtgt                                                                                    <  2:2377303/150‑1 (MQ=255)
                                                                      cATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTtctc                                                                                <  1:2420562/150‑1 (MQ=255)
                                                                      cATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGGTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTtctc                                                                                <  2:4680294/150‑1 (MQ=255)
                                                                       agggCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTcc                                                                               <  1:5935082/148‑1 (MQ=255)
                                                                        tGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTcca                                                                              >  1:3456767/1‑150 (MQ=255)
                                                                        tGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTcca                                                                              <  1:8562027/150‑1 (MQ=255)
                                                                                aCCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGc                                                                      >  2:798801/1‑150 (MQ=255)
                                                                                aCCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGc                                                                      >  2:1803702/1‑150 (MQ=255)
                                                                                  cccGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCAt                                                                    >  1:6471496/1‑150 (MQ=255)
                                                                                  cccGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCAt                                                                    >  2:7760730/1‑150 (MQ=255)
                                                                                     gTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGct                                                                 <  1:4947818/150‑1 (MQ=255)
                                                                                            gATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGc                                                          >  2:5895251/1‑150 (MQ=255)
                                                                                               gtggtgCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTg                                                       <  2:6744120/150‑1 (MQ=255)
                                                                                                    gCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCa                                                  <  1:6935171/150‑1 (MQ=255)
                                                                                                      ccGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAgg                                                >  2:1107301/1‑150 (MQ=255)
                                                                                                      ccGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAgg                                                >  1:672962/1‑150 (MQ=255)
                                                                                                       cGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGa                                               >  2:1663019/1‑150 (MQ=255)
                                                                                                        ggAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGaa                                              >  2:4550974/1‑150 (MQ=255)
                                                                                                        ggAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGaa                                              <  1:4016612/150‑1 (MQ=255)
                                                                                                          aaaaaGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAgg                                            >  1:2686903/1‑150 (MQ=255)
                                                                                                            aaaGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTg                                          <  2:3851055/150‑1 (MQ=255)
                                                                                                              aGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGAc                                        <  1:2938797/150‑1 (MQ=255)
                                                                                                              aGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGAc                                        <  2:7059462/150‑1 (MQ=255)
                                                                                                               gCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGAcg                                       >  1:3687199/1‑150 (MQ=255)
                                                                                                                  cgttcgttGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACgcgc                                    >  1:1097893/1‑150 (MQ=255)
                                                                                                                  cgttcgttGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACgcgc                                    >  2:3319019/1‑150 (MQ=255)
                                                                                                                   gttcgttGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGcc                                   >  1:4520252/1‑150 (MQ=255)
                                                                                                                   gttcgttGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGcc                                   <  2:5770978/150‑1 (MQ=255)
                                                                                                                      cgttGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCaa                                <  1:139332/150‑1 (MQ=255)
                                                                                                                       gttGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAg                               <  2:3382715/150‑1 (MQ=255)
                                                                                                                        ttGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAgg                              >  2:2165404/1‑150 (MQ=255)
                                                                                                                          gAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAg                                                 >  1:769913/1‑129 (MQ=255)
                                                                                                                             ttgttCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTc                         >  2:8714638/1‑150 (MQ=255)
                                                                                                                              tgttCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCg                        <  1:6343230/150‑1 (MQ=255)
                                                                                                                              tgttCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCg                        <  2:4135343/150‑1 (MQ=255)
                                                                                                                               gttCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCga                       <  2:2086717/150‑1 (MQ=255)
                                                                                                                               gttCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCga                       <  1:6630114/150‑1 (MQ=255)
                                                                                                                               gttCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCga                       <  1:1032145/150‑1 (MQ=255)
                                                                                                                                ttCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCgag                      <  2:5953799/150‑1 (MQ=255)
                                                                                                                                 tCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCgaga                     >  2:3534162/1‑150 (MQ=255)
                                                                                                                                  cGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGaa                    <  2:4397814/150‑1 (MQ=255)
                                                                                                                                  cGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGaa                    <  1:1762249/150‑1 (MQ=255)
                                                                                                                                      cTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTAGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTg                <  1:6113187/150‑1 (MQ=255)
                                                                                                                                        ggCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTGTc              <  1:656616/150‑1 (MQ=255)
                                                                                                                                          gatgCTTGCCGCCGCCGGCAAGCCTCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTGTCGa            >  2:2898672/4‑150 (MQ=255)
                                                                                                                                          cAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTGTCGa            >  2:7284040/1‑150 (MQ=255)
                                                                                                                                           aGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTgccctc           >  1:2595127/1‑145 (MQ=255)
                                                                                                                                            tgCTTGCCGCCGCCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTGTCGAGc          >  1:2796871/2‑150 (MQ=255)
                                                                                                                                             gCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTGTCGAGcc         <  1:4089837/150‑1 (MQ=255)
                                                                                                                                             gCTTGCCGCCGCCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTGTCGAGcc         >  1:2308163/1‑150 (MQ=255)
                                                                                                                                               ttGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTc                         <  1:5772174/132‑1 (MQ=255)
                                                                                                                                               ttGCCGCCGCCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTGTCGAGCCGt       >  2:1037826/1‑150 (MQ=255)
                                                                                                                                                 gCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTGTCGAGCCGTTc     <  1:122176/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCCTGCCGCGGGCGATCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTGTCGAGCCGTTCGGT  >  NC_002947/5512452‑5512749

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