Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,667,950 A→C 7.6% I304L (ATC→CTC)  ffh → signal recognition particle GTPase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,667,9500AC 7.6% 380.1 / 22.5 146I304L (ATC→CTC) ffhsignal recognition particle GTPase subunit
Reads supporting (aligned to +/- strand):  ref base A (66/69);  new base C (7/4);  total (73/73)
Fisher's exact test for biased strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.31e-01

Only 100 of 150 total aligned reads displayed.

CCGAGGCGGTGCTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGACAAG  >  NC_002947/1667802‑1668099
                                                                                                                                                    |                                                                                                                                                     
ccGAGGCGGTGCTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGAt                                                                                                                                                      >  1:1185975/1‑150 (MQ=255)
ccGAGGCGGTGCTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGAt                                                                                                                                                      <  1:930866/150‑1 (MQ=255)
 cGAGGCGGTGCTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATc                                                                                                                                                     <  2:763789/150‑1 (MQ=255)
 cGAGGCGGTGCTGCGCGGGCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTCCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATc                                                                                                                                                     <  2:2381433/150‑1 (MQ=255)
 cGAGGCGGGGCTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATc                                                                                                                                                     <  2:774221/150‑1 (MQ=255)
     gCGGTGCTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCACCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGc                                                                                                                                                 <  1:2969985/150‑1 (MQ=255)
      cGGTGCTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCTGTATGGGTGAGAAGACCGAGGCCCTCTAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGATGTGCTCAGCCTGATCGAGCa                                                                                                                                                <  2:562882/150‑1 (MQ=255)
          gctgcGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAgg                                                                                                                                              >  1:298607/1‑148 (MQ=255)
           ctgcGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCg                                                                                                                                           >  2:1226048/1‑150 (MQ=255)
                   tCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGAcc                                                                                                                                   >  2:785805/1‑150 (MQ=255)
                      gTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATc                                                                                                                                >  2:1763259/1‑150 (MQ=255)
                      gTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATc                                                                                                                                <  1:1974879/150‑1 (MQ=255)
                         cGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGAc                                                                                                                             <  2:497942/150‑1 (MQ=255)
                          gTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATc                                                                                                                                <  2:2681081/146‑1 (MQ=255)
                            gCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGACGAGCAGACCATCGAcaag                                                                                                                          >  2:2023441/1‑150 (MQ=255)
                                  aCCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGGATCCTCGGCATGGGCGACGTGCTAAGCCTGCTGGGGGAGGGCGCACAGACCATCGGcaaggccaag                                                                                                                    >  1:2065241/1‑150 (MQ=255)
                                  aCCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGAcaaggccaag                                                                                                                    >  1:2767831/1‑150 (MQ=255)
                                  aCCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGAcaaggccaag                                                                                                                    >  2:2032831/1‑150 (MQ=255)
                                    cGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGAcaaggccaaggc                                                                                                                  <  1:1745358/150‑1 (MQ=255)
                                    cGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGAcaaggccaaggc                                                                                                                  <  2:469448/150‑1 (MQ=255)
                                     ggCACACCCATCAAGTTCATACGTATGGTTGATAAGACCGATGCCCTCGTGTCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACaaggccaaggcc                                                                                                                 <  2:2333701/150‑1 (MQ=255)
                                      gCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCg                                                                                                                >  2:469726/1‑150 (MQ=255)
                                      gCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCg                                                                                                                <  2:1185975/150‑1 (MQ=255)
                                       cAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGa                                                                                                               <  1:2032831/150‑1 (MQ=255)
                                          aCCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCa                                                                                                                                  <  2:1968308/128‑1 (MQ=255)
                                          aCCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCa                                                                                                                                  >  1:1968308/1‑128 (MQ=255)
                                           ccGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAg                                                                                                                                               <  1:2210383/114‑1 (MQ=255)
                                           ccGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAg                                                                                                                                               >  2:2210383/1‑114 (MQ=255)
                                             gATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGTCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCt                                                                                                         >  1:1996914/1‑150 (MQ=255)
                                             gATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAg                                                                                                                                         >  1:2817080/1‑118 (MQ=255)
                                             gATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAg                                                                                                                                         <  2:2817080/118‑1 (MQ=255)
                                              aTCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGa                                                                                                                                                       >  1:1225783/1‑103 (MQ=255)
                                              aTCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTg                                                                                                        <  1:1763259/150‑1 (MQ=255)
                                               tCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCAGACAAGCTgg                                                                                                       >  2:251391/1‑150 (MQ=255)
                                                        tCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAAc                                                                                              <  1:908200/150‑1 (MQ=255)
                                                                 ttGAGAAGACCGATGCCATCGAGCCGTTCCAACCCGACCCCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCGCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaagg                                                                                     <  2:1265437/149‑1 (MQ=255)
                                                                 gTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaagg                                                                                     >  2:2599231/1‑150 (MQ=255)
                                                                 gTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaagg                                                                                     >  2:2674174/1‑150 (MQ=255)
                                                                 gTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaagg                                                                                     >  2:2696643/1‑150 (MQ=255)
                                                                     gaagaCCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGa                                                                                                 >  1:264369/1‑134 (MQ=255)
                                                                      aagaCCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGCCAAGGACAAGGCGGACAAGCTGGCCAAGAAACTGAAGaagggcaag                                                                                >  1:2213237/1‑150 (MQ=255)
                                                                      aagaCCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaagggcaag                                                                                <  2:2245585/150‑1 (MQ=255)
                                                                        gaCCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaagggcaaggg                                                                              >  1:2965287/1‑150 (MQ=255)
                                                                            gAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAGCGGCAAAGCCAACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTc                                                                          <  1:2833954/150‑1 (MQ=255)
                                                                            gAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTc                                                                          >  2:2006391/1‑150 (MQ=255)
                                                                               gCCCTCGGGCCGTTCCACCCCCACCGCGTCGCTTCGCGCATCCTCGGTATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCGGGCCAAGAAACTGAAGTAGGGCAAGGGCTTCGAc                                                                       <  2:1579020/150‑1 (MQ=255)
                                                                               gCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGAc                                                                       <  2:298006/150‑1 (MQ=255)
                                                                                  cTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCGCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTGCGACCTc                                                                    >  2:2584912/1‑150 (MQ=255)
                                                                                  cTCGAGCCGTTCCACACCGACCGCGTCGCGTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTc                                                                    >  1:729245/1‑150 (MQ=255)
                                                                                     gAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGaa                                                                 <  1:1610404/150‑1 (MQ=255)
                                                                                     gAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGaa                                                                 <  1:2102765/150‑1 (MQ=255)
                                                                                     gAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGaa                                                                 <  1:2696643/150‑1 (MQ=255)
                                                                                     gAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGaa                                                                 <  1:737697/150‑1 (MQ=255)
                                                                                     gAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGaa                                                                 <  1:469726/150‑1 (MQ=255)
                                                                                       gCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGa                                                               <  1:2791956/150‑1 (MQ=255)
                                                                                            tCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTcc                                                          <  1:1226048/150‑1 (MQ=255)
                                                                                               aCCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGACGGGCCAGGGCTTCGACCTCGACGACTTCcgcg                                                       >  2:953820/1‑150 (MQ=255)
                                                                                               aCCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCcgcg                                                       >  2:2873545/1‑150 (MQ=255)
                                                                                                ccccGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGAcaaggccaag                                                                                                                    >  2:1668464/1‑88 (MQ=255)
                                                                                                 cccGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAg                                                                >  2:2114093/1‑139 (MQ=255)
                                                                                                 cccGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGAc                                                     <  2:2447673/150‑1 (MQ=255)
                                                                                                 cccGACCGCCTAGATTAGCACATCCGCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGAc                                                     >  2:756049/1‑150 (MQ=255)
                                                                                                      ccGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCt                                                <  1:572604/150‑1 (MQ=255)
                                                                                                         cgTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTgca                                             <  1:2023441/150‑1 (MQ=255)
                                                                                                              cTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCaaa                                        >  2:607210/1‑150 (MQ=255)
                                                                                                              cTTCGCGCATCCTCGGCATGGGCCACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCaaa                                        <  2:2065241/150‑1 (MQ=255)
                                                                                                               ttCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAAt                                       >  1:1967440/1‑150 (MQ=255)
                                                                                                                     cATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATgaagaa                                 <  1:2575990/150‑1 (MQ=255)
                                                                                                                      atcatcGGCATGGGCGACGTGCACAGGCTGATCGAGCAGGCCGAGCAGACCATCGACAAGACCAAAGCCGACAAGCTGGCCAAGAAAGTGAAGAAGGGCAAGGGATTCGACCTCGAAGACTTCCGAGACCAGCTGCAGCAAATgaagaaa                                >  2:879925/1‑149 (MQ=255)
                                                                                                                          tCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAGTGAAGAACATgg                            >  1:3045223/1‑150 (MQ=255)
                                                                                                                          tCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATgg                            >  2:911574/1‑150 (MQ=255)
                                                                                                                            ggCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGGCAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGggc                          <  2:1553073/150‑1 (MQ=255)
                                                                                                                             gCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGggcg                         <  1:911574/150‑1 (MQ=255)
                                                                                                                                 gggCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaagggc                                                                                   >  2:2909208/1‑88 (MQ=255)
                                                                                                                                 gggCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaagggc                                                                                   <  1:2909208/88‑1 (MQ=255)
                                                                                                                                  ggCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTc                    >  2:20089/1‑150 (MQ=255)
                                                                                                                                   gCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTgaa                                                                                          >  1:268211/1‑79 (MQ=255)
                                                                                                                                      ccGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTcggcg                <  1:756049/149‑1 (MQ=255)
                                                                                                                                       tctGCTCGGCCTGCTCGATCAGGCTGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGTAGCAAATGAAGAACATGGGCGGCCTcggcgg               >  2:10829/3‑150 (MQ=255)
                                                                                                                                        gTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCggcggc              >  2:368959/1‑150 (MQ=255)
                                                                                                                                        gTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCggcggc              >  2:376560/1‑150 (MQ=255)
                                                                                                                                           ctcggCCTGCTCGATCAGGCTGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGaa                                                                 <  1:539723/92‑1 (MQ=38)
                                                                                                                                           cTCGGCCTGCTCGATCAGGCTGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTg           >  1:455263/1‑150 (MQ=255)
                                                                                                                                           cTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCt                                                >  2:2779014/1‑113 (MQ=255)
                                                                                                                                            tcggCCTGCTCGATCAGGCTGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTc                                                                          >  1:1412836/4‑86 (MQ=38)
                                                                                                                                            tcggCCTGCTCGATCAGGCTGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATgaag                                   <  1:1770436/122‑1 (MQ=255)
                                                                                                                                            tcggCCTGCTCGATCAGGCTGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATgaag                                   >  2:1770436/4‑125 (MQ=255)
                                                                                                                                            tCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGTCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGa          <  1:879925/150‑1 (MQ=255)
                                                                                                                                            tCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGa          >  2:2699769/1‑150 (MQ=255)
                                                                                                                                            tCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGa          >  2:1971304/1‑150 (MQ=255)
                                                                                                                                              gatcTGATCGAGCAGGCCGAGCAGACCATAGACAAGGCCAAGGCCAACGAGAGGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATg        <  2:1617880/147‑1 (MQ=255)
                                                                                                                                              aGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATg        <  1:2699769/150‑1 (MQ=255)
                                                                                                                                                ccTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGa      >  1:1113819/1‑150 (MQ=255)
                                                                                                                                                ccTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGa      >  1:763986/1‑150 (MQ=255)
                                                                                                                                                ccTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGa      >  2:888871/1‑150 (MQ=255)
                                                                                                                                                ccTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGa      >  1:1796447/1‑150 (MQ=255)
                                                                                                                                                   gATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTcgg                  >  2:536237/1‑135 (MQ=255)
                                                                                                                                                   gATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGACaa   >  2:2873979/1‑150 (MQ=255)
                                                                                                                                                    aTCGAGCAGGCCGAGCAGACCATCGACTGGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGagatc  >  1:1617880/1‑146 (MQ=255)
                                                                                                                                                    aTCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGACAAg  >  1:1376338/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCGAGGCGGTGCTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGACAAG  >  NC_002947/1667802‑1668099

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