Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,172,911 C→G 5.2% A329G (GCG→GGG)  PP_5407 → TnsD‑like transposition protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,172,9110CG 5.2% 1110.9 / 35.8 328A329G (GCG→GGG) PP_5407TnsD‑like transposition protein
Reads supporting (aligned to +/- strand):  ref base C (155/156);  new base G (9/8);  total (164/164)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

Only 100 of 328 total aligned reads displayed.

CCCAAATGACCCGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAGACGTCTGT  >  NC_002947/6172764‑6173059
                                                                                                                                                   |                                                                                                                                                    
cccAAATGACCCGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGa                                                                                                                                                    <  1:4951995/150‑1 (MQ=255)
cccAAATGACCCGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGa                                                                                                                                                    <  1:6298181/150‑1 (MQ=255)
 ccAAATGACCCGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAg                                                                                                                                                   >  1:479480/1‑150 (MQ=255)
 ccAAATGACCCGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAg                                                                                                                                                   >  2:1088838/1‑150 (MQ=255)
       gACCCGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTTGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTaaa                                                                                                                                             <  1:4103972/150‑1 (MQ=255)
       gACCCGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTaaa                                                                                                                                             <  2:3184795/150‑1 (MQ=255)
          ccGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACtt                                                                                                                                          >  2:2341078/1‑150 (MQ=255)
          ccGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACtt                                                                                                                                          <  1:6217389/150‑1 (MQ=255)
           cGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCGCTGGGAGGCTagtcact                                                                                                                                         <  1:865089/150‑7 (MQ=255)
                  cccGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAaca                                                                                                                                  <  1:2057176/150‑1 (MQ=255)
                   ccGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAacac                                                                                                                                 >  2:5667123/1‑150 (MQ=255)
                      tGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACActtt                                                                                                                              <  2:5812130/150‑1 (MQ=255)
                       ggTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACActttc                                                                                                                             <  2:3368711/150‑1 (MQ=255)
                        gTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACActttct                                                                                                                            <  2:3556594/150‑1 (MQ=255)
                           aCTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACActcattagg                                                                                                                         >  2:4598332/1‑143 (MQ=255)
                           aCTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTa                                                                                                                         >  2:4425076/1‑150 (MQ=255)
                            cTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTaa                                                                                                                        <  1:1228908/150‑1 (MQ=255)
                            cTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTaa                                                                                                                        <  1:2796804/150‑1 (MQ=255)
                            cTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTaa                                                                                                                        <  2:5442687/150‑1 (MQ=255)
                             tGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAg                                                                                                                       >  1:6769335/1‑150 (MQ=255)
                              gCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGc                                                                                                                      <  1:162033/150‑1 (MQ=255)
                                 aTCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTa                                                                                                                   >  2:2075148/1‑150 (MQ=255)
                                        gAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCTAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGc                                                                                                            >  2:3631669/1‑150 (MQ=255)
                                        gAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGc                                                                                                            >  1:3161185/1‑150 (MQ=255)
                                            cATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCaa                                                                                                        <  1:5659203/150‑1 (MQ=255)
                                             aTCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCaag                                                                                                       <  2:6709734/150‑1 (MQ=255)
                                                tGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCaaggaa                                                                                                    >  2:897984/1‑150 (MQ=255)
                                                   gTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTCTCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACTCTTTCTTTAAGCCTATTTCAGCCCAaggaagga                                                                                                 <  2:6528542/150‑1 (MQ=255)
                                                   gTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGACGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTACACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAaggaagga                                                                                                 <  1:309483/150‑1 (MQ=255)
                                                    tCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGgag                                                                                                >  1:4186825/1‑150 (MQ=255)
                                                    tCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGgag                                                                                                >  1:5832344/1‑150 (MQ=255)
                                                    tCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGgag                                                                                                >  2:5069450/1‑150 (MQ=255)
                                                      tGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGaa                                                                                              >  1:3245402/1‑150 (MQ=255)
                                                       gATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAAc                                                                                             >  2:472854/1‑150 (MQ=255)
                                                               ggCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGccc                                                                                     >  2:2984003/1‑150 (MQ=255)
                                                                  tATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGa                                                                                  >  1:1902909/1‑150 (MQ=255)
                                                                    ttCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGaaa                                                                                >  1:2927255/1‑150 (MQ=255)
                                                                      cGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGACAcc                                                                              >  1:4485142/1‑150 (MQ=255)
                                                                      cGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAAcc                                                                              <  1:3384602/150‑1 (MQ=255)
                                                                       gggAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAAccc                                                                             >  2:6046868/1‑150 (MQ=255)
                                                                             ttAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGcccc                                                                                    <  2:1340109/137‑1 (MQ=255)
                                                                             ttAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCaaaaaa                                                                       >  2:2775312/1‑150 (MQ=255)
                                                                             ttAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCaaaaaa                                                                       >  1:2340192/1‑150 (MQ=255)
                                                                             ttAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCaaaaaa                                                                       <  2:3928235/150‑1 (MQ=255)
                                                                             ttAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCaaaaaa                                                                       >  1:605812/1‑150 (MQ=255)
                                                                             ttAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCaaaaaa                                                                       >  2:1539927/1‑150 (MQ=255)
                                                                              tAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAAc                                                                      <  1:2984003/150‑1 (MQ=255)
                                                                               aGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACt                                                                     >  1:6534924/1‑150 (MQ=255)
                                                                               aGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACt                                                                     >  2:4090212/1‑150 (MQ=255)
                                                                                gAATCATTCGTTAAGGCCTATCATCATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCGAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACACAAAAAACTa                                                                    <  2:5952917/150‑1 (MQ=255)
                                                                                gAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTa                                                                    <  1:1396462/150‑1 (MQ=255)
                                                                                 aaTCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTAt                                                                   >  1:2525391/1‑150 (MQ=255)
                                                                                  aTCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTAtt                                                                  <  2:5632023/150‑1 (MQ=255)
                                                                                   tCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTAttt                                                                 >  1:2848833/1‑150 (MQ=255)
                                                                                     aTTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTaa                                                               <  2:2624851/150‑1 (MQ=255)
                                                                                     aTTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTaa                                                               <  2:1953108/150‑1 (MQ=255)
                                                                                     aTTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTaa                                                               <  2:2849255/150‑1 (MQ=255)
                                                                                      ttCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAAt                                                              <  2:6291974/150‑1 (MQ=255)
                                                                                      ttCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAAt                                                              >  2:729244/1‑150 (MQ=255)
                                                                                      ttCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAAt                                                              >  1:3691734/1‑150 (MQ=255)
                                                                                      ttCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAAt                                                              >  1:4421168/1‑150 (MQ=255)
                                                                                        cGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGa                                                            <  2:6160061/150‑1 (MQ=255)
                                                                                         gTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGAc                                                           <  2:2410506/150‑1 (MQ=255)
                                                                                           tAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACtt                                                         >  2:2717001/1‑150 (MQ=255)
                                                                                             aGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGa                                                       <  1:1015571/150‑1 (MQ=255)
                                                                                             aGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGa                                                       >  1:116885/1‑150 (MQ=255)
                                                                                             aGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGa                                                       <  1:1486115/150‑1 (MQ=255)
                                                                                              ggCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGaa                                                      <  2:4173519/150‑1 (MQ=255)
                                                                                               gCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAg                                                     <  1:344706/150‑1 (MQ=255)
                                                                                               gCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAg                                                     <  1:2730770/150‑1 (MQ=255)
                                                                                                 cTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGga                                                   >  2:5219942/1‑150 (MQ=255)
                                                                                                 cTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGga                                                   <  2:2508731/150‑1 (MQ=255)
                                                                                                 cTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGga                                                   >  1:4872122/1‑150 (MQ=255)
                                                                                                    tCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGaa                                                >  1:6603656/1‑150 (MQ=255)
                                                                                                      aaCATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGaaaa                                              <  1:3815239/150‑1 (MQ=255)
                                                                                                          tGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAAtcctc                                         <  2:728012/150‑1 (MQ=255)
                                                                                                           gCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAAtcctc                                         <  1:6406851/150‑1 (MQ=255)
                                                                                                            ccGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAAtcctcc                                        <  1:2423586/150‑1 (MQ=255)
                                                                                                             cgcgAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCcg                                       >  2:1493698/1‑150 (MQ=255)
                                                                                                              gcgAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCcgc                                      <  1:1824051/150‑1 (MQ=255)
                                                                                                                        agcaACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACtgatga                            <  1:1498832/150‑1 (MQ=255)
                                                                                                                         gcaACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGAc                           >  2:6801627/1‑150 (MQ=255)
                                                                                                                             cTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGgc                       <  2:765406/150‑1 (MQ=255)
                                                                                                                             cTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGgc                       <  1:5200928/150‑1 (MQ=255)
                                                                                                                             cTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGgc                       <  2:3356712/150‑1 (MQ=255)
                                                                                                                             cTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGgc                       <  1:6127414/150‑1 (MQ=255)
                                                                                                                                  aGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAg                  >  2:1630187/1‑150 (MQ=255)
                                                                                                                                  aGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACATTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTAGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAg                  >  2:3417138/1‑150 (MQ=255)
                                                                                                                                   gTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGaa                                   >  2:3721294/1‑132 (MQ=255)
                                                                                                                                    agtttAGCCTCGCTGGGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCa                >  1:5444314/5‑150 (MQ=255)
                                                                                                                                     gtttAGCCTCGCTGGGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCaa               >  2:4533757/4‑150 (MQ=255)
                                                                                                                                      tttAGCCTCGCTGGGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCaaa              >  1:102326/3‑150 (MQ=255)
                                                                                                                                      tttAGCCTCGCTGGGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCaaa              >  2:3793756/3‑150 (MQ=255)
                                                                                                                                          gCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAg          >  1:4290369/1‑150 (MQ=255)
                                                                                                                                            cTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAGAc        <  2:2896387/150‑1 (MQ=255)
                                                                                                                                             tCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAGACg       >  2:187098/1‑150 (MQ=255)
                                                                                                                                                cAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAGACGTCt    <  2:479480/150‑1 (MQ=255)
                                                                                                                                                cAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAGACGTCt    <  1:5979641/150‑1 (MQ=255)
                                                                                                                                                  gCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAGACGTCtgt  <  1:363190/150‑1 (MQ=255)
                                                                                                                                                  gCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAGACGTCtgt  <  1:2009265/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCCAAATGACCCGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTTCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAGACGTCTGT  >  NC_002947/6172764‑6173059

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: