Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,384,143 T→C 5.7% L117P (CTG→CCG) ‡
A63A (GCA→GCG
PP_5639 →
PP_3855 ←
hypothetical protein
hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,384,1430TC 5.7% 789.3 / 22.6 243L117P (CTG→CCG) ‡
A63A (GCA→GCG
PP_5639
PP_3855
hypothetical protein
hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (133/96);  new base C (7/7);  total (140/103)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 253 total aligned reads displayed.

AGACCCATACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCTGTTGTTTTCG  >  NC_002947/4383994‑4384290
                                                                                                                                                     |                                                                                                                                                   
aGACCCATACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGctct                                                                                                                                                     >  2:2208609/1‑150 (MQ=255)
  aCCCATACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGc                                                                                                                                                   <  2:3902423/150‑1 (MQ=255)
  aCCCATACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGc                                                                                                                                                   <  1:708042/150‑1 (MQ=255)
    ccATACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCg                                                                                                                                                 >  1:1340116/1‑150 (MQ=255)
      aTACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCAGAGCCTGCGCCGGa                                                                                                                                               <  1:2232596/150‑1 (MQ=255)
      aTACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTAAGGCTCCGGCAGAGCCTGCGCCGGa                                                                                                                                               <  2:2314025/150‑1 (MQ=255)
          gagaTTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGcctgcaggc                                                                                                                                           >  2:4456935/1‑143 (MQ=255)
          gagaTTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCt                                                                                                                                           >  1:4874909/1‑150 (MQ=255)
            gaTTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTc                                                                                                                                         <  1:2185886/150‑1 (MQ=255)
            gaTTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCAGAGCCTGCGCCGGAGCCtga                                                                                                                                         <  1:78529/150‑3 (MQ=255)
            gaTTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCAGAGCCTGCGCCGGAGCCtga                                                                                                                                         <  1:2329938/150‑3 (MQ=255)
             aTTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCg                                                                                                                                        <  1:4467229/150‑1 (MQ=255)
               tCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGGTCGGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCa                                                                                                                                      <  2:5445095/150‑1 (MQ=255)
                  gATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATg                                                                                                                                   >  2:3249275/1‑150 (MQ=255)
                   aTCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGa                                                                                                                                  >  2:1346067/1‑150 (MQ=255)
                      gACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCg                                                                                                                               >  2:1832860/1‑150 (MQ=255)
                      gACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCg                                                                                                                               >  2:3145642/1‑150 (MQ=255)
                         aTGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGaca                                                                                                                            <  2:3981362/150‑1 (MQ=255)
                          tGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGacac                                                                                                                           >  2:4933667/1‑150 (MQ=255)
                              tGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTa                                                                                                                       >  2:2389458/1‑150 (MQ=255)
                              tGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTa                                                                                                                       >  2:2118340/1‑150 (MQ=255)
                              tGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTa                                                                                                                       <  1:338424/150‑1 (MQ=255)
                              tGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTa                                                                                                                       >  1:1951017/1‑150 (MQ=255)
                                    tGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTg                                                                                                                 >  1:3707755/1‑150 (MQ=255)
                                       cGTTGCTGACTCAATTGTTCTGCCCTACACCAGCTATCAGCTTTACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCTACACCTTAGCGCTGTaa                                                                                                              <  2:2730734/150‑1 (MQ=255)
                                           gCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCt                                                                                                          <  1:5352641/150‑1 (MQ=255)
                                            cTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCtt                                                                                                         <  2:2991599/150‑1 (MQ=255)
                                             tGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTc                                                                                                        >  2:4726091/1‑150 (MQ=255)
                                              gACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGGTACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCa                                                                                                       <  1:4590652/150‑1 (MQ=255)
                                              gACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCa                                                                                                       <  1:2708679/150‑1 (MQ=255)
                                               acacAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCt                                                                                                          >  1:335948/1‑146 (MQ=255)
                                               acacAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCt                                                                                                          <  2:335948/146‑1 (MQ=255)
                                                 acaATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTa                                                                                                    >  1:890384/1‑150 (MQ=255)
                                                 acaATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTa                                                                                                    >  2:5216138/1‑150 (MQ=255)
                                                   aaTTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTAcc                                                                                                  >  2:4871648/1‑150 (MQ=255)
                                                   aaTTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTAcc                                                                                                  >  2:4379681/1‑150 (MQ=255)
                                                         tCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCg                                                                                            <  2:549329/150‑1 (MQ=255)
                                                         tCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCg                                                                                            <  1:1005040/150‑1 (MQ=255)
                                                          cTGCCCTACACCAGCTCTCATCAATACAAGCTAGGCAATAGGGTGATACCAAGCGTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTCGCGCGGTACTCCTTCAGTTCCTGGTCgg                                                                                           <  2:3518541/150‑1 (MQ=255)
                                                           tGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGt                                                                                          <  2:641534/150‑1 (MQ=255)
                                                            gCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAg                                                                                                                      <  1:135242/121‑1 (MQ=255)
                                                            gCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGtt                                                                                         <  2:1904095/150‑1 (MQ=255)
                                                              ccTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTaa                                                                                       <  1:1211302/150‑1 (MQ=255)
                                                               cTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAAt                                                                                      <  1:105710/150‑1 (MQ=255)
                                                               cTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAAt                                                                                      >  1:919576/1‑150 (MQ=255)
                                                               cTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAAt                                                                                      >  2:4093545/1‑150 (MQ=255)
                                                               cTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAAt                                                                                      <  1:4025894/150‑1 (MQ=255)
                                                               cTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGCCGGTTAAt                                                                                      >  1:658502/1‑150 (MQ=255)
                                                                     cAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTgc                                                                                >  1:2150408/1‑150 (MQ=255)
                                                                     cAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTgc                                                                                >  2:601951/1‑150 (MQ=255)
                                                                     cAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTgc                                                                                >  1:3214097/1‑150 (MQ=255)
                                                                      aGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGt                                                                                              <  2:1954968/135‑1 (MQ=255)
                                                                        cTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCa                                                                             <  2:4317542/150‑1 (MQ=255)
                                                                         tATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCGTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTACTCCTTCAGTACCTGGTCGGGTCATCCCTGCGCAg                                                                            >  2:3716324/1‑150 (MQ=255)
                                                                           tCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAGTGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTc                                                                          >  2:3315739/1‑150 (MQ=255)
                                                                           tCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTc                                                                          >  2:2042132/1‑150 (MQ=255)
                                                                             aGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGAt                                                                                                                                 <  2:3115335/93‑1 (MQ=255)
                                                                             aGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACAGCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCgg                                                                        >  2:3840966/1‑150 (MQ=255)
                                                                             aGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCt                                                                                                          >  2:3376268/1‑116 (MQ=255)
                                                                             aGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTctg                                                                        >  1:3596654/1‑148 (MQ=255)
                                                                              gCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGAAGGCTCTGCCGGAGCCTACGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCggg                                                                       >  1:1733106/1‑150 (MQ=255)
                                                                                aaTACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGaa                                                                     >  2:1017628/1‑150 (MQ=255)
                                                                                aaTACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGaa                                                                     <  1:246941/150‑1 (MQ=255)
                                                                                  tACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACt                                                                   <  2:862779/150‑1 (MQ=255)
                                                                                     aaGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGAc                                                                >  2:1271906/1‑150 (MQ=255)
                                                                                        cTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCa                                                             <  2:254833/150‑1 (MQ=255)
                                                                                        cTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCa                                                             >  2:2558007/1‑150 (MQ=255)
                                                                                        cTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCa                                                             <  2:1176173/150‑1 (MQ=255)
                                                                                           ggCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGa                                                          >  1:4723148/1‑150 (MQ=255)
                                                                                             cAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATc                                                        >  1:3912531/1‑150 (MQ=255)
                                                                                                 atGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATcgcgc                                                    <  1:3477801/150‑1 (MQ=255)
                                                                                                   ggTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCgg                                                  <  1:524121/150‑1 (MQ=255)
                                                                                                        tACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACg                                             >  2:4267328/1‑150 (MQ=255)
                                                                                                           cAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGt                                          >  1:351838/1‑150 (MQ=255)
                                                                                                               cTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCa                                      >  1:1223365/1‑150 (MQ=255)
                                                                                                               cTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCa                                      >  1:4071111/1‑150 (MQ=255)
                                                                                                                ttCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAg                                     <  1:3776249/150‑1 (MQ=255)
                                                                                                                  cAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGACCTGACTCAGGATCGCGCGGGCACGGGTATCAGcc                                   >  2:1460841/1‑150 (MQ=255)
                                                                                                                     cccGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGccccc                                >  1:4275220/1‑150 (MQ=255)
                                                                                                                      ccGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGACCTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCcct                               >  2:936896/1‑150 (MQ=255)
                                                                                                                      ccGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCcct                               >  1:5080537/1‑150 (MQ=255)
                                                                                                                            aTGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCtgtg                         >  1:539594/1‑150 (MQ=255)
                                                                                                                              gCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTgga                       <  1:3249275/150‑1 (MQ=255)
                                                                                                                              gCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTgga                       <  1:3019534/150‑1 (MQ=255)
                                                                                                                               ccTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTggag                      >  1:1527373/1‑150 (MQ=255)
                                                                                                                                cTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTggagc                     >  1:2771925/1‑150 (MQ=255)
                                                                                                                                cTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTggagc                     >  1:2637711/1‑150 (MQ=255)
                                                                                                                                     gCTCCGGCAGAGCCTGCGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTggagcggagc                >  2:4397276/1‑150 (MQ=255)
                                                                                                                                     gCTCCGGCAGAGCCTGCGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTggagcggagc                >  1:5545755/1‑150 (MQ=255)
                                                                                                                                     gCTCCGGCAGAGCCTGCGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCAGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTggagcggagc                >  1:462692/1‑150 (MQ=255)
                                                                                                                                      cTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTCGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGgagcggagcg               <  2:1776598/150‑1 (MQ=255)
                                                                                                                                       tCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGagcggagcgg              >  1:4384699/1‑150 (MQ=255)
                                                                                                                                       tCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGagcggagcgg              >  1:1477266/1‑150 (MQ=255)
                                                                                                                                       tCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGagcggagcgg              >  1:4118729/1‑150 (MQ=255)
                                                                                                                                         cGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGa                                                                 <  2:4643885/97‑1 (MQ=255)
                                                                                                                                                 gCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCTGTTGtttt    <  2:84230/150‑1 (MQ=255)
                                                                                                                                                  ctctGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCTGTTGTTTTc   <  2:4897364/150‑1 (MQ=255)
                                                                                                                                                   tctGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGCGCGGCTGTTGTTTTCg  >  1:4942900/1‑150 (MQ=255)
                                                                                                                                                   tctGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCTGTTGTTTTCg  >  2:2929726/1‑150 (MQ=255)
                                                                                                                                                   tctGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCAGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCTGTTGTTTTCg  >  2:1712935/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGACCCATACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCTGTTGTTTTCG  >  NC_002947/4383994‑4384290

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: