Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,569,312 G→T 7.5% Q623H (CAG→CAT) ‡ PP_3150 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,569,3120GT 7.5% 243.1 / 7.9 80Q623H (CAG→CAT) ‡PP_3150hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (37/37);  new base T (3/3);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

ACACGCCGGAGCACTTTGCCAAAGCCAAGAGCCGTTACCGCTGGGACAAGGCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAG  >  NC_002947/3569170‑3569460
                                                                                                                                              |                                                                                                                                                    
acacGCCGGAGCACTTTGCCAAAGCCAAGAGCCGTTACCGCTGGGACAAGGCGCCGAAGCTGTTCCTGGGGCGAGTACGGGTTGTGAAGAACATGGAGTTTGCCGGGGTGATGGTGTCGCGGGTGCTGATGGTGGCGCCACAGGTGGCg                                                                                                                                                >  4:326800/1‑149 (MQ=255)
 cacGCCGGAGCACTTTGCCAAAGCCAAGAGCCGTTACCGCTGGGACAAGGCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGc                                                                                                                                                 <  7:275286/147‑1 (MQ=255)
 cacGCCGGAGCACTTTGCCAAAGCCAAGAGCCGTTACCGCTGGGACAAGGCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCgg                                                                                                                                               >  6:180492/1‑149 (MQ=255)
 cacGCCGGAGCACTTTGCCAAAGCCAAGAGCCGTTACCGCTGGGACAAGGCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCgc                                                                                                                                               >  4:352289/1‑148 (MQ=255)
 cacGCCGGAGCACTTTGCCAAAGCCAAGAGCCGTTACCGCTGGGACAAGGCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCg                                                                                                                                                <  3:344189/148‑1 (MQ=255)
  acGCCGGAGCACTTTGCCAAAGCCAAGAGCCGTTACTGCTGGGACAAGGCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCggg                                                                                                                                              >  3:176806/1‑149 (MQ=255)
    gCCGGAGCACTTTGCCAAAGCCAAGAGCCGTTACCGCTGGGACAAGGCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACCTGTGGCgccac                                                                                                                                            <  3:255848/149‑5 (MQ=255)
    gCCGGAGCACTTTGCCAAAGCCAAGAGCCGTTACCGCTGGGACAAGGCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGt                                                                                                                                             >  3:378815/1‑148 (MQ=255)
                   cAAAGCCAAGAGCCGTTACCGCTGGGACAAGGCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACggg                                                                                                                             >  7:122094/1‑149 (MQ=255)
                    aaaGCCAAGAGCCGTTACCGCTGGGACAAGGCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGGGCTGATGGTGGCGCCACAGGTGGCGGGTTGCCGCAATGAACGGgc                                                                                                                            >  8:249117/1‑149 (MQ=255)
                    aaaGCCAAGAGCCGTTACCGCTGGGACAAGGCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCgg                                                                                                                                               <  2:217913/130‑1 (MQ=255)
                    aaaGCCAAGAGCCGTTACCGCTGGGACAAGGCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCgg                                                                                                                                               >  1:217913/1‑130 (MQ=255)
                    aaaGCCAAGAGCCGTTACCGCTGGGACAAGGCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGgc                                                                                                                            >  1:55206/1‑149 (MQ=255)
                    aaaGCCAAGAGCCGTTACCGCTGGGACAAGGCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGgc                                                                                                                            >  2:257139/1‑149 (MQ=255)
                      aGCCAACAGCCGTTACCGCTGGGACAAGGCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGgcgc                                                                                                                          <  3:326800/149‑1 (MQ=255)
                           agagCCGTTACCGCTGGGACAAGGCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATg                                                                                                                     <  1:226327/149‑1 (MQ=255)
                           agagCCGTTACCGCTGGGACAAGGCGCCGAAGCTGTTCCTGGGGAGAGTACGGGATGTGAAGAACATTGAGTTTGCCCGGGTGATGCTGTCGCAGGAGCTGAAGGTGGCGCCACAGGTGGCGGGTAACCGGAAAGAACGGGCGCTTATg                                                                                                                     >  8:269746/1‑149 (MQ=255)
                            gagCCGTTACTGCTGGGACAAGGCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGc                                                                                                                    <  4:176806/149‑1 (MQ=255)
                                     ccgcgcggACAAGGCGCCGAAGCTGGTCCTGGGGCGAGTACGGTATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACAGCAATGAACGGGCGCGTATGCTGGCgggg                                                                                                            <  7:269746/142‑1 (MQ=255)
                                        cTGGGACAAGGCGCCGAAGCTGTTCGTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGt                                                                                                        <  8:122094/149‑1 (MQ=255)
                                        cTGGGAACAGGCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAAGATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGt                                                                                                        <  3:118049/149‑1 (MQ=255)
                                            gACAAGGCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCAACAGGTGGCGGGTAACCGCAATGAACCGGCGCGTATGCTGGCGGGGTTGTATAc                                                                                                    >  5:392860/1‑149 (MQ=255)
                                             aCAAGGCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATAcc                                                                                                   <  1:257139/149‑1 (MQ=255)
                                                 ggCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCgg                                                                                                                                               >  2:326120/1‑101 (MQ=255)
                                                 ggCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCgg                                                                                                                                               <  1:326120/101‑1 (MQ=255)
                                                  gcgcCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTg                                                                                              >  6:73348/1‑149 (MQ=255)
                                                     ccGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTc                                                                                           >  6:269766/1‑149 (MQ=255)
                                                     ccGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTc                                                                                           <  5:257266/149‑1 (MQ=255)
                                                        aaGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGGTGTATACCGCCTGTTcgcg                                                                                        <  3:111726/149‑1 (MQ=255)
                                                          gCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCggg                                                                                      >  5:111356/1‑149 (MQ=255)
                                                            tGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGgcg                                                                                    >  7:86935/1‑149 (MQ=255)
                                                             gTTCCTGGGGCGAGTACGGGATGTGAAGAACATGCAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGgcgc                                                                                   >  1:140741/1‑149 (MQ=255)
                                                                ccTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCgg                                                                                 >  1:302615/1‑148 (MQ=255)
                                                                ccTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGc                                                                                <  6:259810/149‑1 (MQ=255)
                                                                 cTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGc                                                                                <  5:180492/148‑1 (MQ=255)
                                                                    gtgCGAGTACGGGATGTGAAGAAGATGGAGGTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTTACCGCAATGAACGGGCGCTTATGCTTGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCAt                                                                              <  5:36422/145‑1 (MQ=255)
                                                                    gggCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGGCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGa                                                                            <  7:249117/149‑1 (MQ=255)
                                                                       cGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAAtt                                                                         <  2:302615/149‑1 (MQ=255)
                                                                         aGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGa                                                                       >  1:39857/1‑149 (MQ=255)
                                                                          gTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGAt                                                                      >  3:297361/1‑149 (MQ=255)
                                                                              gggATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACAGCATGTTCGCGGGCGCGGCATGAATTGATCGTg                                                                  >  7:183498/1‑149 (MQ=255)
                                                                                   gtgAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCgg                                                              <  2:140741/148‑1 (MQ=255)
                                                                                     gaagaaCATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCggggg                                                           >  3:83675/1‑149 (MQ=255)
                                                                                      aagaaCATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCggggg                                                           <  1:235452/148‑1 (MQ=255)
                                                                                       agaaCATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCGGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCgggggg                                                          <  4:378815/148‑1 (MQ=255)
                                                                                        gaaCATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGgc                                                                                                                            >  4:114460/1‑81 (MQ=255)
                                                                                        gaaCATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGgc                                                                                                                            <  3:114460/81‑1 (MQ=255)
                                                                                        gaaCATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGa                                                        <  2:359441/149‑1 (MQ=255)
                                                                                             tGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTg                                                   >  2:133134/1‑149 (MQ=255)
                                                                                             tGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTg                                                   >  1:204657/1‑149 (MQ=255)
                                                                                               gAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGa                                                 <  5:163157/149‑1 (MQ=255)
                                                                                                        cGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCaa                                        <  2:39857/149‑1 (MQ=255)
                                                                                                          tgTGATGGTGTCGCCGCAGCAGATGGTGGCGCCACAGGTGGCGGGTAACCGAAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAgg                                      <  8:273275/148‑1 (MQ=255)
                                                                                                            tGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGAc                                    >  3:43603/1‑149 (MQ=255)
                                                                                                            tGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGAc                                    >  3:311297/1‑149 (MQ=255)
                                                                                                              aTGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGAcc                                   <  5:369275/148‑1 (MQ=255)
                                                                                                                ggTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAg                                 <  5:73348/148‑1 (MQ=255)
                                                                                                                 gtgtCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGt                               >  1:247676/1‑149 (MQ=255)
                                                                                                                  tgtCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGtt                              <  4:161569/149‑1 (MQ=255)
                                                                                                                    tCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCg                            >  3:322487/1‑149 (MQ=255)
                                                                                                                         ggAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGAtt                       <  4:297361/149‑1 (MQ=255)
                                                                                                                          gAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGTGGATGCTGGACTGGTTCAAGGACCAGGTTCGGGATTg                      <  8:86935/149‑1 (MQ=255)
                                                                                                                          gAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTg                      <  3:352289/149‑1 (MQ=255)
                                                                                                                              tGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAag                   >  1:240753/1‑148 (MQ=255)
                                                                                                                              tGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAag                   >  1:372093/1‑148 (MQ=255)
                                                                                                                                 tggtggCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAagcagc               >  2:268294/1‑149 (MQ=255)
                                                                                                                                 gttacccgcCACCTGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAgg                                      <  5:209611/120‑1 (MQ=255)
                                                                                                                                 gttacccgcCACCTGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAgg                                      >  6:209611/7‑126 (MQ=255)
                                                                                                                                  ggtggCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCt              >  3:110590/1‑149 (MQ=255)
                                                                                                                                     ggCTCCACTGTTGGCGGGTAACCTCAATGAACGGTTGCGTATGCTGGCTGGGTTGTATATCGACTGTGCGCGTGTGCGGGATGAATTGATAGTGGCTGGGGGGATGCTGGACTGGGTCAAGGACGAGGTTCTGGATTGAAGCAGCTgg            <  8:183498/148‑1 (MQ=255)
                                                                                                                                     cccgcCACCTGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTggg           >  8:21156/3‑149 (MQ=255)
                                                                                                                                      gcgcCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTgggg          >  6:114207/1‑149 (MQ=255)
                                                                                                                                      ccgcCACCTGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGaa                    >  1:194691/2‑139 (MQ=255)
                                                                                                                                      ccgcCACCTGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGaa                    <  2:194691/138‑1 (MQ=255)
                                                                                                                                       cgcCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGc         <  4:110590/149‑1 (MQ=255)
                                                                                                                                         ccACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCg       >  3:305354/1‑149 (MQ=255)
                                                                                                                                           acaGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAg                                       >  3:155653/1‑115 (MQ=255)
                                                                                                                                           acaGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAg                                       <  4:155653/115‑1 (MQ=255)
                                                                                                                                             aGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCaaa   <  4:311297/149‑1 (MQ=255)
                                                                                                                                             aGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCaaa   <  4:43603/149‑1 (MQ=255)
                                                                                                                                              ggTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGGTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGGGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCGGGGGCCGCAAag  <  2:158345/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
ACACGCCGGAGCACTTTGCCAAAGCCAAGAGCCGTTACCGCTGGGACAAGGCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAAG  >  NC_002947/3569170‑3569460

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