Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,706 G→C 5.4% L39V (CTC→GTC)  parA ← chromosome partition protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,7060GC 5.4% 271.8 / 4.1 85L39V (CTC→GTC) parAchromosome partition protein
Reads supporting (aligned to +/- strand):  ref base G (42/38);  new base C (2/3);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGGCCGGCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGCCTTGCGGCGC  >  NC_002947/1557‑1854
                                                                                                                                                     |                                                                                                                                                    
gTTGGCCGGCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAg                                                                                                                                                      <  2:1209285/150‑1 (MQ=255)
        gcagcaACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATca                                                                                                                                              >  2:942725/1‑150 (MQ=255)
         cagcaACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATcag                                                                                                                                             <  1:278726/150‑1 (MQ=255)
           gcaACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCTGATCGACCTCGATccgca                                                                                                                                           <  2:334336/150‑1 (MQ=255)
            caACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCTGATCGACCTCGATccgcag                                                                                                                                          <  1:415019/150‑1 (MQ=255)
                  gAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTAGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACacg                                                                                                                                    >  1:540922/1‑150 (MQ=255)
                       cTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCt                                                                                                                               <  1:625137/150‑1 (MQ=255)
                        tCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCtt                                                                                                                              >  2:674979/1‑150 (MQ=255)
                        tCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCtt                                                                                                                              >  1:769514/1‑150 (MQ=255)
                                     tagtGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGcagccagc                                                                                                                 >  2:720160/1‑150 (MQ=255)
                                     tagtGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGAAGGTAGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGcagccagc                                                                                                                 >  1:740750/1‑150 (MQ=255)
                                       gtGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGa                                                                                                               >  2:1003912/1‑150 (MQ=255)
                                            tAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAggcgg                                                                                                          <  1:674979/150‑1 (MQ=255)
                                             agccgcggCCAAGTCGCCCTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCCCTGCCCATGGTGGCGTTCCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAggcggc                                                                                                         <  1:228945/144‑1 (MQ=255)
                                               ccTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCga                                                                                                       <  1:1111331/150‑1 (MQ=255)
                                                cTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCgag                                                                                                      >  1:123889/1‑150 (MQ=255)
                                                cTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCgag                                                                                                      >  1:494680/1‑150 (MQ=255)
                                                     ccAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGa                                                                                                 >  1:648400/1‑150 (MQ=255)
                                                       aaGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATa                                                                                               <  2:494680/150‑1 (MQ=255)
                                                         gTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACa                                                                                             >  1:1197017/1‑150 (MQ=255)
                                                            gCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAggt                                                                                          <  1:838415/150‑1 (MQ=255)
                                                            gCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAggt                                                                                          <  1:331525/150‑1 (MQ=255)
                                                                  cccGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGt                                                                                    >  1:1354089/1‑150 (MQ=255)
                                                                   ccGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATAGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTc                                                                                   <  1:460602/150‑1 (MQ=255)
                                                                    cGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGGCGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCt                                                                                  >  1:1248894/1‑150 (MQ=255)
                                                                    cGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTgg                                                                                                                            >  1:974791/1‑108 (MQ=255)
                                                                    cGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTgg                                                                                                                            <  2:974791/108‑1 (MQ=255)
                                                                    cGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCt                                                                                  <  2:655310/150‑1 (MQ=255)
                                                                     gATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCtt                                                                                 <  1:1255030/150‑1 (MQ=255)
                                                                       tGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGc                                                                               <  2:772886/150‑1 (MQ=255)
                                                                        gAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGcc                                                                              <  2:1197017/150‑1 (MQ=255)
                                                                        gAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGcc                                                                              <  2:740750/150‑1 (MQ=255)
                                                                          gCAGGTCGTAGAGCGGATGCTCGAGCTCGTGCTTGTCCACAGCGCTGCCCATGGTGGCGTTGCCCTGCGGAGCGAGGTCGATCAGCAGGGCACGACGGTTGGTCGCTGCGGGCGGGGCGGCGGGATTGATGCAGGTGGTtggctggcgcg                                                                            >  2:544431/1‑147 (MQ=255)
                                                                          gCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTgg                                                                                                                            <  2:1299254/102‑1 (MQ=255)
                                                                          gCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTgg                                                                                                                            >  1:1299254/1‑102 (MQ=255)
                                                                          gCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTg                                                                                                  <  1:1012865/128‑1 (MQ=255)
                                                                          gCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTg                                                                                                  >  2:1012865/1‑128 (MQ=255)
                                                                             ggTAGTAGACCGAATGCTCGAGCTCGCGCTTGTCCACAGGGCTGCCCATGGGGTCGTGGCCCTGCGGATCGAGGTCGAGCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCacac                                                                         >  2:235247/1‑150 (MQ=255)
                                                                              gTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTg                                                                                >  1:1057202/1‑142 (MQ=255)
                                                                              gTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTg                                                                                <  2:1057202/142‑1 (MQ=255)
                                                                                cGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCAcaccac                                                                      <  2:123889/150‑1 (MQ=255)
                                                                                 gTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACaccacc                                                                     >  1:389646/1‑150 (MQ=255)
                                                                                  tAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTTCCAACACCACCt                                                                    >  1:1002444/1‑150 (MQ=255)
                                                                                     aCCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTc                                                                 <  2:1354089/150‑1 (MQ=255)
                                                                                      ccGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGAt                                                                                                    >  2:713893/1‑114 (MQ=255)
                                                                                      ccGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGAt                                                                                                    <  1:713893/114‑1 (MQ=255)
                                                                                       cGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTg                                                               >  1:130443/1‑150 (MQ=255)
                                                                                       cGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTg                                                               <  1:462623/150‑1 (MQ=255)
                                                                                       cGAATGCTCGAGCTCGTGCTAGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTg                                                               <  2:540922/150‑1 (MQ=255)
                                                                                             cTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGTTCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTcgc                                                         <  1:1165224/150‑1 (MQ=255)
                                                                                             cTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTcgc                                                         >  1:779742/1‑150 (MQ=255)
                                                                                              tCGAACTAGTGCTTGTCCTTACCGGTGCTCATCATTGCCGGTCCCTTCGGATCGAGGTCGATCAGCATGACACGACGCTTGGTCGCAGCCAGCGAGGCAGCGAGATTGAGACAAGTGGTTGTCTTGCCCACACCACCTTTCTGGTTcgcg                                                        >  2:1020498/1‑150 (MQ=255)
                                                                                               cGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGa                                                       <  2:1329632/150‑1 (MQ=255)
                                                                                                gAGCTCGTGCTTGTCCACACCGCTGCCTATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGAt                                                      >  1:2514/1‑150 (MQ=255)
                                                                                                gAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAggtgg                                                                                        >  1:1421217/1‑116 (MQ=255)
                                                                                                gAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAggtgg                                                                                        <  2:1421217/116‑1 (MQ=255)
                                                                                                             tCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTa                                         >  2:1218948/1‑150 (MQ=255)
                                                                                                               cacaCCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCttt                                                                                                                             <  1:508535/64‑2 (MQ=255)
                                                                                                               cacaCCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCttt                                                                                                                             >  2:508535/1‑63 (MQ=255)
                                                                                                               cacaCCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGcagc                                                                                                                     <  2:1224062/72‑1 (MQ=255)
                                                                                                               cacaCCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGcagc                                                                                                                     >  1:1224062/1‑72 (MQ=255)
                                                                                                               cacaCCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGc                                       >  2:869561/1‑150 (MQ=255)
                                                                                                                acacCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTg                                                               >  2:1326414/1‑125 (MQ=255)
                                                                                                                acacCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTg                                                               <  1:1326414/125‑1 (MQ=255)
                                                                                                                acacCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGcc                                      >  1:340223/1‑150 (MQ=255)
                                                                                                                acacCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGcc                                      >  2:527885/1‑150 (MQ=255)
                                                                                                                  acCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATAcc                                            <  1:477148/142‑1 (MQ=255)
                                                                                                                  acCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATAcc                                            >  2:477148/1‑142 (MQ=255)
                                                                                                                    cGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCAttg                                  <  2:779742/150‑1 (MQ=255)
                                                                                                                     gCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGAtt                                                     >  1:1201487/1‑130 (MQ=255)
                                                                                                                     gCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGAtt                                                     <  2:1201487/130‑1 (MQ=255)
                                                                                                                      cTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCAttgtt                                <  1:1003912/150‑1 (MQ=255)
                                                                                                                       tGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCAttgttg                               >  2:341895/1‑150 (MQ=255)
                                                                                                                          ccATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGa                            <  2:648400/150‑1 (MQ=255)
                                                                                                                            aTGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAgt                          <  1:720160/150‑1 (MQ=255)
                                                                                                                               gtggCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCtt                                                                                 <  2:262964/92‑1 (MQ=255)
                                                                                                                               gtggCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCtt                                                                                 >  1:262964/1‑92 (MQ=255)
                                                                                                                                 ggCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTcc                     <  1:942725/150‑1 (MQ=255)
                                                                                                                                    gTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAg                  >  2:1101249/1‑150 (MQ=255)
                                                                                                                                         ccTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATg             >  2:323015/1‑150 (MQ=255)
                                                                                                                                           tGCTGATCGACCTCGATCCGCAGGGCACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGcc           >  2:533484/1‑150 (MQ=255)
                                                                                                                                           tGCTGATCGACCTCGATCCGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAg                                                                                            >  1:1409870/1‑69 (MQ=25)
                                                                                                                                           tGCTGATCGACCTCGATCCGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAg                                                                                            <  2:1409870/69‑1 (MQ=25)
                                                                                                                                                 tCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGcc           >  1:17936/1‑144 (MQ=255)
                                                                                                                                                 tCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGcc           <  2:17936/144‑1 (MQ=255)
                                                                                                                                                    aGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGCCTTGCGgcgc  >  2:956225/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTTGGCCGGCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGCCTTGCGGCGC  >  NC_002947/1557‑1854

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