Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,367,854 C→A 14.9% S224S (TCG→TCT ftsH ← integral membrane ATP‑dependent zinc metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,367,8540CA14.9% 476.5 / 61.5 173S224S (TCG→TCTftsHintegral membrane ATP‑dependent zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base C (69/78);  new base A (10/16);  total (79/94)
Fisher's exact test for biased strand distribution p-value = 5.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 180 total aligned reads displayed.

TCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAA  >  NC_002947/5367705‑5368003
                                                                                                                                                     |                                                                                                                                                     
tcgtGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTcc                                                                                                                                                       >  2:1736673/1‑150 (MQ=255)
   ttgccgccgCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGaa                                                                                                                                                    <  1:3530286/147‑1 (MQ=255)
    gACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAAc                                                                                                                                                   <  1:445190/150‑1 (MQ=255)
    gACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAAc                                                                                                                                                   <  1:3652638/150‑1 (MQ=255)
       cgccgcCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCgg                                                                                                                                                >  2:2596943/1‑150 (MQ=255)
       cgccgcCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCgg                                                                                                                                                >  2:4950577/1‑150 (MQ=255)
       cgccgcCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCgg                                                                                                                                                >  1:633396/1‑150 (MQ=255)
        gccgccAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGa                                                                                                                                               >  2:2791114/1‑150 (MQ=255)
          cgccCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGaaa                                                                                                                                             <  1:3445244/150‑1 (MQ=255)
             ccATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATgg                                                                                                                                          <  2:2549723/150‑1 (MQ=255)
             ccATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATgg                                                                                                                                          <  2:4324980/150‑1 (MQ=255)
              cATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGt                                                                                                                                         <  1:858809/150‑1 (MQ=255)
              aagaccGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGt                                                                                                                                         <  2:793100/146‑1 (MQ=255)
              aagaccGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGt                                                                                                                                         >  1:3976061/5‑150 (MQ=255)
               aTGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTg                                                                                                                                        >  1:1022562/1‑150 (MQ=255)
                gaccGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTga                                                                                                                                       <  2:3314582/148‑1 (MQ=255)
                 gCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTgaa                                                                                                                                      <  1:2391080/150‑1 (MQ=255)
                 gCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTgaa                                                                                                                                      <  1:2531581/150‑1 (MQ=255)
                  ccGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTgaag                                                                                                                                     >  1:3725085/1‑150 (MQ=255)
                     gcaccACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTGAAGAAc                                                                                                                                  >  1:2542015/4‑150 (MQ=255)
                      cgcCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACg                                                                                                                                 <  2:1465568/150‑1 (MQ=255)
                      cgcCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACg                                                                                                                                 <  2:2448675/150‑1 (MQ=255)
                         cACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTa                                                                                                                              >  2:2022813/1‑150 (MQ=255)
                          aCGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTAc                                                                                                                             >  2:5014065/1‑150 (MQ=255)
                             aTGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACttt                                                                                                                          >  1:1791622/1‑150 (MQ=255)
                               ggCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTGAAGAACGGCACattgg                                                                                                                        <  2:171785/150‑3 (MQ=255)
                               ggCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTGAAGAACGGCACattgg                                                                                                                        <  2:651794/150‑3 (MQ=255)
                                     cAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTttcgcttcgc                                                                                                                  <  1:2155473/150‑1 (MQ=255)
                                     cAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATATCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTttcgcttcgc                                                                                                                  >  2:1233204/1‑150 (MQ=255)
                                       aCGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGccc                                                                                                                >  2:1204282/1‑150 (MQ=255)
                                        cGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCcg                                                                                                               >  1:1661353/1‑150 (MQ=255)
                                           cGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGa                                                                                                            <  1:1475714/150‑1 (MQ=255)
                                           cGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTGAAGAACGGCACATTGGCTTCGCCagcaa                                                                                                            >  1:785225/1‑148 (MQ=255)
                                                   tcgtcgATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTg                                                                                                    <  2:4514282/150‑1 (MQ=255)
                                                    cgtcgATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGATCATGTCGCGCAGACGGCTGGCACCGACACACACGAACATTTCCACGAAGTCCGAACCGGAAATGTTGTAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTgg                                                                                                   <  1:281692/150‑1 (MQ=255)
                                                     gtcgATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGc                                                                                                  >  1:4129704/1‑150 (MQ=255)
                                                     gtcgATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGACCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGc                                                                                                  >  2:3502237/1‑150 (MQ=255)
                                                       cgATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCca                                                                                                >  1:1376413/1‑150 (MQ=255)
                                                       cgATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCca                                                                                                <  1:1456270/150‑1 (MQ=255)
                                                          tGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCcagca                                                                                             <  2:4739005/150‑1 (MQ=255)
                                                            caGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAgg                                                                                           <  1:360793/149‑1 (MQ=255)
                                                             aGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAggg                                                                                          <  1:4522051/150‑1 (MQ=255)
                                                              gatgatGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGt                                                                                         >  1:3512140/1‑150 (MQ=255)
                                                              gatgatGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGt                                                                                         >  1:4715463/1‑150 (MQ=255)
                                                              gatgatGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGt                                                                                         >  2:4403448/1‑150 (MQ=255)
                                                                   tgCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtac                                                                                    <  2:3006371/150‑1 (MQ=255)
                                                                   tgCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtac                                                                                    <  1:4949851/150‑1 (MQ=255)
                                                                      aCGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtaccgg                                                                                 <  1:1428559/150‑1 (MQ=255)
                                                                      aCGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtaccgg                                                                                 >  1:2398501/1‑150 (MQ=255)
                                                                      aCGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtaccgg                                                                                 >  2:338492/1‑150 (MQ=255)
                                                                        ggggCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtaccggtc                                                                               >  2:3646281/1‑149 (MQ=255)
                                                                         gggCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtaccggtac                                                                              >  2:104426/1‑150 (MQ=255)
                                                                         gggCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtaccggtac                                                                              <  2:1145580/150‑1 (MQ=255)
                                                                            cGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtaccggtaccgg                                                                           >  2:3323232/1‑150 (MQ=255)
                                                                            cGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtaccggtaccgg                                                                           <  1:1916443/150‑1 (MQ=255)
                                                                              tGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGt                                                                         <  2:416242/150‑1 (MQ=255)
                                                                              tGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGt                                                                         <  1:3389958/150‑1 (MQ=255)
                                                                                 ttcttGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTggg                                                                      <  2:421296/150‑1 (MQ=255)
                                                                                     tGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGa                                                                  <  2:3672688/150‑1 (MQ=255)
                                                                                         cTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCAt                                                              >  1:2089593/1‑150 (MQ=255)
                                                                                          tGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATc                                                             >  1:1948280/1‑150 (MQ=255)
                                                                                           gCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCa                                                            <  1:2569921/150‑1 (MQ=255)
                                                                                             tCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGc                                                          >  1:4519538/1‑150 (MQ=255)
                                                                                             tCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGc                                                          >  1:2694002/1‑150 (MQ=255)
                                                                                               gAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGcac                                                        >  1:453355/1‑150 (MQ=255)
                                                                                               gAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGcac                                                        <  2:2993755/150‑1 (MQ=255)
                                                                                                 aCATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGCTGGcc                                                                                                       >  2:1524247/1‑101 (MQ=255)
                                                                                                 aCATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGcc                                                                                                       <  1:1524247/101‑1 (MQ=255)
                                                                                                      tCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACcgcgc                                                 <  2:2623926/150‑1 (MQ=255)
                                                                                                        gcgcACACGGCTGGCACCGATACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCgg                                               <  1:4946581/150‑1 (MQ=255)
                                                                                                        gcgcACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGTCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCgg                                               >  2:3822224/1‑150 (MQ=255)
                                                                                                        gcgcACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCgg                                               >  1:3324025/1‑150 (MQ=255)
                                                                                                          gcACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGa                                             <  1:3091246/150‑1 (MQ=255)
                                                                                                              acGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCTCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACg                                         >  2:894829/1‑150 (MQ=255)
                                                                                                              acGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACg                                         >  1:4222226/1‑150 (MQ=255)
                                                                                                               cGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACgg                                        <  1:1136976/150‑1 (MQ=255)
                                                                                                                   tGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCAc                                    >  1:1032888/1‑150 (MQ=255)
                                                                                                                    ggCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCAcc                                   <  2:184773/150‑1 (MQ=255)
                                                                                                                           acacCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGgcgc                            >  2:4567862/1‑150 (MQ=255)
                                                                                                                             acCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTg                          >  1:1205397/1‑150 (MQ=255)
                                                                                                                               ccACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGa                        >  1:4974365/1‑150 (MQ=255)
                                                                                                                               ccACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGa                        >  1:4634556/1‑150 (MQ=255)
                                                                                                                               ccACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGa                        >  1:327936/1‑150 (MQ=255)
                                                                                                                                cACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGaa                       >  2:671120/1‑150 (MQ=255)
                                                                                                                                cACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGaa                       <  2:4035363/150‑1 (MQ=255)
                                                                                                                                 aCGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAAc                      >  1:1691141/1‑150 (MQ=255)
                                                                                                                                   gAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACtt                    >  1:2742133/1‑150 (MQ=255)
                                                                                                                                      cATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGcc                 <  1:1447651/150‑1 (MQ=255)
                                                                                                                                        tttCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCg               >  1:4753191/1‑150 (MQ=255)
                                                                                                                                        tttCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCg               >  2:607794/1‑150 (MQ=255)
                                                                                                                                         ttCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCgg              >  1:1163757/1‑150 (MQ=255)
                                                                                                                                             aCGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATcg          >  1:432459/1‑150 (MQ=255)
                                                                                                                                             aCGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATcg          >  2:1619047/1‑150 (MQ=255)
                                                                                                                                                 aGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCa      <  2:4767669/150‑1 (MQ=255)
                                                                                                                                                  gTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAg     <  1:2225246/150‑1 (MQ=255)
                                                                                                                                                  gTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAg     <  1:1558275/150‑1 (MQ=255)
                                                                                                                                                   tCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAgg    <  2:364114/150‑1 (MQ=255)
                                                                                                                                                     cGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGATCTTGCCCGGATCGCGCAGGaa  <  2:1816501/150‑1 (MQ=255)
                                                                                                                                                     cGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGaa  <  2:5036132/150‑1 (MQ=255)
                                                                                                                                                     cGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGaa  <  1:3391544/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAA  >  NC_002947/5367705‑5368003

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: