Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 280,200 A→C 6.3% L17R (CTG→CGG) ‡ sctS ← sulfur compound ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947280,2000AC 6.3% 558.7 / 21.1 175L17R (CTG→CGG) ‡sctSsulfur compound ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (77/86);  new base C (4/7);  total (82/93)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

Only 100 of 184 total aligned reads displayed.

GCGTGCCGCGGAAGAACGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGCCAACAGTTCCG  >  NC_002947/280051‑280349
                                                                                                                                                     |                                                                                                                                                     
gcgTGCCGCGGAAGAACGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTca                                                                                                                                                       <  2:703925/150‑1 (MQ=255)
gcgTGCCGCGGAAGAACGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTca                                                                                                                                                       <  1:3342105/150‑1 (MQ=255)
 cgTGCCGCGGAAGAACGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTcag                                                                                                                                                      <  1:5108523/150‑1 (MQ=255)
  gTGCCGCGGAAGAACGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTcagc                                                                                                                                                     <  2:4861130/150‑1 (MQ=255)
      cgcgGAAGAACGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGa                                                                                                                                                 <  2:3324487/150‑1 (MQ=255)
       gcgGAAGAACGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGaa                                                                                                                                                <  1:210264/150‑1 (MQ=255)
           aagaaCGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGgc                                                                                                                                            <  1:3895866/150‑1 (MQ=255)
              aaCGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGcgcg                                                                                                                                         <  2:3142862/150‑1 (MQ=255)
              aaCGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGcgcg                                                                                                                                         >  2:2263417/1‑150 (MQ=255)
              aaCGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGcgcg                                                                                                                                         >  2:4079858/1‑150 (MQ=255)
               aCGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCCTGCTGAAGGGCGCgg                                                                                                                                        <  2:2851810/150‑1 (MQ=255)
                cGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGa                                                                                                                                       <  2:2467476/150‑1 (MQ=255)
                cGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGa                                                                                                                                       >  2:1970496/1‑150 (MQ=255)
                 gACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAg                                                                                                                                      >  1:124601/1‑150 (MQ=255)
                  acacGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGt                                                                                                                                     <  2:4592964/150‑1 (MQ=255)
                    acGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCg                                                                                                                                   >  2:3743803/1‑150 (MQ=255)
                     cGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCCTGCTGAAGGGCGCgggttata                                                                                                                                  <  2:4933145/150‑7 (MQ=255)
                      gTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGAc                                                                                                                                 >  2:2562831/1‑150 (MQ=255)
                        aGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTc                                                                                                                                    >  1:4081190/1‑145 (MQ=255)
                          aCCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACg                                                                                                                             >  1:5010563/1‑150 (MQ=255)
                           cccGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGa                                                                                                                            <  2:2365994/150‑1 (MQ=255)
                              gcgcCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGtt                                                                                                                         <  2:3089867/150‑1 (MQ=255)
                                 ccagcctgTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCa                                                                                                                      <  1:1994982/150‑1 (MQ=255)
                                 ccagcctgTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCa                                                                                                                      <  1:2934337/150‑1 (MQ=255)
                                 ccagccagTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCCGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCa                                                                                                                      <  2:1849344/150‑1 (MQ=255)
                                   agccagTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAgg                                                                                                                    <  1:423155/150‑1 (MQ=255)
                                   agccagTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAgg                                                                                                                    <  1:4565767/150‑1 (MQ=255)
                                             gCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGa                                                                                                          >  2:4280391/1‑150 (MQ=255)
                                                gCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCa                                                                                                       >  2:1609076/1‑150 (MQ=255)
                                                  aGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCAtt                                                                                                     <  2:4359824/150‑1 (MQ=255)
                                                   gAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCAttt                                                                                                    <  2:4223399/150‑1 (MQ=255)
                                                    aaTCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTg                                                                                                   <  2:3221642/150‑1 (MQ=255)
                                                      tCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGt                                                                                                 <  2:570099/150‑1 (MQ=255)
                                                       cTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTg                                                                                                <  2:2867484/150‑1 (MQ=255)
                                                        ttCGACAGCTGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGa                                                                                               >  2:2311543/1‑150 (MQ=255)
                                                        ttCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCAttt                                                                                                    >  2:4703235/1‑145 (MQ=255)
                                                         tCGACAGCCGCATCAGCGCCAGGGCACAGCCCAGCACCAGACCGAAGAACATGCCGCCGACCCTGAGCAGCACCGTATAACCCGCGCCCTTCTGCCGGATGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGAc                                                                                              >  1:4988481/1‑150 (MQ=255)
                                                         tCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGAc                                                                                              >  1:3553852/1‑150 (MQ=255)
                                                             cAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCCTGCTGAAGGGCGCgg                                                                                                                                        >  1:4000574/1‑104 (MQ=38)
                                                               gCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGc                                                                                        >  2:4042394/1‑150 (MQ=255)
                                                                     tCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGaa                                                                                  <  2:2947941/150‑1 (MQ=255)
                                                                     tCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCGTCAGCCGGAAGGGCGCGGAGGCGACACCGCGGGGCAGGCTTTCAGCGAGCATGTGGTGACCGCGGCACCGac                                                                                  >  2:4171209/1‑149 (MQ=255)
                                                                       aGCGCCAGGGCAAAGCCCATCACCAGACAGAAGAACATGCCGCCGACACTGAGCAACACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGt                                                                                <  2:4252177/150‑1 (MQ=255)
                                                                       aGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGt                                                                                >  1:4496941/1‑150 (MQ=255)
                                                                           ccAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACtt                                                                            >  1:4660909/1‑150 (MQ=255)
                                                                            cAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTtc                                                                           >  2:1489279/1‑150 (MQ=255)
                                                                            cAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTtc                                                                           >  1:1216301/1‑150 (MQ=255)
                                                                             aGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGgcgc                                                                                                                                          >  1:4569639/1‑86 (MQ=255)
                                                                             aGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTtct                                                                          <  2:127803/150‑1 (MQ=255)
                                                                               ggCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTtctct                                                                        <  2:513546/150‑2 (MQ=255)
                                                                               ggCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCg                                                                        <  2:2073704/150‑1 (MQ=255)
                                                                               ggCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCg                                                                        <  2:2809434/150‑1 (MQ=255)
                                                                                gCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCgg                                                                       >  2:4784875/1‑150 (MQ=255)
                                                                                 cAAAGCCCAGCACCAGACCAAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGa                                                                      >  1:4161895/1‑150 (MQ=255)
                                                                                     gCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGc                                                                  >  1:2285272/1‑150 (MQ=255)
                                                                                      cccAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCt                                                                 >  1:2405292/1‑150 (MQ=255)
                                                                                        cAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTg                                                               <  1:2562831/150‑1 (MQ=255)
                                                                                         aGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAg                                                                                 <  1:2113983/131‑1 (MQ=255)
                                                                                            aCCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGaa                                                                                  <  1:474621/127‑1 (MQ=255)
                                                                                             ccAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGt                                                                                                 >  1:1337057/1‑111 (MQ=255)
                                                                                             ccAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAg                                                          >  2:659174/1‑150 (MQ=255)
                                                                                                gACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGt                                                       >  2:3285490/1‑150 (MQ=255)
                                                                                                   cGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGcc                                                    >  1:864982/1‑150 (MQ=255)
                                                                                                    gaagaaCATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGt                                                                                                 >  2:1895058/1‑104 (MQ=255)
                                                                                                    gaagaaCATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGt                                                                                                 <  1:1895058/104‑1 (MQ=255)
                                                                                                            tGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCttt                                           >  1:2871645/1‑150 (MQ=255)
                                                                                                              ccgccgACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGa                                         <  2:835810/150‑1 (MQ=255)
                                                                                                                    acacTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGt                                   >  2:5280012/1‑150 (MQ=255)
                                                                                                                    acacTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGt                                   <  2:4821041/150‑1 (MQ=255)
                                                                                                                     cacTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTc                                  <  2:589542/150‑1 (MQ=255)
                                                                                                                        tGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGAt                               >  1:4774989/1‑150 (MQ=255)
                                                                                                                         gAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATg                              <  2:3976767/150‑1 (MQ=255)
                                                                                                                         gAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATg                              <  2:1212528/150‑1 (MQ=255)
                                                                                                                          agcagcACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGtt                                                                                                                         <  1:582280/58‑1 (MQ=255)
                                                                                                                             agcaCCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCt                          <  2:2405292/150‑1 (MQ=255)
                                                                                                                              gcaCCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCtt                         <  1:3232007/150‑1 (MQ=255)
                                                                                                                               caCCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCttg                        >  1:485354/1‑150 (MQ=255)
                                                                                                                               caCCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCttg                        >  1:29020/1‑150 (MQ=255)
                                                                                                                                aCCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCttgt                       >  2:5058151/1‑150 (MQ=255)
                                                                                                                                 ccGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTg                                                               <  2:4138688/109‑1 (MQ=255)
                                                                                                                                  cGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCttgttg                     >  1:5308109/1‑150 (MQ=255)
                                                                                                                                  cGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCttgttg                     >  1:3019534/1‑150 (MQ=255)
                                                                                                                                    tcgactccGCGCCCTTCCTGCTGAAGGGCGCGGAGTCGACAACCAGTTGCACGCTTTAAGCGATAAATTGGTGACATCGGCACCGAAGTACATCTCGGACAGCTTGGCCAGCGTGCCATAAGCTTTTAGCTTGTCGATGGCCATGTTTAt                   >  1:339993/7‑150 (MQ=255)
                                                                                                                                      taACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATgg                 <  1:4079858/150‑1 (MQ=255)
                                                                                                                                      taACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATgg                 >  2:2840588/1‑150 (MQ=255)
                                                                                                                                      gactccGCGCCCTTCCTGCTGAAGGGCGCGGGGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGcc                                                    <  1:3141874/111‑1 (MQ=38)
                                                                                                                                       aaCCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGc                >  1:4398754/1‑150 (MQ=255)
                                                                                                                                       aaCCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGc                <  2:2407873/150‑1 (MQ=255)
                                                                                                                                         tccGCGCCCTTCCTGCTGAAGGGCGCGGGGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTg              >  2:2894526/2‑150 (MQ=255)
                                                                                                                                         cccGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTg              >  1:730417/1‑150 (MQ=255)
                                                                                                                                           cgcgCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGcc            >  1:3079851/1‑150 (MQ=255)
                                                                                                                                           cgcgCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGcc            >  2:4268203/1‑150 (MQ=255)
                                                                                                                                           cgcgCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGcc            <  2:1341742/150‑1 (MQ=255)
                                                                                                                                               cccgacagcagGAAGGGCGCGGAGTCGACACCGAGTTGCAGGCTTTCAGCGCTCATTTGGCGACATCGGCACCGAAGTACTTCTCGGACACCGTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGATGATGGCTGCCAACa        <  2:5433514/145‑1 (MQ=255)
                                                                                                                                               cccTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGCCAACa        <  1:1609076/150‑1 (MQ=255)
                                                                                                                                                ccTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGCCAACAg       >  2:86914/1‑150 (MQ=255)
                                                                                                                                                  ttCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGCCAACAGtt     >  2:377229/1‑150 (MQ=255)
                                                                                                                                                  ttCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGCCAACAGtt     >  2:2485499/1‑150 (MQ=255)
                                                                                                                                                   tCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGCCAACAGTTc    <  1:3743803/150‑1 (MQ=255)
                                                                                                                                                     agcagGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGCCAACAGTTCCg  >  1:1713318/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGTGCCGCGGAAGAACGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGCCAACAGTTCCG  >  NC_002947/280051‑280349

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