Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,397,891 T→C 11.3% intergenic (+106/‑80) PP_1214 → / → ruvC transcriptional regulator/crossover junction endodeoxyribonuclease RuvC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,397,8910TC11.3% 367.6 / 29.4 136intergenic (+106/‑80)PP_1214/ruvCtranscriptional regulator/crossover junction endodeoxyribonuclease RuvC
Reads supporting (aligned to +/- strand):  ref base T (80/36);  new base C (10/5);  total (94/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

Only 100 of 142 total aligned reads displayed.

CAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGACCGCCC  >  NC_002947/1397743‑1398040
                                                                                                                                                      |                                                                                                                                                      
cAATGCCCAGATTTCCGACGAGATCATGGAATACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATgtggtg                                                                                                                                                      <  2:2927634/150‑1 (MQ=255)
cAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATgtggtg                                                                                                                                                      <  2:4594265/150‑1 (MQ=255)
cAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGGGAACCCGCGCCGACAGGGGTTACGCCCCATTAAAATgtggtg                                                                                                                                                      >  2:4713267/1‑150 (MQ=255)
  aTGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCa                                                                                                                                                    <  2:3673342/150‑1 (MQ=255)
  aTGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCa                                                                                                                                                    <  2:3673309/150‑1 (MQ=255)
   tGCCCATATTTCCGACGAGATCATGGAAAACCTCGGCTTAGGTGATTTGGCATTTTGTGCCAACTCTGAGAGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGTTTAAACCCGCTACTACAGGTATTTCGCCCATTTCAAATGTGGTGCAg                                                                                                                                                   <  1:765026/150‑1 (MQ=255)
     cccAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTa                                                                                                                                                 <  1:3673228/150‑1 (MQ=255)
     cccAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGGAAACCCGCGCCGACAGGGATGACCCCCCATTCAAAGGTGGTGccgga                                                                                                                                                 >  2:2502923/1‑146 (MQ=255)
      ccAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTTAACCCGCTCCAACAGGTATTACGCCCCATTCAAATGTGGTGCCGTAc                                                                                                                                                >  2:3463578/1‑150 (MQ=255)
       cAGATATCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAATACCGCTCCGAACGGTAATTAGCCCCATTCAAATGGGGTGCAGTAcc                                                                                                                                               >  2:4303954/1‑150 (MQ=255)
        aGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCt                                                                                                                                              <  2:5612935/150‑1 (MQ=255)
           tttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTa                                                                                                                                           <  2:3255665/150‑1 (MQ=21)
           tttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTa                                                                                                                                           <  2:1897019/150‑1 (MQ=21)
            ttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAg                                                                                                                                          <  1:345264/150‑1 (MQ=255)
            ttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAg                                                                                                                                          <  1:4592392/150‑1 (MQ=255)
                gacgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGc                                                                                                                                      <  1:5462700/150‑1 (MQ=255)
                gacgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGc                                                                                                                                      <  2:4138137/150‑1 (MQ=18)
                  cgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCgg                                                                                                                                    <  2:3394817/150‑1 (MQ=255)
                  cgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCgg                                                                                                                                    <  2:2782954/150‑1 (MQ=255)
                  cgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCgg                                                                                                                                    <  1:743076/150‑1 (MQ=255)
                  cgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCgg                                                                                                                                    <  2:5089330/150‑1 (MQ=255)
                        cATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTAc                                                                                                                              <  2:1742333/150‑1 (MQ=255)
                        cATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTAc                                                                                                                              <  2:880129/150‑1 (MQ=255)
                         aTGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTAcc                                                                                                                             <  1:4719153/150‑1 (MQ=255)
                         aTGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTAcc                                                                                                                             <  1:4719280/150‑1 (MQ=255)
                           ggAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACACGTATTAAGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGAGGGTTTACCcg                                                                                                                           >  2:4215675/1‑150 (MQ=255)
                           ggAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCTACAGGGATTACGCCCAGTTAAAATGTGGTGCAGTACCGGGAGGAGCGGGTTTAACcg                                                                                                                           >  2:4215517/1‑150 (MQ=38)
                            gAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCTACAGGTATTACGCCCAGTTAAAAAGTGGTGCAGTACCAGTAGGAGCGGGGTTAACcgc                                                                                                                          >  2:1273044/1‑150 (MQ=38)
                            gAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCGTTCGCGGGTGATGCCGCGCCAACCGGTATTAGGGCCCGTGGAAAGGTGGTGGGGGAGCGGTTGGAGGGGGGTTAGCcgc                                                                                                                          >  1:3216237/1‑150 (MQ=255)
                                 ccTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTAACTGTAGGAGCGGGTTTACCCGCGAAAg                                                                                                                     >  2:4210017/1‑150 (MQ=255)
                                                               cAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTTAACCCGCGCCAACACGGATTACCCCCAATTACAAAGTGGGGCAGTAACCGTAGGTGCCGGGTTAACCGCCCAAGGGCCCGCGCGGTTTTTAGGGATGcgaaaat                                                                                       >  2:287406/1‑144 (MQ=255)
                                                               cAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCAACACGTAATACCCCCATTTCAAATGTGGTGCAGTACCGGTAGGAGCGGGGTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGAAGCATCAtt                                                                                       >  2:1264932/1‑150 (MQ=37)
                                                                      gAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCGTTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGGATCATTTTTTCGa                                                                                >  2:2242375/1‑150 (MQ=255)
                                                                       aGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCCCACCACAAGGTATTACGCCCCTTTCAAATGTGGGGCAGTAACGGGAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCCGTTTGTAGAGATGCCTAATTTTTTCGAc                                                                               >  2:2948915/1‑150 (MQ=37)
                                                                         cGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCCCTCCTACAGGTATTACGCCCCAGTCAAATGTGGTGCACTAACGGTAGGAGCCGGGTTAACCGCGAAAGGGCCCGCGCACGTTGTAGAGATGCAACAAGTTTTCGACgg                                                                             >  2:3730950/1‑150 (MQ=25)
                                                                           ttttATCGCGCGTGTGCCGGCCCTTGCGGGTTTAACCCGGCCCCCACCGGTATACCGCCCTTTCCAATTGTGGTGCTTTACCTTTAGGCGGGGTTTAACCCGCGAATGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGc                                                                           <  2:4354331/150‑1 (MQ=255)
                                                                              gATCGCGCCTGTGCCGGCCCTTTCGCGGTTAAACCCGCCCCTACTGTTATTACGCCCTTTTCAAATGTGGTGCTGTACCTGTTGGAGGGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCttt                                                                        <  2:1846021/150‑1 (MQ=37)
                                                                                tcgcgcGTGGGCCGGCCCTTGGGCGGGTACCCCGGCTCCTCCTGTTTTCCCGCCAATTGCAATGGTGGTGCCGTCCGTGTTGGGGGGGTTTTACCCGGTAACGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTAt                                                                      <  2:2567521/148‑1 (MQ=255)
                                                                                tcgcgcCGGGGCCGGCCGTTGCGGGGGTAACCCGGCGCGTCCCGGTATTCCGCCCTTTTCAATGGTGGTGCTGTCCGGGTCGGGGGGGTTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTAt                                                                      <  2:5615009/148‑1 (MQ=255)
                                                                                    gcgtgggccgccCCTGTCGGGTCTAACCCGGCTCCACCCGGTATACCCCCCATTTCAATGGGGTGGCCTCACCTTTAGGGGGGGTTTCACCCGCGAATGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGAct                                                                  <  2:1606156/144‑1 (MQ=255)
                                                                                      cTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTAAAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCCGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGActct                                                                >  2:5545014/1‑150 (MQ=255)
                                                                                             gccccTTGCGCGGTTAACCCCGCGCCTACCGTTATACCCCACATTTCAATGGTGGTGCCGACCCGGTTGGAGCGGTTTACCCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTa                                                         <  2:5281642/148‑1 (MQ=255)
                                                                                              gCCCTTGCGGGGGTAAACCCGGCTCTAACGGTTTTACCCCCATTTAAATTTGGGGGGCCGGCCGTGTGGGAGGGTTTTCACCCGCGAACGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTa                                                         <  2:1860829/150‑1 (MQ=255)
                                                                                                ccTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCCTTTCAAATGTGGTGCAGTACCTGTGGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCCGTTTTTAGATATGCCTCATTTTTTCGACGGGGTTTATGACTCGGATTATAGGt                                                      >  2:756205/1‑150 (MQ=25)
                                                                                                         ggTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTGGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGCTTCTAGGTATCGAcccc                                             >  2:4773723/1‑150 (MQ=255)
                                                                                                          gTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCg                                            >  1:1693442/1‑150 (MQ=255)
                                                                                                           tAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCg                                                                                 >  1:2507321/1‑112 (MQ=255)
                                                                                                           tAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCg                                                                                 <  2:2507321/112‑1 (MQ=255)
                                                                                                           tAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCgg                                           >  2:175295/1‑150 (MQ=255)
                                                                                                              aCCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTc                                        >  2:4487617/1‑150 (MQ=255)
                                                                                                              aCCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTc                                        >  1:3183683/1‑150 (MQ=255)
                                                                                                                tcGCTCCTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc                                      <  1:4154527/149‑1 (MQ=18)
                                                                                                                ccGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc                                      >  1:3409097/1‑150 (MQ=255)
                                                                                                                ccGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc                                      >  1:3549983/1‑150 (MQ=255)
                                                                                                                ccGCTCCTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAACGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc                                      >  1:1882002/1‑150 (MQ=18)
                                                                                                                 cGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgcg                                     >  1:2636855/1‑150 (MQ=255)
                                                                                                                  gCTCCTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCGAGGTATCGACCCCGGCTCgcgc                                    >  2:2909495/1‑150 (MQ=18)
                                                                                                                    tCCTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGAAGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCAt                                  >  2:2563398/1‑150 (MQ=18)
                                                                                                                     ccTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATc                                 >  1:258441/1‑150 (MQ=255)
                                                                                                                      cTACAGTTATTCCGCCAATTTAAATTGTGGTGCAGCCCCTTTTGGAGGGGTTTAACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCa                                <  2:2165777/150‑1 (MQ=255)
                                                                                                                      cTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCa                                >  1:2964950/1‑150 (MQ=255)
                                                                                                                      cTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCa                                >  2:879690/1‑150 (MQ=255)
                                                                                                                      cTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCa                                >  2:4147370/1‑150 (MQ=255)
                                                                                                                       tACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAc                               >  1:5450429/1‑150 (MQ=21)
                                                                                                                       tACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAc                               >  2:2658925/1‑150 (MQ=21)
                                                                                                                       tACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAc                               >  1:3965967/1‑150 (MQ=21)
                                                                                                                       tACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAc                               >  1:2292607/1‑150 (MQ=21)
                                                                                                                       tACAGGTACTGCACCACATTTGAAATGGGAGTAATACCTGTAGGAGCGGGTTTACCCGCGAATGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAc                               >  1:3452668/1‑150 (MQ=21)
                                                                                                                        aCAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAcc                              >  2:2199441/1‑150 (MQ=255)
                                                                                                                        aCAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAcc                              >  1:1049564/1‑150 (MQ=255)
                                                                                                                         cAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATAACCg                             >  1:4950502/1‑150 (MQ=255)
                                                                                                                          aGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCgg                            >  1:2450205/1‑150 (MQ=255)
                                                                                                                          aGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCgg                            >  1:4515405/1‑150 (MQ=255)
                                                                                                                          aGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCgg                            >  2:613422/1‑150 (MQ=255)
                                                                                                                          aGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCgg                            >  1:4515657/1‑150 (MQ=255)
                                                                                                                           ggTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGc                           >  1:280952/1‑150 (MQ=255)
                                                                                                                            gTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCt                          >  1:4950151/1‑150 (MQ=255)
                                                                                                                            gTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCt                          >  2:2399730/1‑150 (MQ=255)
                                                                                                                               ttACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACGGGCTATg                       >  1:46050/1‑150 (MQ=255)
                                                                                                                                  cGCCCATTTCAATTGTGGTGCAGTCCCTGTCGGTGCGGGTTCACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCg                    <  2:4497820/150‑1 (MQ=255)
                                                                                                                                   gCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTCGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt                   >  1:4073859/1‑150 (MQ=255)
                                                                                                                                   gCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt                   >  2:4950455/1‑150 (MQ=255)
                                                                                                                                   gCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt                   >  1:1434410/1‑150 (MQ=255)
                                                                                                                                   gCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt                   >  1:215867/1‑150 (MQ=255)
                                                                                                                                   gCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt                   >  2:2578665/1‑150 (MQ=255)
                                                                                                                                   gCCCATTTCAAATGTGGTGCAGTACCGGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt                   >  2:2578688/1‑150 (MQ=255)
                                                                                                                                   accacATTTGAAATGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGt                   >  1:576119/5‑150 (MQ=37)
                                                                                                                                    cccATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTg                  <  1:2025211/150‑1 (MQ=255)
                                                                                                                                    cccATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTg                  >  1:4702389/1‑150 (MQ=255)
                                                                                                                                    cccATATCAAAAGGCGAGCAGTACCTGGAGGAGCCGGTTTACCCGCGAAAGGGCCGGCGCAGAGTGTAGAGATGCATAATATTTTAGACGCGCATTCTGACTCAGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTg                  >  2:2561693/1‑150 (MQ=255)
                                                                                                                                     ccATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTgg                 >  1:2024596/1‑150 (MQ=255)
                                                                                                                                     ccATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTgg                 >  2:4613762/1‑150 (MQ=255)
                                                                                                                                     ccATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTgg                 >  1:663356/1‑150 (MQ=255)
                                                                                                                                         tttCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTAc              >  1:1935811/1‑150 (MQ=255)
                                                                                                                                          ttCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACg             >  2:1839605/1‑150 (MQ=255)
                                                                                                                                          ttCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACg             >  2:1027137/1‑150 (MQ=255)
                                                                                                                                             aaaTGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCa          >  1:5344716/1‑150 (MQ=255)
                                                                                                                                             aaaTGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCa          >  1:2886836/1‑150 (MQ=255)
                                                                                                                                                     gtgCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGACCGcc   >  2:2059345/1‑150 (MQ=255)
                                                                                                                                                      tgCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGACCGccc  >  1:2031124/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
CAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGACCGCCC  >  NC_002947/1397743‑1398040

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