Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,306,064 C→T 5.0% A79V (GCC→GTC)  PP_3774 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,306,0640CT 5.0% 1152.1 / 28.6 356A79V (GCC→GTC) PP_3774hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (174/164);  new base T (5/13);  total (179/177)
Fisher's exact test for biased strand distribution p-value = 5.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 370 total aligned reads displayed.

GAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACATCACCCTGGGGTTTCGGC  >  NC_002947/4305916‑4306213
                                                                                                                                                    |                                                                                                                                                     
gAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGcc                                                                                                                                                      <  2:3996606/150‑1 (MQ=255)
gAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGcc                                                                                                                                                      >  1:3331515/1‑150 (MQ=255)
  aGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCACAGGGGCTTCGAGCAAAAGCACCCACAGGACCCGCTGCGAATGGACGGAAACGCAAGACTGGGCGCTGGCCAGCGTCCa                                                                                                                                                    <  1:6220681/150‑1 (MQ=255)
  aGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGGCGCTGGCCAGCGTCCa                                                                                                                                                    <  2:4279193/150‑1 (MQ=255)
  aGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGGCGCTGGCCAGCGTCCa                                                                                                                                                    <  1:2755974/150‑1 (MQ=255)
    cGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGGCGCTGGCCAGCGTCCAgt                                                                                                                                                  <  2:3183237/150‑2 (MQ=255)
         gCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCACAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGGCGCTGGCCAGCGTCCAgtcttgc                                                                                                                                             <  1:3530157/150‑7 (MQ=255)
          cccGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCATAGGGGCTTCGAGCAAAAGCACCCCCAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACt                                                                                                                                            <  2:3713575/150‑1 (MQ=255)
          cccGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACt                                                                                                                                            <  1:4634792/150‑1 (MQ=255)
           ccGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTg                                                                                                                                           >  2:2345583/1‑150 (MQ=255)
              cTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGaa                                                                                                                                        >  2:2396537/1‑150 (MQ=255)
                   ccAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCAGGCCAGCGCCCAGGATCACTGGAATCGcc                                                                                                                                   >  2:4895269/1‑150 (MQ=255)
                    cAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCa                                                                                                                                  >  1:4773988/1‑150 (MQ=255)
                    cAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCa                                                                                                                                  >  2:4729383/1‑150 (MQ=255)
                    cAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCa                                                                                                                                  <  2:3631396/150‑1 (MQ=255)
                     aGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAg                                                                                                                                 <  1:3275321/150‑1 (MQ=255)
                     aGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAg                                                                                                                                 >  2:2615364/1‑150 (MQ=255)
                           gTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATg                                                                                                                           <  1:3420364/150‑1 (MQ=255)
                            tGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAATGGGCTTCGAGCAAAAGAACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGa                                                                                                                          >  2:1971978/1‑150 (MQ=255)
                              aCGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACg                                                                                                                        >  1:6489099/1‑150 (MQ=255)
                                ggTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCt                                                                                                                      <  2:1603711/150‑1 (MQ=255)
                                 gTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTg                                                                                                                     <  1:4236634/150‑1 (MQ=255)
                                  tGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGa                                                                                                                    <  1:3311985/150‑1 (MQ=255)
                                    cTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGTGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGATGAcc                                                                                                                  <  2:6678740/150‑1 (MQ=255)
                                     tCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGAccc                                                                                                                 >  2:2697781/1‑149 (MQ=255)
                                     tCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCa                                                                                                                 >  1:2084982/1‑150 (MQ=255)
                                     tCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCa                                                                                                                 >  1:4334691/1‑150 (MQ=255)
                                     tCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCa                                                                                                                 >  1:1006751/1‑150 (MQ=255)
                                      cAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCaa                                                                                                                >  1:1114144/1‑150 (MQ=255)
                                      cAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCaa                                                                                                                <  2:4060383/150‑1 (MQ=255)
                                         cGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGAGCACTGGACTCGCCCGACACTGACGCTGACCACggg                                                                                                             >  1:6745822/1‑149 (MQ=255)
                                            tGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCaa                                                                                                          <  2:3679267/150‑1 (MQ=255)
                                             ggCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAAc                                                                                                         <  2:6377405/150‑1 (MQ=255)
                                                  gATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGgc                                                                                                    <  2:4820807/150‑1 (MQ=255)
                                                  gATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGGCAACTCGgc                                                                                                    >  1:1826718/1‑150 (MQ=255)
                                                   aTCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGgcg                                                                                                   >  2:5523583/1‑150 (MQ=255)
                                                    tCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGgcgc                                                                                                  >  1:5220304/1‑150 (MQ=255)
                                                     ccAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCa                                                                                                 >  1:2212205/1‑150 (MQ=255)
                                                      cAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAg                                                                                                >  2:5064843/1‑150 (MQ=255)
                                                       aGTTCTCAGCCCCAACCAAATGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGTACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGAACATGGTCAACTCTGCGCAGa                                                                                               >  2:5575354/1‑150 (MQ=255)
                                                        gTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGaa                                                                                              <  2:324905/150‑1 (MQ=255)
                                                          tcgaaGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAgg                                                                                            >  1:4150299/5‑150 (MQ=255)
                                                            tcAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAggtg                                                                                          >  2:833347/1‑150 (MQ=255)
                                                            tcAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAggtg                                                                                          >  2:5711682/1‑150 (MQ=255)
                                                            tcAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAggtg                                                                                          <  1:170442/150‑1 (MQ=255)
                                                             cAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAggtgg                                                                                         <  1:5033074/150‑1 (MQ=255)
                                                              aGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAggtggt                                                                                        >  1:760277/1‑150 (MQ=255)
                                                                ccccAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCa                                                                                      <  2:303344/150‑1 (MQ=255)
                                                                ccccAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCa                                                                                      >  1:3358913/1‑150 (MQ=255)
                                                                       caaAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCa                                                                               <  2:3574323/150‑1 (MQ=255)
                                                                       caaAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCa                                                                               <  1:2583853/150‑1 (MQ=255)
                                                                       caaAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCa                                                                               <  1:4508259/150‑1 (MQ=255)
                                                                       caaAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCa                                                                               <  2:1841220/150‑1 (MQ=255)
                                                                         aaGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGa                                                                             <  2:5392653/150‑1 (MQ=255)
                                                                         aaGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGa                                                                             <  1:3077341/150‑1 (MQ=255)
                                                                            gggCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGc                                                                          <  1:4995016/150‑1 (MQ=255)
                                                                              gCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCAc                                                                        >  1:819022/1‑150 (MQ=255)
                                                                                ttCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGt                                                                      >  2:2544921/1‑150 (MQ=255)
                                                                                ttCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGt                                                                      >  2:5684132/1‑150 (MQ=255)
                                                                                   gAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTAcg                                                                   >  2:1902464/1‑150 (MQ=255)
                                                                                     gCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTAcgcg                                                                 >  1:6617623/1‑150 (MQ=255)
                                                                                      cAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCCGAAGGTGGTCATCGCCCAGATGCACGTACGCgg                                                                >  2:2746234/1‑150 (MQ=255)
                                                                                      cAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCgg                                                                <  2:3867543/150‑1 (MQ=255)
                                                                                      cAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCgg                                                                <  2:6526112/150‑1 (MQ=255)
                                                                                       aaaaGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCggg                                                               >  2:5783437/1‑150 (MQ=255)
                                                                                         aaGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCa                                                                               <  1:6359354/132‑1 (MQ=255)
                                                                                         aaGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGga                                                             <  2:1118049/150‑1 (MQ=255)
                                                                                         aaGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGga                                                             <  2:1372797/150‑1 (MQ=255)
                                                                                             aCCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGgacgac                                                         <  1:95346/150‑1 (MQ=255)
                                                                                                 acaGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGcc                                                     >  1:5545341/1‑150 (MQ=255)
                                                                                                     gACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCccgcc                                                 >  2:1654586/1‑150 (MQ=255)
                                                                                                     gACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCccgcc                                                 >  2:4501567/1‑150 (MQ=255)
                                                                                                       cacGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCgg                                               <  2:1006751/150‑1 (MQ=255)
                                                                                                        acGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGt                                              <  2:4759554/150‑1 (MQ=255)
                                                                                                          gCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTaa                                            >  1:5953944/1‑150 (MQ=255)
                                                                                                           cTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTaaa                                           <  2:6085367/150‑1 (MQ=255)
                                                                                                            tGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGCCGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAg                                          <  2:4306490/150‑1 (MQ=255)
                                                                                                            tGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAg                                          <  1:4887263/150‑1 (MQ=255)
                                                                                                               gAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGgtg                                       >  1:2960664/1‑150 (MQ=255)
                                                                                                                    gACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTg                                  <  2:1871477/150‑1 (MQ=255)
                                                                                                                     aCGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGTATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTgg                                 >  2:3082893/1‑150 (MQ=255)
                                                                                                                     aCGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTgg                                 <  2:4021779/150‑1 (MQ=255)
                                                                                                                     aCGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTgg                                 >  1:1927649/1‑150 (MQ=255)
                                                                                                                     aCGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTATgg                                 >  1:1466289/1‑150 (MQ=255)
                                                                                                                       ggAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGaa                               <  2:6743389/150‑1 (MQ=255)
                                                                                                                              tgatccTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGc                        >  1:2862647/6‑150 (MQ=255)
                                                                                                                              cAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGc                        >  1:6148007/1‑150 (MQ=255)
                                                                                                                              cAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGc                        >  2:6158461/1‑150 (MQ=255)
                                                                                                                               aaGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCa                       <  1:382084/150‑1 (MQ=255)
                                                                                                                               aaGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCa                       >  2:3416041/1‑150 (MQ=255)
                                                                                                                                 gACTGGACGCTGGCCAGCGCCCAGGATCACAGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACCCGGGGACGACGGCCCGCCCGTAAAGGTGTTCTGGAACAAGGGCAAc                     >  1:5493701/1‑150 (MQ=255)
                                                                                                                                  aCTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACa                    <  1:699854/150‑1 (MQ=255)
                                                                                                                                      gACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACATCAc                >  1:135678/1‑150 (MQ=255)
                                                                                                                                          cTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACATCACCCTg            <  1:4895768/150‑1 (MQ=255)
                                                                                                                                          cTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACATCACCCTg            >  2:2814142/1‑150 (MQ=255)
                                                                                                                                              ccAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACATCACCCTGGGGt        >  2:5966442/1‑150 (MQ=255)
                                                                                                                                              ccAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACATCACCCTGGGGt        <  1:2396537/150‑1 (MQ=255)
                                                                                                                                                 gcgcCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACATCACCCTGGGGTTTc     <  2:6628977/150‑1 (MQ=255)
                                                                                                                                                    cccAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACATCACCCTGGGGTTTCGgc  <  2:5152630/150‑1 (MQ=255)
                                                                                                                                                    cccAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACATCACCCTGGGGTTTCGgc  <  2:1457613/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACATCACCCTGGGGTTTCGGC  >  NC_002947/4305916‑4306213

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