Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 280,202 C→G 6.3% L16L (CTG→CTC) ‡ sctS ← sulfur compound ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947280,2020CG 6.3% 574.4 / 21.3 175L16L (CTG→CTC) ‡sctSsulfur compound ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (78/86);  new base G (4/7);  total (82/93)
Fisher's exact test for biased strand distribution p-value = 5.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 181 total aligned reads displayed.

GTGCCGCGGAAGAACGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGCCAACAGTTCCGG  >  NC_002947/280053‑280350
                                                                                                                                                     |                                                                                                                                                    
    cgcgGAAGAACGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGACCAACCTGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGa                                                                                                                                                  <  2:2062574/150‑1 (MQ=255)
    cgcgGAAGAACGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGa                                                                                                                                                  <  2:3324487/150‑1 (MQ=255)
    cgcgGAAGAACGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGa                                                                                                                                                  >  1:338699/1‑150 (MQ=255)
     gcgGAAGAACGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGaa                                                                                                                                                 <  1:210264/150‑1 (MQ=255)
     gcgGAAGAACGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGaa                                                                                                                                                 <  1:863026/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)
             aCGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCCTGCTGAAGGGCGCgg                                                                                                                                         <  2:5209315/150‑1 (MQ=255)
             aCGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCCTGCTGAAGGGCGCgg                                                                                                                                         <  1:3029096/150‑1 (MQ=255)
              cGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGa                                                                                                                                        <  2:2467476/150‑1 (MQ=255)
               gACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAg                                                                                                                                       >  1:124601/1‑150 (MQ=255)
                  acGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCg                                                                                                                                    >  2:3743803/1‑150 (MQ=255)
                   cgatcccccGCGCCAGCCTTTTCAGCGGCATAATCTTCGACAGCCGCATCAGCGCCAGGGCAACGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATACCCCGCGCCCTTCAGCTGGATGGGCGCGGAGTCGa                                                                                                                                   <  1:2521549/144‑1 (MQ=255)
                   cGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGa                                                                                                                                   <  2:4145082/150‑1 (MQ=255)
                   cGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGAAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCTGGATGGGCGCGGAGTCGa                                                                                                                                   <  1:2564886/150‑1 (MQ=255)
                    gTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGAc                                                                                                                                  >  2:2562831/1‑150 (MQ=255)
                      aGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTc                                                                                                                                     <  2:4081190/145‑1 (MQ=255)
                        aCCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACg                                                                                                                              >  1:5010563/1‑150 (MQ=255)
                         cccGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGa                                                                                                                             <  2:2365994/150‑1 (MQ=255)
                          ccGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAg                                                                                                                            >  2:147463/1‑150 (MQ=255)
                          ccGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAg                                                                                                                            >  1:629519/1‑150 (MQ=255)
                            gcgcCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGtt                                                                                                                          <  2:3089867/150‑1 (MQ=255)
                               ccagcctgTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCa                                                                                                                       <  1:2934337/150‑1 (MQ=255)
                               ccagccagTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCCGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCa                                                                                                                       <  2:1849344/150‑1 (MQ=255)
                               ccagccagTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCa                                                                                                                       <  2:603042/150‑1 (MQ=255)
                                cagccagTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAg                                                                                                                      >  1:5195823/1‑150 (MQ=255)
                                 agccagTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAgg                                                                                                                     <  1:423155/150‑1 (MQ=255)
                                 agccagTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAgg                                                                                                                     <  1:4565767/150‑1 (MQ=255)
                                       ttCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCa                                                                                                               <  2:5010563/150‑1 (MQ=255)
                                          aGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCg                                                                                                            >  1:5328527/1‑150 (MQ=255)
                                              gCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCa                                                                                                        >  2:1609076/1‑150 (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)
                                                       tCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGAc                                                                                               >  1:4988475/1‑150 (MQ=255)
                                                           cAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCCTGCTGAAGGGCGCgg                                                                                                                                         >  1:4000574/1‑104 (MQ=38)
                                                           cAGCCGCATCAGCCCCAGGCCAAAGCCCAGCACCACACCGAACAACATCCCGCCGACACTGAGCAGCACCGTAAAACCCGCGCCATACAGCAGCAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCg                                                                                           <  1:3668097/150‑1 (MQ=255)
                                                             gCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGc                                                                                         >  2:4042394/1‑150 (MQ=255)
                                                                   tCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGaa                                                                                   <  2:2947941/150‑1 (MQ=255)
                                                                     aGCGCCAGGGCAAAGCCCATCACCAGACAGAAGAACATGCCGCCGACACTGAGCAACACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGt                                                                                 <  2:4252177/150‑1 (MQ=255)
                                                                     aGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGt                                                                                 >  1:4496941/1‑150 (MQ=255)
                                                                      gcgcCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTa                                                                                <  2:5328527/150‑1 (MQ=255)
                                                                         ccAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACtt                                                                             >  1:4660909/1‑150 (MQ=255)
                                                                          cAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTtc                                                                            >  1:1216301/1‑150 (MQ=255)
                                                                           aGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGgcgc                                                                                                                                           <  2:4569639/86‑1 (MQ=255)
                                                                           aGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGgcgc                                                                                                                                           >  1:4569639/1‑86 (MQ=255)
                                                                           aGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTtct                                                                           <  2:2854865/150‑1 (MQ=255)
                                                                             ggCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCg                                                                         <  2:4660909/150‑1 (MQ=255)
                                                                             ggCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCg                                                                         <  2:4867380/150‑1 (MQ=255)
                                                                             ggCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCg                                                                         <  2:2809434/150‑1 (MQ=255)
                                                                              gCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCgg                                                                        >  2:4784875/1‑150 (MQ=255)
                                                                               cAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGa                                                                       <  1:3444798/150‑1 (MQ=255)
                                                                               cAAAGCCCAGCACCAGACCAAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGa                                                                       >  1:4161895/1‑150 (MQ=255)
                                                                                   gCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGc                                                                   >  1:2285272/1‑150 (MQ=255)
                                                                                    cccAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCt                                                                  >  1:2405292/1‑150 (MQ=255)
                                                                                       aGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAg                                                                                  >  2:2113983/1‑131 (MQ=255)
                                                                                       aGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAg                                                                                  <  1:2113983/131‑1 (MQ=255)
                                                                                         cACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGcc                                                             <  2:5130074/150‑1 (MQ=255)
                                                                                          aCCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGaa                                                                                   <  1:474621/127‑1 (MQ=255)
                                                                                           ccAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGt                                                                                                  >  1:1337057/1‑111 (MQ=255)
                                                                                           ccAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGt                                                                                                  <  2:1337057/111‑1 (MQ=255)
                                                                                            cAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGc                                                          <  2:1239343/150‑1 (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)
                                                                                                     caaCATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCa                                                 <  1:3120858/149‑1 (MQ=255)
                                                                                                          tGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCttt                                            >  2:157293/1‑150 (MQ=255)
                                                                                                            ccgccgACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGa                                          <  2:835810/150‑1 (MQ=255)
                                                                                                                  acacTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGt                                    <  2:4821041/150‑1 (MQ=255)
                                                                                                                      tGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGAt                                >  1:4774989/1‑150 (MQ=255)
                                                                                                                       gAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATg                               <  2:3976767/150‑1 (MQ=255)
                                                                                                                        agcagcACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGtt                                                                                                                          >  2:582280/1‑58 (MQ=255)
                                                                                                                          cagcaCCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGcc                            >  1:3196911/1‑150 (MQ=255)
                                                                                                                           agcaCCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCt                           <  2:2405292/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                                                                >  1:4138688/1‑109 (MQ=255)
                                                                                                                                cGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCttgttg                      >  1:3019534/1‑150 (MQ=255)
                                                                                                                                cGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCttgttg                      >  2:2429827/1‑150 (MQ=255)
                                                                                                                                   ataACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATg                   >  1:1910224/1‑150 (MQ=255)
                                                                                                                                    taACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATgg                  >  2:2840588/1‑150 (MQ=255)
                                                                                                                                    taACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATgg                  <  1:4079858/150‑1 (MQ=255)
                                                                                                                                     aaCCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGc                 >  1:4398754/1‑150 (MQ=255)
                                                                                                                                       tccGCGCCCTTCCTGCTGAAGGGCGCGGGGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTg               >  2:2894526/2‑150 (MQ=255)
                                                                                                                                        ccGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGc              <  2:4289914/150‑1 (MQ=255)
                                                                                                                                        ccGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGc              >  2:690618/1‑150 (MQ=255)
                                                                                                                                         cgcgCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGcc             >  2:4268203/1‑150 (MQ=255)
                                                                                                                                             cccTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGCCAACa         <  1:1609076/150‑1 (MQ=255)
                                                                                                                                              ccTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGCCAACAg        >  2:86914/1‑150 (MQ=255)
                                                                                                                                                ttCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGCCAACAGtt      >  2:418165/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)
                                                                                                                                                     |                                                                                                                                                    
GTGCCGCGGAAGAACGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCAGCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGCCAACAGTTCCGG  >  NC_002947/280053‑280350

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: