Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,569,311 A→C 7.9% Q623P (CAG→CCG) ‡ PP_3150 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,569,3110AC 7.9% 219.0 / 8.4 76Q623P (CAG→CCG) ‡PP_3150hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (36/33);  new base C (3/3);  total (39/37)
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.93e-01

AAAGGCTACACGCCGGAGCACTTTGCCAAAGCCAAGAGCCGTTACCGCTGGGACAAGGCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAA  >  NC_002947/3569163‑3569459
                                                                                                                                                    |                                                                                                                                                    
aaaGGCTACACGCCGGAGCACTTTGCCAAAGCCAAGAGCCGTTACCGCTGGGACAAGGCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCcaca                                                                                                                                                      <  4:299944/149‑1 (MQ=255)
       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                                                                                                                                              <  1:326120/101‑1 (MQ=255)
                                                        ggCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCgg                                                                                                                                              >  2:326120/1‑101 (MQ=255)
                                                         gcgcCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTg                                                                                             >  6:73348/1‑149 (MQ=255)
                                                            ccGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTc                                                                                          <  5:257266/149‑1 (MQ=255)
                                                            ccGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTc                                                                                          >  6:269766/1‑149 (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                                                  >  1:204657/1‑149 (MQ=255)
                                                                                                    tGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTg                                                  >  2:133134/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                                     >  6:209611/7‑126 (MQ=255)
                                                                                                                                        gttacccgcCACCTGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAgg                                     <  5:209611/120‑1 (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)
                                                                                                                                                    |                                                                                                                                                    
AAAGGCTACACGCCGGAGCACTTTGCCAAAGCCAAGAGCCGTTACCGCTGGGACAAGGCGCCGAAGCTGTTCCTGGGGCGAGTACGGGATGTGAAGAACATGGAGTTTGCCCGGGTGATGGTGTCGCCGGAGCTGATGGTGGCGCCACAGGTGGCGGGTAACCGCAATGAACGGGCGCGTATGCTGGCGGGGTTGTATACCGCCTGTTCGCGGGCGCGGCATGAATTGATCGTGCCGGGGGGGATGCTGGACTGGGTCAAGGACCAGGTTCGGGATTGAAGCAGCTGGGGCCGCAAA  >  NC_002947/3569163‑3569459

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