Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,301,859 A→G 5.8% intergenic (+33/‑65) gloA → / → PP_3767 lactoylglutathione lyase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,301,8590AG 5.8% 1060.4 / 34.4 341intergenic (+33/‑65)gloA/PP_3767lactoylglutathione lyase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (149/172);  new base G (12/8);  total (161/180)
Fisher's exact test for biased strand distribution p-value = 2.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.58e-01

Only 100 of 345 total aligned reads displayed.

CTGGAAGTGCCGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGCCCTATGCC  >  NC_002947/4301710‑4302007
                                                                                                                                                     |                                                                                                                                                    
  ggAAGTGCCGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACt                                                                                                                                                    <  1:6697705/150‑1 (MQ=255)
   gAAGTGCCGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTc                                                                                                                                                   <  2:7451939/150‑1 (MQ=255)
      gTGCCGTTCCAGAAACGCCTGCAGGACTGCCGCATGAACCTCCTTGCCTTCGTCAAGGACCCGGACGTTTATTGGGTCGAAGTGATCCAGCCTACCTAGCTCGAAGGCTAAGCGTTAGTCGCGGCGCGCTCGTCTGGCCGGTAAATCCCg                                                                                                                                                <  2:607048/150‑1 (MQ=255)
         ccGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAGTTCCCGGcc                                                                                                                                             <  2:8818545/150‑1 (MQ=255)
         ccGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAGTTCCCGGcc                                                                                                                                             <  1:4216830/150‑1 (MQ=255)
          cGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAGTTCCCGGCct                                                                                                                                            <  1:7439681/150‑1 (MQ=255)
          cGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGAGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAGTTCCCGGCct                                                                                                                                            <  1:3712828/150‑1 (MQ=255)
          cGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCctc                                                                                                                                           >  1:2617033/1‑150 (MQ=255)
              ccAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTg                                                                                                                                        <  2:7757230/150‑1 (MQ=255)
              ccAGAAACGCCTGCAGGACGGCCGCATGAACCACCCGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTg                                                                                                                                        <  1:7551875/150‑1 (MQ=255)
                aGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCt                                                                                                                                      >  2:6853702/1‑150 (MQ=255)
                  aaaCGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCtgt                                                                                                                                    <  1:1344677/150‑1 (MQ=255)
                  aaaCGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCtgt                                                                                                                                    >  2:6536859/1‑150 (MQ=255)
                   aaCGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCtgtg                                                                                                                                   <  1:8626859/150‑1 (MQ=255)
                   aaCGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCtgtg                                                                                                                                   <  2:3966489/150‑1 (MQ=255)
                      gCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGACGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGta                                                                                                                                <  1:5946618/150‑1 (MQ=255)
                      gCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGta                                                                                                                                >  1:3056885/1‑150 (MQ=255)
                      gCCCGCAGGACGGCCGCATGAACCCCCTGGCCGTCGTCACGGACCCGGACGGATATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGta                                                                                                                                <  1:6418956/150‑1 (MQ=255)
                        cTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGtata                                                                                                                              >  1:6979068/1‑150 (MQ=255)
                        cTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGtata                                                                                                                              <  2:4325978/150‑1 (MQ=255)
                        cTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGtata                                                                                                                              <  1:4339912/150‑1 (MQ=255)
                         tGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTatat                                                                                                                             <  2:582260/150‑1 (MQ=255)
                         tGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTatat                                                                                                                             >  1:7838492/1‑150 (MQ=255)
                         tGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTatat                                                                                                                             <  2:8203153/150‑1 (MQ=255)
                          gCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATg                                                                                                                            >  1:722825/1‑150 (MQ=255)
                          gCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGGGATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCGTGCTGTGGTATATg                                                                                                                            >  1:2168178/1‑150 (MQ=255)
                           cAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGa                                                                                                                           <  1:4699019/150‑1 (MQ=255)
                           cAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGa                                                                                                                           <  2:5643891/150‑1 (MQ=255)
                            aGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGaa                                                                                                                          <  1:3519188/150‑1 (MQ=255)
                             ggACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAg                                                                                                                         <  1:2354596/150‑1 (MQ=255)
                                   ccGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACt                                                                                                                   >  2:2221786/1‑150 (MQ=255)
                                   ccGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACt                                                                                                                   <  2:5173798/150‑1 (MQ=255)
                                      cATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCt                                                                                                                >  2:4453250/1‑150 (MQ=255)
                                       aTGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCtt                                                                                                               <  2:5075880/150‑1 (MQ=255)
                                       aTGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCtt                                                                                                               <  2:3305179/150‑1 (MQ=255)
                                           accaccTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTcaca                                                                                                           >  2:1338747/1‑150 (MQ=255)
                                           accaccTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTcaca                                                                                                           >  1:8799327/1‑150 (MQ=255)
                                              accTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATAc                                                                                                        <  1:2678921/150‑1 (MQ=255)
                                              accTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATAc                                                                                                        >  1:1881268/1‑150 (MQ=255)
                                              accTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATAc                                                                                                        <  2:7638035/150‑1 (MQ=255)
                                                  ggCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACcaca                                                                                                    <  2:8132635/150‑1 (MQ=255)
                                                    ccTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACCTACCACAGc                                                                                                  <  2:2991191/150‑1 (MQ=255)
                                                        cGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAgg                                                                                              <  2:2444070/150‑1 (MQ=255)
                                                              ggACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGa                                                                                        <  2:1914546/150‑1 (MQ=255)
                                                               gACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGtata                                                                                                                              >  2:325487/1‑111 (MQ=255)
                                                               gACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGtata                                                                                                                              <  1:325487/111‑1 (MQ=255)
                                                               gACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGAc                                                                                       <  1:1683710/150‑1 (MQ=255)
                                                                    ggACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGc                                                                                  >  2:7197027/1‑150 (MQ=255)
                                                                        ggTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAAtc                                                                              <  1:6903121/150‑1 (MQ=255)
                                                                        ggTTATTGGGTAGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCAAATACCACAGCGACGAAAGGACGATGCAAtc                                                                              >  2:4687268/1‑150 (MQ=255)
                                                                            aTTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCtt                                                                          <  1:5913384/150‑1 (MQ=255)
                                                                               gggTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCAc                                                                       <  2:7053584/150‑1 (MQ=255)
                                                                                  tCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACtgt                                                                    <  1:5402473/150‑1 (MQ=255)
                                                                                  tCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACtgt                                                                    <  1:7970123/150‑1 (MQ=255)
                                                                                  tCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACtgt                                                                    >  1:6231393/1‑150 (MQ=255)
                                                                                         gATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATAcc                                                             >  1:5918079/1‑150 (MQ=255)
                                                                                          aTCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTCAGGACGATGCGATCTCTTCACTGTGAATACcg                                                            >  2:2664752/1‑150 (MQ=255)
                                                                                            ccAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACcgcg                                                          >  2:7844218/1‑150 (MQ=255)
                                                                                            ccAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACcgcg                                                          >  2:876503/1‑150 (MQ=255)
                                                                                                ccAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGAGGCAATCTCTTCACTGTGAATCCCGCGAGca                                                      >  2:490297/1‑150 (MQ=255)
                                                                                                    ccGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGAGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACAcc                                                  <  2:858238/150‑1 (MQ=255)
                                                                                                     cGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCa                                                 >  1:144949/1‑150 (MQ=255)
                                                                                                     cGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCa                                                 <  1:1099316/150‑1 (MQ=255)
                                                                                                     cGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCa                                                 >  1:8365556/1‑150 (MQ=255)
                                                                                                     cGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCa                                                 >  1:3977070/1‑150 (MQ=255)
                                                                                                       aGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATc                                               <  2:2499838/150‑1 (MQ=255)
                                                                                                          tCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACt                                            <  1:3337705/150‑1 (MQ=255)
                                                                                                          tCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACt                                            >  2:6355048/1‑150 (MQ=255)
                                                                                                          tCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACt                                            >  2:1461265/1‑150 (MQ=255)
                                                                                                           cGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTg                                           >  2:3150307/1‑150 (MQ=255)
                                                                                                            gAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGc                                          >  1:4149573/1‑150 (MQ=255)
                                                                                                                 cTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCg                                     >  1:5188925/1‑150 (MQ=255)
                                                                                                                  tAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGt                                    >  2:339714/1‑150 (MQ=255)
                                                                                                                  tAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTACCGGCCTCTTGCTGTGGTATATTAAGGTAACTGCTTCACATACCAAAGAGAGGTAAGGAAGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGt                                    >  1:660550/1‑150 (MQ=255)
                                                                                                                    aGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGa                                  >  2:7249728/1‑150 (MQ=255)
                                                                                                                    aGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGa                                  >  1:8174193/1‑150 (MQ=255)
                                                                                                                        ttCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGcc                              <  2:4212129/150‑1 (MQ=255)
                                                                                                                        ttCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGcc                              <  2:5427668/150‑1 (MQ=255)
                                                                                                                           gtcgCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCa                           >  1:5711439/1‑150 (MQ=255)
                                                                                                                               cggcgCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCCCTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATccc                       <  1:1587473/150‑1 (MQ=255)
                                                                                                                                      cTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCAc                >  2:8445323/1‑150 (MQ=255)
                                                                                                                                         gTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCAcccc             >  1:250067/1‑150 (MQ=255)
                                                                                                                                          tCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCt            <  2:1523211/150‑1 (MQ=255)
                                                                                                                                          agaGGCCGGGAGTTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCt            >  2:1691584/3‑150 (MQ=255)
                                                                                                                                           gaGGCCGGGAGTTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTg           >  2:7306527/2‑150 (MQ=255)
                                                                                                                                            aGGCCGGGAGTTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGc          >  1:2687663/1‑150 (MQ=255)
                                                                                                                                            aGGCCGGGAGTTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGc          >  1:4914758/1‑150 (MQ=255)
                                                                                                                                            aGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGc          <  1:4570414/150‑1 (MQ=255)
                                                                                                                                            aGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGc          >  2:1213002/1‑150 (MQ=255)
                                                                                                                                             ggCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGcc         <  2:3635631/150‑1 (MQ=255)
                                                                                                                                             ggCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGcc         <  2:6463964/150‑1 (MQ=255)
                                                                                                                                              gCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGccc        >  2:1052415/1‑150 (MQ=255)
                                                                                                                                               ccGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGCCCt       >  2:4634095/1‑150 (MQ=255)
                                                                                                                                               ccGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGCCCt       <  1:4507941/150‑1 (MQ=255)
                                                                                                                                                cGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGGCGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGCCCTa      >  1:993468/1‑150 (MQ=255)
                                                                                                                                                cGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGCCCTa      <  2:4864375/150‑1 (MQ=255)
                                                                                                                                                cGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGCCCTa      >  1:7106410/1‑150 (MQ=255)
                                                                                                                                                cGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGCCCTa      >  2:2336720/1‑150 (MQ=255)
                                                                                                                                                  ggAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGCCCTATg    <  2:3057635/150‑1 (MQ=255)
                                                                                                                                                    aaCTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGCCCTATgcc  >  1:50843/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTGGAAGTGCCGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGCCCTATGCC  >  NC_002947/4301710‑4302007

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