Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 543,083 A→G 8.6% V440V (GTA→GTG) ‡ rpoC → DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947543,0830AG 8.6% 482.9 / 24.5 163V440V (GTA→GTG) ‡rpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base A (73/76);  new base G (8/6);  total (81/82)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 166 total aligned reads displayed.

GCGACCACCATCAAGGCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAA  >  NC_002947/542934‑543229
                                                                                                                                                     |                                                                                                                                                  
gCGACCACCATCAAGGCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTa                                                                                                                                                    <  2:1724249/150‑1 (MQ=255)
gCGACCACCATCAAGGCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTa                                                                                                                                                    <  2:4390286/150‑1 (MQ=255)
     caccaTCAAGGCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGAt                                                                                                                                               <  1:2575696/150‑1 (MQ=255)
        caTCAAGGCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGa                                                                                                                                            >  2:1326643/1‑150 (MQ=255)
        caTCAAGGCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGa                                                                                                                                            <  2:5381113/150‑1 (MQ=255)
        caTCAAGGCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGa                                                                                                                                            >  2:535038/1‑150 (MQ=255)
         aTCAAGGCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGaa                                                                                                                                           <  1:4063333/150‑1 (MQ=255)
         aTCAAGGCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGaa                                                                                                                                           <  2:3372169/150‑1 (MQ=255)
            aaGGCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGt                                                                                                                                        >  1:565261/1‑150 (MQ=255)
             aGGCTGCCAAGAAGATGGTGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTaa                                                                                                                                      <  2:2984819/150‑1 (MQ=255)
              ggCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTaa                                                                                                                                      <  2:432226/150‑1 (MQ=255)
              ggCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTaa                                                                                                                                      <  2:3710825/150‑1 (MQ=255)
              ggCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTaa                                                                                                                                      <  2:2484803/150‑1 (MQ=255)
              ggCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTaa                                                                                                                                      <  1:4560977/150‑1 (MQ=255)
              ggCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCACTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTaa                                                                                                                                      <  2:1050204/150‑1 (MQ=255)
                cTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAgg                                                                                                                                    <  1:2245767/150‑1 (MQ=255)
                cTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAgg                                                                                                                                    <  1:1645600/150‑1 (MQ=255)
                cTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAgg                                                                                                                                    <  1:4583059/150‑1 (MQ=255)
                    cAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTAt                                                                                                                                >  2:199609/1‑150 (MQ=255)
                    cAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGAACCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTAt                                                                                                                                <  2:612303/150‑1 (MQ=255)
                     aagaagATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTAtt                                                                                                                               >  2:2641480/1‑150 (MQ=255)
                     aagaagATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTAtt                                                                                                                               <  1:4292537/150‑1 (MQ=255)
                      agaagaTGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTc                                                                                                                              >  2:4978819/1‑150 (MQ=255)
                       gaagaTGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCa                                                                                                                             >  2:666674/1‑150 (MQ=255)
                            tGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGc                                                                                                                        >  1:4844075/1‑150 (MQ=255)
                              gTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCAc                                                                                                                      <  2:1125510/150‑1 (MQ=255)
                                 gAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCg                                                                                                                   >  1:468047/1‑150 (MQ=255)
                                  aGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGc                                                                                                                  <  2:2114693/150‑1 (MQ=255)
                                    cgcgAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTg                                                                                                                <  1:3898103/150‑1 (MQ=255)
                                      cgAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGt                                                                                                              >  1:1252228/1‑150 (MQ=255)
                                               agagGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGtgctgc                                                                                                     <  1:4143358/150‑1 (MQ=255)
                                               agagGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGtgctgc                                                                                                     <  2:2956380/150‑1 (MQ=255)
                                                 agGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGt                                                                                                   <  1:480898/150‑1 (MQ=255)
                                                     gtgGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACaa                                                                                               <  2:933172/150‑1 (MQ=255)
                                                      tgGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAAc                                                                                              >  1:3447820/1‑150 (MQ=255)
                                                         gACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGcc                                                                                           <  1:4284630/150‑1 (MQ=255)
                                                         gACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGcc                                                                                           <  2:1000074/150‑1 (MQ=255)
                                                         gACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGcc                                                                                           <  1:2367640/150‑1 (MQ=255)
                                                         gACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGcc                                                                                           >  2:950073/1‑150 (MQ=255)
                                                          aCGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCg                                                                                          <  1:320535/150‑1 (MQ=255)
                                                            gTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGAc                                                                                        <  2:1152328/150‑1 (MQ=255)
                                                                tcGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCg                                                                                    >  2:4005353/1‑150 (MQ=255)
                                                                     gAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGt                                                                               >  1:3763299/1‑150 (MQ=255)
                                                                            ttCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACCACGCCGACTTCGACGGTGACCAGa                                                                        >  1:2487981/1‑150 (MQ=255)
                                                                             tCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGAt                                                                       <  2:3547321/150‑1 (MQ=255)
                                                                             tCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGAt                                                                       >  2:5443944/1‑150 (MQ=255)
                                                                               gcgAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATgg                                                                     >  2:584507/1‑150 (MQ=255)
                                                                                cgAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGc                                                                    >  1:5207216/1‑150 (MQ=255)
                                                                                      ccccGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCa                                                              <  2:2877739/150‑1 (MQ=255)
                                                                                       cccGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCAc                                                             >  2:2740346/1‑150 (MQ=255)
                                                                                       cccGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCAc                                                             >  1:146504/1‑150 (MQ=255)
                                                                                        ccGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACg                                                            >  2:117429/1‑150 (MQ=255)
                                                                                         cGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGt                                                           >  1:1473806/1‑150 (MQ=255)
                                                                                           tACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGc                                                         <  1:2032114/150‑1 (MQ=255)
                                                                                             cttctgaaCCGTGCGCCTACGCTGCACCGTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTcct                                                       <  2:2783369/144‑2 (MQ=255)
                                                                                             cTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCg                                                       <  2:1252228/150‑1 (MQ=255)
                                                                                                cTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTg                                                          <  1:263985/144‑1 (MQ=255)
                                                                                                 tCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGa                                                   >  1:745893/1‑150 (MQ=255)
                                                                                                  gaaCCGTGCGCCTACGCTGCACCGTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTgtc                                                  <  1:181501/149‑3 (MQ=255)
                                                                                                  gaaCCGTGCGCCTACGCTGCACCGTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTgtc                                                  <  1:2236796/149‑3 (MQ=255)
                                                                                                   aaCCGTGCGCCTACGCTGCACCGTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGAt                                                                       >  2:2044473/1‑128 (MQ=255)
                                                                                                   aaCCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGAcc                                                 <  2:2932088/150‑1 (MQ=255)
                                                                                                     ccGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCt                                               <  2:4645536/150‑1 (MQ=255)
                                                                                                         gCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGaa                                           <  1:3763278/150‑1 (MQ=255)
                                                                                                           aCCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGc                                         >  2:2580462/1‑150 (MQ=255)
                                                                                                           aCCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGc                                         >  1:2065677/1‑150 (MQ=255)
                                                                                                           aCCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGc                                         >  2:4651166/1‑150 (MQ=255)
                                                                                                           aCCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGc                                         >  2:5144729/1‑150 (MQ=255)
                                                                                                            ccGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGcc                                        >  1:1373312/1‑150 (MQ=255)
                                                                                                            ccGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGcc                                        >  2:1773843/1‑150 (MQ=255)
                                                                                                                cgctgcACCGTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGc                                    >  2:1749512/6‑150 (MQ=255)
                                                                                                                  cTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTc                                  >  2:329348/1‑150 (MQ=255)
                                                                                                                    gcACCGTCTTGGTATTCAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGa                                >  2:4482745/2‑150 (MQ=255)
                                                                                                                          tCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAgcgc                           >  2:2836620/1‑149 (MQ=255)
                                                                                                                          tCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAgcgc                           <  1:2836670/149‑1 (MQ=255)
                                                                                                                              ggTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTgcg                      >  2:534798/1‑150 (MQ=255)
                                                                                                                              ggTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTgcg                      >  2:2898145/1‑150 (MQ=255)
                                                                                                                               gTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTgcgc                     <  1:2699985/150‑1 (MQ=255)
                                                                                                                                 attcAGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCtg                   >  2:3452350/4‑150 (MQ=255)
                                                                                                                                  tCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCtga                  <  1:535038/150‑1 (MQ=255)
                                                                                                                                    cAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCtgatg                >  1:5236690/1‑150 (MQ=255)
                                                                                                                                    cAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCtgatg                <  1:2641480/150‑1 (MQ=255)
                                                                                                                                    cAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCtgatg                >  1:2149421/1‑150 (MQ=255)
                                                                                                                                    cAGGCGTTTGAACCCGTACTGATCGAAGCTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAgcgc                           <  2:325619/139‑1 (MQ=255)
                                                                                                                                    cAGGCGTTTGAACCCGTACTGATCGAAGCTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAgcgc                           >  1:325619/1‑139 (MQ=255)
                                                                                                                                     aGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGTCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCtgatga               >  1:2904274/1‑150 (MQ=255)
                                                                                                                                     aGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCtgatga               >  1:4445135/1‑150 (MQ=255)
                                                                                                                                     aGGCATTCGAACCCACGCTGATCGAAGGCAAGGCCATTCAGCTGCATCCGCTCGTCTGTACAGCGTTCAACGCTGACTTCGACGGTGACCAGATGGCTGTTCACGTTCCTCTGTCGCTTGAAGCTCAGCTTGAAGCACGCGTGCtgatga               >  2:3516309/1‑150 (MQ=255)
                                                                                                                                      ggCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACg                                                                                 <  1:2060979/83‑1 (MQ=255)
                                                                                                                                      ggCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACg                                                                                 >  2:2060979/1‑83 (MQ=255)
                                                                                                                                         gTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTc           >  1:5431610/1‑150 (MQ=255)
                                                                                                                                           ttGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGa         <  2:4198731/150‑1 (MQ=255)
                                                                                                                                           ttGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGa         <  1:183891/150‑1 (MQ=255)
                                                                                                                                               aCCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACcaa     <  1:50691/150‑1 (MQ=255)
                                                                                                                                               aCCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACcaa     >  1:2483320/1‑150 (MQ=255)
                                                                                                                                               aCCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACcaa     >  2:949120/1‑150 (MQ=255)
                                                                                                                                                 ccGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAg                                                                         >  2:4024463/1‑80 (MQ=255)
                                                                                                                                                 ccGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAg                                                                         <  1:4024463/80‑1 (MQ=255)
                                                                                                                                                  cGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACcaacaa  >  1:3188364/1‑150 (MQ=255)
                                                                                                                                                   gTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGaa                               >  2:4621075/1‑120 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCGACCACCATCAAGGCTGCCAAGAAGATGGTCGAGCGCGAGCTGCCAGAGGTGTGGGACGTTCTCGCTGAAGTGATTCGCGAACACCCCGTACTGCTCAACCGTGCACCGACCCTTCACCGTCTGGGTATCCAGGCGTTTGAACCCGTACTGATCGAAGGTAAGGCTATTCAGCTGCACCCGCTGGTCTGTGCTGCGTACAACGCCGACTTCGACGGTGACCAGATGGCCGTTCACGTGCCGCTGACCCTGGAAGCCCAGCTCGAAGCGCGTGCGCTGATGATGTCGACCAACAA  >  NC_002947/542934‑543229

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: