Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 806,331 G→A 6.2% L334L (CTC→CTT PP_0693 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947806,3310GA 6.2% 288.7 / 7.0 97L334L (CTC→CTTPP_0693hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (42/49);  new base A (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.87e-01

CCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAGCTCTTCCTTGCG  >  NC_002947/806183‑806480
                                                                                                                                                    |                                                                                                                                                     
ccAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCCGCTCGTAAAATTCGTCTTTGGTCGTTATCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGa                                                                                                                                                      <  1:1498647/150‑1 (MQ=255)
ccAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGa                                                                                                                                                      <  2:1144325/150‑1 (MQ=255)
ccAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGa                                                                                                                                                      <  1:2168591/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                                                                                                                          >  2:799864/1‑150 (MQ=255)
                            ttCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTg                                                                                                                          <  1:283528/150‑1 (MQ=255)
                            ttCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTg                                                                                                                          <  1:1633342/150‑1 (MQ=255)
                            ttCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCACCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTCAAGCACCgct                                                                                                                                              >  1:1195237/1‑128 (MQ=255)
                            ttCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCACCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTCAAGCACCgct                                                                                                                                              <  2:1195237/130‑3 (MQ=255)
                            ttCGGCACTGAACTGCCGTGACTGCTCCAGCGCGCCGGCCAGCTCCTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCTAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGTTTGCCTAGCCAGCGCTg                                                                                                                          <  2:728787/150‑1 (MQ=255)
                              cGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTTGGCGTCTTTGGCCGCCAGCAGGTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGa                                                                                                                        <  2:1754454/150‑1 (MQ=255)
                               ggCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTCCTGGAGCACCTACAGCTCGCCTAGCCAGCGCTGGAt                                                                                                                       >  2:2071705/1‑150 (MQ=255)
                                    tGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGc                                                                                                                                                    >  2:1812569/1‑116 (MQ=255)
                                    tGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGc                                                                                                                                                    <  1:1812569/116‑1 (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                                                                                                                           <  2:20053/132‑1 (MQ=255)
                                             gTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCt                                                                                                                           >  1:20053/1‑132 (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                                                                                  >  2:1592468/1‑142 (MQ=255)
                                                                            gtGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTgg                                                                                  <  1:1592468/142‑1 (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                                                                                                                                              <  1:394662/73‑3 (MQ=255)
                                                                                     tttGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTCAAGCACCgct                                                                                                                                              >  2:394662/1‑71 (MQ=255)
                                                                                        ggCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGTACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGtt                                                              >  1:1237950/1‑150 (MQ=255)
                                                                                        ggCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGtt                                                              <  2:575352/150‑1 (MQ=255)
                                                                                        ggCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGtt                                                              <  2:254338/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                                                                           <  1:1781493/116‑1 (MQ=255)
                                                                                                             cAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTgccgc                                                                           >  2:1781493/1‑116 (MQ=255)
                                                                                                             cAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGcc                                         >  2:1868567/1‑150 (MQ=255)
                                                                                                              aGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGtt                                                              <  2:2133609/128‑1 (MQ=255)
                                                                                                              aGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGtt                                                              >  1:2133609/1‑128 (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)
                                                                                                                                                    gAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGc                                          <  2:2023011/110‑1 (MQ=255)
                                                                                                                                                    gAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGc                                          >  1:2023011/1‑110 (MQ=255)
                                                                                                                                                    gAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAGCTCTTCCTTgcg  >  2:1782365/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGCACCGACAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAGCTCTTCCTTGCG  >  NC_002947/806183‑806480

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