Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 806,330 T→C 6.2% K335E (AAG→GAG)  PP_0693 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947806,3300TC 6.2% 292.9 / 6.6 97K335E (AAG→GAG) PP_0693hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (41/49);  new base C (3/3);  total (45/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.49e-01

TTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAG  >  NC_002947/806181‑806468
                                                                                                                                                     |                                                                                                                                          
ttCCAGGTCGCGTACGCGCTCGTCGAGGTCTTAGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCtt                                                                                                                                            >  1:1171559/1‑150 (MQ=255)
  ccAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCCGCTCGTAAAATTCGTCTTTGGTCGTTATCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGa                                                                                                                                          <  1:1498647/150‑1 (MQ=255)
  ccAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGa                                                                                                                                          <  1:2168591/150‑1 (MQ=255)
  ccAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGa                                                                                                                                          <  2:1144325/150‑1 (MQ=255)
    aGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGc                                                                                                                                        >  1:1633990/1‑150 (MQ=255)
      gTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCAc                                                                                                                                      >  1:1866795/1‑150 (MQ=255)
      gTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCAc                                                                                                                                      <  1:1421542/150‑1 (MQ=255)
       tCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGTCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCAcc                                                                                                                                     <  1:2105231/150‑1 (MQ=255)
       tCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCAcc                                                                                                                                     >  1:1762172/1‑150 (MQ=255)
         gcgTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGa                                                                                                                                   >  2:2181910/1‑150 (MQ=255)
         gcgTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGa                                                                                                                                   <  1:437761/150‑1 (MQ=255)
           gTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACa                                                                                                                                 <  2:1633990/150‑1 (MQ=255)
            tACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAg                                                                                                                                >  1:575352/1‑150 (MQ=255)
                   tcgtcgAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCAcc                                                                                                                                     <  2:2365083/138‑1 (MQ=255)
                   tcgtcgAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCAcc                                                                                                                                     >  1:2365083/1‑138 (MQ=255)
                     gtcgAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAg                                                                                                                       >  1:1688291/1‑150 (MQ=255)
                      tcgAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGc                                                                                                                      <  2:1742407/150‑1 (MQ=255)
                             cTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCt                                                                                                               <  2:1086354/150‑1 (MQ=255)
                              ttCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTg                                                                                                              <  1:1633342/150‑1 (MQ=255)
                              ttCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTg                                                                                                              <  1:283528/150‑1 (MQ=255)
                              ttCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTg                                                                                                              >  2:799864/1‑150 (MQ=255)
                              ttCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCACCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTCAAGCACCgct                                                                                                                                  <  2:1195237/130‑3 (MQ=255)
                              ttCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCACCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTCAAGCACCgct                                                                                                                                  >  1:1195237/1‑128 (MQ=255)
                              ttCGGCACTGAACTGCCGTGACTGCTCCAGCGCGCCGGCCAGCTCCTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCTAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGTTTGCCTAGCCAGCGCTg                                                                                                              <  2:728787/150‑1 (MQ=255)
                                cGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTTGGCGTCTTTGGCCGCCAGCAGGTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGa                                                                                                            <  2:1754454/150‑1 (MQ=255)
                                 ggCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTCCTGGAGCACCTACAGCTCGCCTAGCCAGCGCTGGAt                                                                                                           >  2:2071705/1‑150 (MQ=255)
                                      tGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGc                                                                                                                                        <  1:1812569/116‑1 (MQ=255)
                                      tGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGc                                                                                                                                        >  2:1812569/1‑116 (MQ=255)
                                         aCTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTaa                                                                                                   <  1:495878/150‑1 (MQ=255)
                                         aCTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGATTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTaa                                                                                                   <  1:1047444/150‑1 (MQ=255)
                                            gCCGGGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCCGTGCTTGAGCACCGACAGCGCGCCTAGCCAGCGCTGGATACCGGTAAttt                                                                                                <  1:1165113/150‑1 (MQ=255)
                                              cGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCa                                                                                              <  2:962002/150‑1 (MQ=255)
                                               gTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCt                                                                                                               >  1:20053/1‑132 (MQ=255)
                                               gTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCt                                                                                                               <  2:20053/132‑1 (MQ=255)
                                                    tGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTcc                                                                                        >  2:2104195/1‑150 (MQ=255)
                                                      ggCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAg                                                                                      <  2:1688291/150‑1 (MQ=255)
                                                       gCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAgg                                                                                     >  1:1777928/1‑150 (MQ=255)
                                                         cAGCGCGACGGTCAGCTCGTCGTTGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGTTTGAGCATCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCATgcg                                                                                   <  2:482685/150‑1 (MQ=255)
                                                           gcgcgcCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCa                                                                                 <  1:799864/150‑1 (MQ=255)
                                                             gcgcCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATc                                                                               <  1:2181910/150‑1 (MQ=255)
                                                                    ccAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCt                                                                        <  2:925725/150‑1 (MQ=255)
                                                                     cAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCCGTTGCTGCTCCAGTACCTTGACGTGATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCACCTCCTg                                                                       >  1:454933/1‑150 (MQ=255)
                                                                     cAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTg                                                                       >  2:38522/1‑150 (MQ=255)
                                                                     cAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTg                                                                       >  1:1751166/1‑150 (MQ=255)
                                                                      aGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCTTgg                                                                      >  1:545888/1‑150 (MQ=255)
                                                                      aGCTCGTCGTGGGCACCTTTTGCCGCCAGAATTTCCTGCTCCAGTACCTTGACGTCATCGACCCGCCCATGCAGTTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTgg                                                                      <  1:1022416/150‑1 (MQ=255)
                                                                        cTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGt                                                                    <  1:1934780/150‑1 (MQ=255)
                                                                             cgtGGGCGTCTTTGGCCGCCAGCAGTTGCTGCCATCGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTgccgc                                                               <  2:1171559/150‑1 (MQ=255)
                                                                              gtGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTgg                                                                      <  1:1592468/142‑1 (MQ=255)
                                                                              gtGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTgg                                                                      >  2:1592468/1‑142 (MQ=255)
                                                                                 ggCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGAAAACTTTCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACg                                                           >  1:1354706/1‑150 (MQ=255)
                                                                                    gTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGgcgc                                                                                  >  1:595957/1‑124 (MQ=255)
                                                                                    gTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGgcgc                                                                                  <  2:595957/124‑1 (MQ=255)
                                                                                     tCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTgg                                                       <  1:1728412/150‑1 (MQ=255)
                                                                                      cTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGt                                                      >  2:110405/1‑150 (MQ=255)
                                                                                       tttGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTCAAGCACCgct                                                                                                                                  >  2:394662/1‑71 (MQ=255)
                                                                                       tttGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTCAAGCACCgct                                                                                                                                  <  1:394662/73‑3 (MQ=255)
                                                                                          ggCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGTACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGtt                                                  >  1:1237950/1‑150 (MQ=255)
                                                                                          ggCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGtt                                                  <  2:254338/150‑1 (MQ=255)
                                                                                          ggCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGtt                                                  <  2:575352/150‑1 (MQ=255)
                                                                                          ggCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGtt                                                  <  2:1751166/150‑1 (MQ=255)
                                                                                            ccgccAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTa                                                                                                    <  2:1822956/98‑1 (MQ=255)
                                                                                            ccgccAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTa                                                                                                    >  1:1822956/1‑98 (MQ=255)
                                                                                                cagcagTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAg                                            >  2:1045301/1‑150 (MQ=255)
                                                                                                        gctgctCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCt                                    <  1:121366/150‑1 (MQ=255)
                                                                                                            ctCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTc                                <  2:2388963/150‑1 (MQ=255)
                                                                                                             tCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTg                                                        >  1:1754622/1‑125 (MQ=255)
                                                                                                             tCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTg                                                        <  2:1754622/125‑1 (MQ=255)
                                                                                                               cAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTgccgc                                                               >  2:1781493/1‑116 (MQ=255)
                                                                                                               cAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTgccgc                                                               <  1:1781493/116‑1 (MQ=255)
                                                                                                               cAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGcc                             >  2:1868567/1‑150 (MQ=255)
                                                                                                                aGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGtt                                                  >  1:2133609/1‑128 (MQ=255)
                                                                                                                aGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGtt                                                  <  2:2133609/128‑1 (MQ=255)
                                                                                                                 gTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCtt                                                                                                                                            <  1:821874/37‑1 (MQ=255)
                                                                                                                 gTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCtt                                                                                                                                            >  2:821874/1‑37 (MQ=255)
                                                                                                                 gTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCAc                           >  1:858708/1‑150 (MQ=255)
                                                                                                                    ccTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGt                        >  1:1281039/1‑150 (MQ=255)
                                                                                                                     cTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTc                                                 <  1:1991490/124‑1 (MQ=255)
                                                                                                                     cTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTc                                                 >  2:1991490/1‑124 (MQ=255)
                                                                                                                     cTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTa                       >  1:1931952/1‑150 (MQ=255)
                                                                                                                     cTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTa                       >  2:1304390/1‑150 (MQ=255)
                                                                                                                      ttGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGcc                                                                                                                          <  2:164325/50‑1 (MQ=255)
                                                                                                                      ttGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGcc                                                                                                                          >  1:164325/1‑50 (MQ=255)
                                                                                                                      ttGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCGACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCAGGGTAg                      >  1:408459/1‑150 (MQ=255)
                                                                                                                          cGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCt                  <  1:1868567/150‑1 (MQ=255)
                                                                                                                              aacGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGtt                                                  <  1:1520703/112‑1 (MQ=255)
                                                                                                                              aacGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGtt                                                  >  2:1520703/3‑114 (MQ=255)
                                                                                                                               tCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGAc             <  2:1354706/150‑1 (MQ=255)
                                                                                                                                  aCCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTACCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTgg          >  1:2293596/1‑150 (MQ=255)
                                                                                                                                  aCCAGCGCAAAGCGGTGCTTGAGCACAGACTGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTgg          <  1:1211081/150‑1 (MQ=255)
                                                                                                                                   ccAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGa         <  1:110405/150‑1 (MQ=255)
                                                                                                                                   ccAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGATCCTGCTCGCCACGGTAGTGCTCGGACTGGa         <  1:1045301/150‑1 (MQ=255)
                                                                                                                                   ccAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCAGAAACGCCAGCAATACGGCACGCTTACCTGCTCGCCACGGTAGTGCTCGGACTGGa         >  1:541005/1‑150 (MQ=255)
                                                                                                                                    cAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTa                       <  2:391542/135‑1 (MQ=255)
                                                                                                                                    cAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTa                       >  1:391542/1‑135 (MQ=255)
                                                                                                                                      gcgcAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGa      >  2:1875478/1‑150 (MQ=255)
                                                                                                                                         cAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCa   >  1:1078933/1‑150 (MQ=255)
                                                                                                                                          aaaGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAg  >  1:1447961/1‑150 (MQ=255)
                                                                                                                                             tcGGTGCTCAAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCCGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCg                <  2:2378335/132‑1 (MQ=255)
                                                                                                                                             tcGGTGCTCAAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCCGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCg                >  1:2378335/2‑133 (MQ=255)
                                                                                                                                                     |                                                                                                                                          
TTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAG  >  NC_002947/806181‑806468

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