Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,021,757 C→A 5.8% intergenic (‑52/‑92) PP_5272 ← / → PP_5273 hypothetical protein/oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,021,7570CA 5.8% 378.8 / 6.6 119intergenic (‑52/‑92)PP_5272/PP_5273hypothetical protein/oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (59/53);  new base A (3/4);  total (62/57)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

Only 100 of 120 total aligned reads displayed.

GGGCCGCAGCAGTCGTGGCAGCCGGGCTCGCACTCGAAAGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCGCAGCACGCCGCCTCCTATTACGCCGCCAGCAGCGCGCCACACC  >  NC_002947/6021617‑6021906
                                                                                                                                            |                                                                                                                                                     
gggCCGCAGCAGTCGTGGCAGCCGGGCTCGCACTCGAAAGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGTATTATAGGCTTATTCGCGGGTTCACCCGCGat                                                                                                                                              <  2:869803/150‑2 (MQ=255)
gggCCGCAGCAGTCGTGGCAGCCGGGCTCGCACTCGAAAGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTGAACCCGCGaa                                                                                                                                              <  1:1497294/150‑1 (MQ=255)
   ccGCAGCAGTCGTGGCAGCCGGGCTCGCACTCGATCGATGCGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGa                                                                                                                                           <  1:1307288/150‑1 (MQ=255)
   ccGCAGCAGTCGTGGCAGCCGGGCTCGCACTCGAAAGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGa                                                                                                                                           <  1:497512/150‑1 (MQ=255)
   ccGCAGCAGTCGTGGCAGCCGGGCTCGCACTCGAAAGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGa                                                                                                                                           <  1:2433366/150‑1 (MQ=255)
   ccGCAGCAGTCGTGGCAGCCGGGCTCGCACTCGAAAGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGa                                                                                                                                           <  1:980572/150‑1 (MQ=255)
   ccGCAGCAGTCGTGGCAGCCGGGCTCGCACTCGAAAGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTGAACCCGCGaataa                                                                                                                                           <  1:2178806/150‑4 (MQ=255)
    cGCAGCAGTCGTGGCAGCCGGGCTCGCACTCGAAAGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTGAACCCGCGaataag                                                                                                                                          >  1:858058/1‑146 (MQ=255)
      cagcagTCGTGGCAGCCGGGCTCGCACTCGAAAGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTa                                                                                                                                        <  1:2229455/150‑1 (MQ=255)
       agcagTCGTGGCAGCCGGGCTCGCACTCGAAAGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTAt                                                                                                                                       >  2:2033956/1‑150 (MQ=255)
       agcagTCGTGGCAGCCGGGCTCGCACTCGAAAGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTGGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTAt                                                                                                                                       <  1:1467807/150‑1 (MQ=255)
       agcagTCGTGGCAGACGGGCTCGCCTTAGAACGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTAt                                                                                                                                       <  2:858157/150‑1 (MQ=255)
           gTCGTGGCAGCCGGGCTCGCACTCGAAAGAGGGGATGGGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAg                                                                                                                                   <  2:622403/150‑1 (MQ=255)
               tGGCAGCCGGGCTCGCACTCGAAAGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCt                                                                                                                               >  1:2416998/1‑150 (MQ=255)
                    gCCGGGCTCGCACTCGAAAGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTg                                                                                                                          <  2:1889848/150‑1 (MQ=255)
                      cGGGCTCGCACTCGAAAGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGt                                                                                                                        >  1:1835595/1‑150 (MQ=255)
                       gggCTCGCACTCGAAAGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGtt                                                                                                                       <  1:2345830/150‑1 (MQ=255)
                       gggCTCGCACTCGAAAGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGtt                                                                                                                       <  2:2416998/150‑1 (MQ=255)
                           tCGCACTCGAAAGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTNccc                                                                                                                   <  2:1506573/150‑1 (MQ=255)
                            cGCACTCGAAAGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGaa                                                                                                                                    >  2:2114354/1‑132 (MQ=255)
                               aCTCGAAAGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGATTGTTCACGAGGTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGcc                                                                                                               >  2:64045/1‑150 (MQ=255)
                                 tCGAAAGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCcgc                                                                                                             >  1:1545604/1‑150 (MQ=255)
                                 tCGAAAGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGAAGGGGTACAAGCTGACCAAGTGATTGGATTATAGGCTTATTCGCGGGGTCACCCGCGAACGAGTATTAAGCCCTCGCTGGGTGCCCGGCcgc                                                                                                             >  1:1556196/1‑150 (MQ=255)
                                   gAAAGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGa                                                                                                           <  2:1949197/150‑1 (MQ=255)
                                    aaaGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGAt                                                                                                          <  2:1809544/150‑1 (MQ=255)
                                    aaaGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGAt                                                                                                          <  1:2136547/150‑1 (MQ=255)
                                     aaGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATg                                                                                                         <  2:1926644/150‑1 (MQ=255)
                                           ggATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAg                                                                                                   <  2:1835595/150‑1 (MQ=255)
                                              tGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGt                                                                                                                        >  1:325098/1‑126 (MQ=255)
                                              tGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGt                                                                                                                        <  2:325098/126‑1 (MQ=255)
                                                 gTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATc                                                                                             >  2:164120/1‑150 (MQ=255)
                                                  ttCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATcc                                                                                            <  2:2402428/150‑1 (MQ=255)
                                                  ttCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATcc                                                                                            <  1:991066/150‑1 (MQ=255)
                                                  ttCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATcc                                                                                            <  1:195908/150‑1 (MQ=255)
                                                      cgcAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGt                                                                                        <  1:296992/150‑1 (MQ=255)
                                                       gcAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTc                                                                                       >  1:327050/1‑150 (MQ=255)
                                                         aGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTccc                                                                                     <  1:569512/150‑1 (MQ=255)
                                                            aaGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTcccccc                                                                                  >  2:543076/1‑150 (MQ=255)
                                                              gTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGccc                                                                                                                   >  1:1600469/1‑115 (MQ=255)
                                                              gTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGccc                                                                                                                   <  2:1600469/115‑1 (MQ=255)
                                                                     ttGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCAttttt                                                                         >  2:488977/1‑150 (MQ=255)
                                                                     ttGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTGAACCcgc                                                                                                                                                 <  1:416052/78‑1 (MQ=255)
                                                                       gtgGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGa                                                                                                                                           >  2:1054407/1‑82 (MQ=255)
                                                                       gtgGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGa                                                                                                                                           <  1:1054407/82‑1 (MQ=255)
                                                                        tgGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAg                                                                      <  2:327050/150‑1 (MQ=255)
                                                                         ggCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAgg                                                                     <  1:736745/150‑1 (MQ=255)
                                                                         ggCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAgg                                                                     <  1:2214411/150‑1 (MQ=255)
                                                                          gCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAggg                                                                    >  1:431360/1‑150 (MQ=255)
                                                                                  aGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCa                                                            >  1:1835609/1‑150 (MQ=255)
                                                                                    gACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGt                                                                                                                                         >  2:468113/1‑71 (MQ=255)
                                                                                       aTGACGGGATACAAGCTGACAAAGCAGGTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGAt                                                       >  1:1520584/1‑150 (MQ=255)
                                                                                         gACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATcc                                                     >  2:428347/1‑150 (MQ=255)
                                                                                         gACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATcc                                                     >  2:402472/1‑150 (MQ=255)
                                                                                          aCGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCca                                                    >  1:210935/1‑150 (MQ=255)
                                                                                               aTACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCcaca                                                  <  1:2022526/147‑1 (MQ=255)
                                                                                               aTACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCcaca                                                  >  2:2022526/1‑147 (MQ=255)
                                                                                                   aaGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAg                                                                      >  1:102962/1‑123 (MQ=255)
                                                                                                   aaGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAg                                                                      <  2:102962/123‑1 (MQ=255)
                                                                                                   aaGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCgca                                           <  1:488977/150‑1 (MQ=255)
                                                                                                   aaGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCgca                                           >  1:1518377/1‑150 (MQ=255)
                                                                                                     gCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCgcagc                                         >  1:539956/1‑150 (MQ=255)
                                                                                                     gCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCgcagc                                         >  2:1583304/1‑150 (MQ=255)
                                                                                                       tGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCCATGATCCACAGCGCGCAGCAc                                       >  1:1561634/1‑150 (MQ=255)
                                                                                                       tGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCGCAGCAc                                       >  1:2208851/1‑150 (MQ=255)
                                                                                                        gACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCGCAGCAcg                                      >  1:2016555/1‑150 (MQ=255)
                                                                                                        gACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCGCAGCAcg                                      >  2:694375/1‑150 (MQ=255)
                                                                                                         aCAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGt                                                                                                                                         >  2:1307598/1‑50 (MQ=255)
                                                                                                         aCAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGt                                                                                                                                         <  1:1307598/50‑1 (MQ=255)
                                                                                                         aCAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCGCAGCAcgc                                     <  1:694375/150‑1 (MQ=255)
                                                                                                          cAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCGCAGCAcgcc                                    <  1:2033956/150‑1 (MQ=255)
                                                                                                              gTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCGCAGCACGCCGcct                                >  2:1871092/1‑150 (MQ=255)
                                                                                                                  ttGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGACGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCGCAGCACGCCGCCTCCTa                            >  2:1074654/1‑150 (MQ=255)
                                                                                                                    ggATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCAttt                                                                           <  2:172242/101‑1 (MQ=255)
                                                                                                                    ggATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCAttt                                                                           >  1:172242/1‑101 (MQ=255)
                                                                                                                     gATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCGCAGCACGCCGCCTCCTATTa                         <  1:1583304/150‑1 (MQ=255)
                                                                                                                      aTTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCGCAGCACGCCGCCTCCTATTAc                        >  1:1958752/1‑150 (MQ=255)
                                                                                                                          taGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCGCAGCACGCCGCCTCCTATTAcgccg                    >  1:171351/1‑150 (MQ=255)
                                                                                                                           aGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCGCAGCACGCCGCCTCCTATTAcgccgc                   <  2:1491229/150‑1 (MQ=255)
                                                                                                                            ggCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCGCAGCACGCCGCCTCTTATTACgccgcc                  <  2:1556196/150‑1 (MQ=255)
                                                                                                                            ggCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCGCAGCACGCCGCCTCCTATTACgccgcc                  <  2:1545604/150‑1 (MQ=255)
                                                                                                                             gCTTATTCGCGGGTTCACCCGCGAACTAGTATGAAGCCCTCGCTGGTTCCGCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCGCAGCACGCCGCCCCCTATTACgccgccc                 >  1:551434/1‑149 (MQ=255)
                                                                                                                             gCTTATTCGCGGGTTCACCCGCGAACGAGTATGACGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCGCAGCACGCCGCCTCCTATTACGCCGCca                 <  1:1074654/150‑1 (MQ=255)
                                                                                                                             gCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATg                                                                                                         >  1:1085009/1‑62 (MQ=255)
                                                                                                                             gCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATg                                                                                                         <  2:1085009/62‑1 (MQ=255)
                                                                                                                             gCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGc                                                                                                  <  1:1261157/69‑1 (MQ=255)
                                                                                                                             gCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGc                                                                                                  >  2:1261157/1‑69 (MQ=255)
                                                                                                                             gCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCa                                                                       >  1:1645717/1‑96 (MQ=255)
                                                                                                                             gCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCa                                                                       <  2:1645717/96‑1 (MQ=255)
                                                                                                                             gCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCGCAGCACGCCGCCTCCTATTACgccgccc                 >  2:1550082/1‑149 (MQ=255)
                                                                                                                             gCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCGCAGCACGCCGCCTCCTATTACGCCGCca                 >  1:105254/1‑150 (MQ=255)
                                                                                                                             gCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCGCAGCACGCCGCCTCCTATTACGCCGCca                 >  1:1816545/1‑150 (MQ=255)
                                                                                                                             gCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTGTTCAGGGCCGCCCCAATGATCCACAGCGCGCAACACGCCGCCTCCTGTTACGCCGCca                 >  1:2112102/1‑150 (MQ=255)
                                                                                                                               ttATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCATAGCGCGCAGCACGCCGCCTCCTATTACGCCGCcagc               >  2:1900501/1‑150 (MQ=255)
                                                                                                                                  ttCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCGCAGCACGCCGCCTCCTATTACGCCGCCagcagc            >  2:1152807/1‑150 (MQ=255)
                                                                                                                                    cgcgGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCGCAGCACGCCGCCTCCTATTACGCCGCCAGCAgcgc          <  2:2016555/150‑1 (MQ=255)
                                                                                                                                     gcgGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCGCAGCACGCCGCCTCCTATTACGCCGCCAGCAgcgcg         >  2:497915/1‑150 (MQ=255)
                                                                                                                                     gcgGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCGCAGCACGCCGCCTCCTATTACGCCGCCAGCAgcgcg         >  1:224841/1‑150 (MQ=255)
                                                                                                                                         gTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCgca                                           >  1:2159222/1‑112 (MQ=255)
                                                                                                                                         gTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCgca                                           <  2:2159222/112‑1 (MQ=255)
                                                                                                                                            cACCCGCGAACGAGTATGAAGCCGTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCGCAGCACGCCGCCTCCTATTACGCCGCCAGCAGCGCGCCTCAcc  >  2:1185461/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                     
GGGCCGCAGCAGTCGTGGCAGCCGGGCTCGCACTCGAAAGAGGGGATGAGTTCGCGCAGGAAGTGGATCTTGTGGCGGTTGCAGGACATGACGGGATACAAGCTGACAAAGTGATTGGATTATAGGCTTATTCGCGGGTTCACCCGCGAACGAGTATGAAGCCCTCGCTGGTTGCCCGGCCGCGATGGCACAGCTTATCCTCGTCCCCCCTCCATTTTTCAGGGCCGCCCCAATGATCCACAGCGCGCAGCACGCCGCCTCCTATTACGCCGCCAGCAGCGCGCCACACC  >  NC_002947/6021617‑6021906

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: