Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,296,910 A→G 6.2% 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,296,9100AG 6.2% 231.1 / 3.6 82intergenic (+33/‑65)gloA/PP_3767lactoylglutathione lyase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/50);  new base G (2/3);  total (28/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

CTGGAAGTGCCGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGCCCTAT  >  NC_002947/4296761‑4297055
                                                                                                                                                     |                                                                                                                                                 
cTGGAAGTGCCGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGaa                                                                                                                                                   >  1:1028455/1‑150 (MQ=255)
cTGGAAGTGCCGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGaa                                                                                                                                                   >  2:1805110/1‑150 (MQ=255)
cTGGAAGTGCCGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGaa                                                                                                                                                   <  1:12809/150‑1 (MQ=255)
 tGGAAGTGCCGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAAc                                                                                                                                                  <  2:1803640/150‑1 (MQ=255)
      gTGCCGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCg                                                                                                                                             <  1:1885208/150‑1 (MQ=255)
       tGCCGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCgg                                                                                                                                            <  2:386326/150‑1 (MQ=255)
        gCCGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAGTTCCCGGc                                                                                                                                           >  2:107554/1‑150 (MQ=255)
         ccGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCAATTTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAGTTCCCGGcc                                                                                                                                          <  2:732027/150‑1 (MQ=255)
            ttCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAGTTCCCGGCctga                                                                                                                                       <  1:107554/150‑3 (MQ=255)
            ttCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCctct                                                                                                                                       <  2:1381424/150‑1 (MQ=255)
             tCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCtt                                                                                                                                      <  1:724300/150‑1 (MQ=255)
                aGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGCGTCCTCCCGCCCCCTGGca                                                                                                                                   >  2:871887/1‑149 (MQ=255)
                 gagaggCCTGCAGGACGGCGGCATGAACCACCTGGGCTGCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCtg                                                                                                                                  <  2:775159/145‑1 (MQ=255)
                   aaCGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCtgtg                                                                                                                                <  1:75496/150‑1 (MQ=255)
                    aCGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTgg                                                                                                                               <  2:45353/150‑1 (MQ=255)
                     cGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGt                                                                                                                              <  1:721221/150‑1 (MQ=255)
                     cGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGt                                                                                                                              <  1:1805110/150‑1 (MQ=255)
                        cTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGtata                                                                                                                           <  2:411319/150‑1 (MQ=255)
                            aGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGGCGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGaa                                                                                                                       <  2:1245443/150‑1 (MQ=255)
                            aGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGaa                                                                                                                       <  2:1695971/150‑1 (MQ=255)
                             ggACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAg                                                                                                                      <  2:1688406/150‑1 (MQ=255)
                             ggACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAg                                                                                                                      <  2:797529/150‑1 (MQ=255)
                              gACGGCCGCATGAACCACCTGGCCTTCGCCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAgg                                                                                                                     >  2:1766482/1‑150 (MQ=255)
                                      cATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCt                                                                                                             >  2:1629909/1‑150 (MQ=255)
                                           accaccTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTcaca                                                                                                        <  1:1294421/150‑1 (MQ=255)
                                           accaccTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTcaca                                                                                                        >  1:1657590/1‑150 (MQ=255)
                                             caccTGGCCTGCGTCAAGGGCCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCTCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATa                                                                                                      <  2:414259/150‑1 (MQ=255)
                                               ccTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATAcc                                                                                                    >  2:671599/1‑150 (MQ=255)
                                                 tGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGACTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACcac                                                                                                  >  2:306217/1‑150 (MQ=255)
                                                  ggCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACcaca                                                                                                 <  1:325759/150‑1 (MQ=255)
                                                     cTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGTTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCg                                                                                              >  1:595364/1‑150 (MQ=255)
                                                     cTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTatat                                                                                                                          >  2:360131/1‑122 (MQ=255)
                                                     cTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTatat                                                                                                                          <  1:360131/122‑1 (MQ=255)
                                                     cTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTaa                                                                                                                  <  1:481157/130‑1 (MQ=255)
                                                     cTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTaa                                                                                                                  >  2:481157/1‑130 (MQ=255)
                                                      ttCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGa                                                                                             <  1:580326/150‑1 (MQ=255)
                                                           cAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGGTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTaa                                                                                        <  1:1561870/150‑1 (MQ=255)
                                                           cAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTaa                                                                                        <  1:1600361/150‑1 (MQ=255)
                                                            aaGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCt                                                                                                                                   >  1:1584981/1‑106 (MQ=255)
                                                            aaGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCt                                                                                                                                   <  2:1584981/106‑1 (MQ=255)
                                                            aaGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAg                                                                                       <  1:1766482/150‑1 (MQ=255)
                                                               gACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGa                                                                                     <  2:278509/149‑1 (MQ=255)
                                                               gACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGa                                                                                     >  1:278509/1‑149 (MQ=255)
                                                                aCCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACg                                                                                   <  1:856444/150‑1 (MQ=255)
                                                                  ccGGACGGTTATTGGGTCGACGTGATGCAACCCCCCGAGCTCCAAGGCTAAGCGTTCGTCGCGGCGCGCGCGGCGGGCCGGGAACTCCCGGCCTCCTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGAt                                                                                 <  1:124105/150‑1 (MQ=255)
                                                                  ccGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGAt                                                                                 <  1:466745/150‑1 (MQ=255)
                                                                   cGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGt                                                                                                                              <  2:458333/104‑1 (MQ=255)
                                                                   cGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGt                                                                                                                              >  1:458333/1‑104 (MQ=255)
                                                                   cGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATg                                                                                >  1:255823/1‑150 (MQ=255)
                                                                    ggACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGc                                                                               <  2:1390234/150‑1 (MQ=255)
                                                                      aCGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCaa                                                                             <  2:298596/150‑1 (MQ=255)
                                                                      aCGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCaa                                                                             <  2:656448/150‑1 (MQ=255)
                                                                      aCGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCaa                                                                             >  2:1154444/1‑150 (MQ=255)
                                                                           tATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATctct                                                                        <  2:1292015/150‑1 (MQ=255)
                                                                            aTTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTaa                                                                                        <  1:237610/133‑1 (MQ=255)
                                                                            aTTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTaa                                                                                        >  2:237610/1‑133 (MQ=255)
                                                                             ttGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAGTTCCCGGCctga                                                                                                                                       >  1:1723001/1‑83 (MQ=255)
                                                                             ttGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAGTTCCCGGCctga                                                                                                                                       <  2:1723001/85‑3 (MQ=255)
                                                                             ttGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTc                                                                      >  1:1613987/1‑150 (MQ=255)
                                                                               gggTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCtt                                                                       <  1:256616/147‑1 (MQ=255)
                                                                               gggTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCtt                                                                       >  2:256616/1‑147 (MQ=255)
                                                                                 gTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACtg                                                                  >  2:1521281/1‑150 (MQ=255)
                                                                                  tCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACtgt                                                                 <  2:595364/150‑1 (MQ=255)
                                                                                           tCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACcgc                                                        <  2:766059/150‑1 (MQ=255)
                                                                                             cAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGa                                                      <  1:1092059/150‑1 (MQ=255)
                                                                                                         cTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCGCTTCACTGTGAATACCGCGAGCACACCATCAc                                          >  1:669742/1‑150 (MQ=255)
                                                                                                          tCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGTGATCTTCCGGCCTCTTGCTGTGGTATATGAAGGGAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACt                                         <  1:1560891/150‑1 (MQ=255)
                                                                                                          tCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACt                                         >  2:1465313/1‑150 (MQ=255)
                                                                                                                gCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAAc                                                                                                                 >  1:1180437/1‑72 (MQ=255)
                                                                                                                gCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAAc                                                                                                                 <  2:1180437/72‑1 (MQ=255)
                                                                                                                    aGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGa                               <  2:1028455/150‑1 (MQ=255)
                                                                                                                    aGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGa                               <  1:1712402/150‑1 (MQ=255)
                                                                                                                      cGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATg                             <  2:255823/150‑1 (MQ=255)
                                                                                                                            tcgCCGCGCGCTCGCCAGCGCGTGAACTCCCTGCCTCTTTCTGTGGTATAGGAAGTTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATgaacgtc                       <  1:1499748/150‑7 (MQ=255)
                                                                                                                                       tCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTAACTAAGCATATTGCTAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCAcc            >  2:529050/1‑150 (MQ=255)
                                                                                                                                        cGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCAccc           <  1:1693568/150‑1 (MQ=255)
                                                                                                                                         gTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGCGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCAccc           <  2:1036614/149‑1 (MQ=255)
                                                                                                                                         gTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCAcccc          >  2:1777292/1‑150 (MQ=255)
                                                                                                                                         gTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCAccc           >  1:1036614/1‑149 (MQ=255)
                                                                                                                                              gCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGAc                 <  2:1821819/138‑1 (MQ=255)
                                                                                                                                              gCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGAc                 >  1:1821819/1‑138 (MQ=255)
                                                                                                                                               ccGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGt                                                                                          >  1:910988/1‑64 (MQ=255)
                                                                                                                                               ccGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGt                                                                                          <  2:910988/64‑1 (MQ=255)
                                                                                                                                                 gggAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGCCCTAt  <  1:671599/150‑1 (MQ=255)
                                                                                                                                                 gggAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGCCCTAt  <  2:851084/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CTGGAAGTGCCGTTCCAGAAACGCCTGCAGGACGGCCGCATGAACCACCTGGCCTTCGTCAAGGACCCGGACGGTTATTGGGTCGAAGTGATCCAGCCAACCGAGCTCGAAGGCTAAGCGTTCGTCGCGGCGCGCTCGTCAGGCCGGGAACTCCCGGCCTCTTGCTGTGGTATATGAAGGTAACTGCTTCACATACCACAGCGAGGTAAGGACGATGCAATCTCTTCACTGTGAATACCGCGAGCACACCATCACTGCCAGCGTGATGCCCCATCCCGACTCACCCCTGCCCTAT  >  NC_002947/4296761‑4297055

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: