Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,338,458 T→A 8.8% intergenic (+34/+108) PP_3812 → / ← PP_5632 lipase/allo‑isoleucine metabolism protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,338,4580TA 8.8% 249.5 / 12.3 88intergenic (+34/+108)PP_3812/PP_5632lipase/allo‑isoleucine metabolism protein
Reads supporting (aligned to +/- strand):  ref base T (44/36);  new base A (4/4);  total (48/40)
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.60e-01

GCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCG  >  NC_002947/4338310‑4338606
                                                                                                                                                    |                                                                                                                                                    
gcgcTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCttt                                                                                                                                                     <  1:298557/150‑1 (MQ=255)
gcgcTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCttt                                                                                                                                                     <  1:4293/150‑1 (MQ=255)
 cgcTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCtttt                                                                                                                                                    <  2:1905789/150‑1 (MQ=255)
      cAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTc                                                                                                                                               <  1:500778/150‑1 (MQ=255)
      cAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAAAAGACCCGCCCAGGCGGGTCTATTGCtca                                                                                                                                               <  2:827904/150‑3 (MQ=255)
       aaCGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCt                                                                                                                                              >  2:779328/1‑150 (MQ=255)
        aCGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTACGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTa                                                                                                                                             >  1:1098165/1‑150 (MQ=255)
            aaTGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTg                                                                                                                                         >  2:1960886/1‑150 (MQ=255)
             aTGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGc                                                                                                                                        >  1:282748/1‑150 (MQ=255)
             aTGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGc                                                                                                                                        <  2:282748/150‑1 (MQ=255)
                  tCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGaaa                                                                                                                                   <  2:656715/150‑1 (MQ=255)
                  tCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGaaa                                                                                                                                   <  1:1800905/150‑1 (MQ=255)
                    ctGAGTATGGGCTGCTGAACCCGTTGAGCCCAGTGCCGGCAGTAAAGGCGGCGCTGCGCCAGACGGGCAACGAACTTAACGTGCATCTGCCGTCAGGCGATGAGCAATCGACCTGCCTGGGCGGGTCTTTTGCTTCTATTTGCGAAAACa                                                                                                                                 >  2:275887/2‑150 (MQ=255)
                    aTGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACa                                                                                                                                 >  2:900657/1‑150 (MQ=255)
                        cTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGccc                                                                                                                             <  1:118979/150‑1 (MQ=255)
                            acttggctCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGc                                                                                                                                        <  2:961537/130‑1 (MQ=255)
                            acttggctCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGc                                                                                                                                        >  1:961537/6‑135 (MQ=255)
                             ggCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAgg                                                                                                                        >  2:649153/1‑150 (MQ=255)
                               ctgctCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCt                                                                                                                      >  1:2318143/1‑150 (MQ=255)
                                  ctCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCa                                                                                                                   <  2:1766803/150‑1 (MQ=255)
                                   tCAACCCGTTGGGCCAAGTGCAGGACCTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCaa                                                                                                                  <  2:879676/150‑1 (MQ=255)
                                   tCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAAGGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCtcaac                                                                                                                  >  2:2093816/1‑146 (MQ=255)
                                   tCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCaa                                                                                                                  >  1:1840191/1‑150 (MQ=255)
                                   tCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCaa                                                                                                                  >  1:1928586/1‑150 (MQ=255)
                                   tCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCaa                                                                                                                  >  1:366911/1‑150 (MQ=255)
                                    cAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAAt                                                                                                                 <  1:1631846/150‑1 (MQ=255)
                                        cggTTGAGACAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAAt                                                                                                             >  1:1927187/3‑150 (MQ=255)
                                        ccGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAAt                                                                                                             <  1:779328/150‑1 (MQ=255)
                                         cGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATc                                                                                                            >  1:1149183/1‑150 (MQ=255)
                                           ttGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTg                                                                                                          >  2:30593/1‑150 (MQ=255)
                                            tGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGc                                                                                                         <  1:1960886/150‑1 (MQ=255)
                                                  aaGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGc                                                                                                         <  2:892453/144‑1 (MQ=255)
                                                  aaGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGc                                                                                                         >  1:892453/1‑144 (MQ=255)
                                                   aGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGt                                                                                                  >  2:2318658/1‑150 (MQ=255)
                                                    gTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGTCAGCGGTg                                                                                                 >  1:1859669/1‑150 (MQ=255)
                                                          gCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTg                                                                                                                                         >  1:1237128/1‑104 (MQ=255)
                                                          gCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTg                                                                                                                                         <  2:1237128/104‑1 (MQ=255)
                                                          gCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAg                                                                                           >  1:899845/1‑150 (MQ=255)
                                                              tAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGTCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGcc                                                                                                                              <  2:572893/111‑1 (MQ=255)
                                                              tAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGcc                                                                                                                              >  1:572893/1‑111 (MQ=255)
                                                                aaGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCgg                                                                                     <  1:30593/150‑1 (MQ=255)
                                                                aaGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCgg                                                                                     <  2:1840191/150‑1 (MQ=255)
                                                                 aGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGc                                                                                    <  2:420715/150‑1 (MQ=255)
                                                                       cgcTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTc                                                                              <  2:561513/150‑1 (MQ=255)
                                                                          tGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTcgcg                                                                           >  2:791951/1‑150 (MQ=255)
                                                                           gCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCgg                                                                          <  1:1478765/150‑1 (MQ=255)
                                                                           acGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCgg                                                                          <  2:1098165/149‑1 (MQ=255)
                                                                             ggCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTg                                                                                                                                                   <  1:2338036/75‑1 (MQ=255)
                                                                             ggCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTg                                                                                                                                                   >  2:2338036/1‑75 (MQ=255)
                                                                                     ggCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCcgc                                                                <  1:203592/150‑1 (MQ=255)
                                                                                     ggCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCcgc                                                                <  2:899845/150‑1 (MQ=255)
                                                                                        aaCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTGTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGGAAACCCGCGaa                                                             >  2:732615/1‑150 (MQ=255)
                                                                                         ccGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGCCCCGCCTGGGCGGGTCTTTTCCTTCTATTGTCCGCAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGGCCGGCCCTTTCGCGGGTTTACCCGCGaaa                                                            <  1:275887/149‑1 (MQ=255)
                                                                                         aCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCTaaa                                                            <  2:42807/150‑1 (MQ=255)
                                                                                         aCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAACCCAGCGCATTGGCCGGCACGTGTGCTTCGATTGGCTGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGaaa                                                            >  1:635165/1‑150 (MQ=255)
                                                                                           gAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGccc                                                                                                                             <  2:1583768/83‑1 (MQ=255)
                                                                                           gAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGccc                                                                                                                             >  1:1583768/1‑83 (MQ=255)
                                                                                           gAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAgg                                                          >  1:2161339/1‑150 (MQ=255)
                                                                                               ttAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCg                                                      >  1:65562/1‑150 (MQ=255)
                                                                                                  aaGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCa                                                   >  2:6549/1‑150 (MQ=255)
                                                                                                    gcgCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGAAAGGGCCg                                                      <  1:2245782/143‑1 (MQ=255)
                                                                                                    gcgCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGAAAGGGCCg                                                      >  2:2245782/3‑145 (MQ=255)
                                                                                                     tGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGGCGGCAAAt                                                >  2:2081820/1‑150 (MQ=255)
                                                                                                        aTCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCt                                             >  2:2038862/1‑150 (MQ=255)
                                                                                                          cTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGc                                           >  1:36476/1‑150 (MQ=255)
                                                                                                           tGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCt                                          >  2:836370/1‑150 (MQ=255)
                                                                                                           tGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCt                                          <  2:36476/150‑1 (MQ=255)
                                                                                                             cAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGATTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCa                                        >  2:935232/1‑150 (MQ=255)
                                                                                                                tGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGcc                                                       <  2:30301/132‑1 (MQ=255)
                                                                                                                tGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGcc                                                       >  1:30301/1‑132 (MQ=255)
                                                                                                                     cGATCAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGCAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGCCCCTTTCGCGGGTACACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCaa                                <  1:649153/150‑1 (MQ=255)
                                                                                                                       aTGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAAc                              >  1:1541106/1‑150 (MQ=255)
                                                                                                                          aGCAAAAGACCCGCCCAGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCgg                                                                          <  2:2225035/103‑1 (MQ=39)
                                                                                                                          aGCAAAAGACCCGCCCAGGCGGGTCTTTTGCTTCTATTTGAGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCgg                                                                          >  1:2225035/1‑103 (MQ=38)
                                                                                                                           gCAAAAGACCCGCCCAGGCGGGTCTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCgg                                                                          >  1:701199/1‑102 (MQ=38)
                                                                                                                           gCAAAAGACCCGCCCAGGCGGGTCTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCgg                                                                          <  2:701199/102‑1 (MQ=38)
                                                                                                                            catccgccccGGCTGCGCGCGTCTTTTGCTTCTATCTGCGGAAACACCCCGCACGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTc                         <  1:2093816/143‑1 (MQ=255)
                                                                                                                            cAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGaaa                                                            >  2:1897260/1‑115 (MQ=255)
                                                                                                                            cAATAGACCCGCCTGGGCGGGTCTTTTGCTTATATATGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGaaa                                                            <  1:1897260/115‑1 (MQ=39)
                                                                                                                             aaaaGACCCGCCCAGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTcc                        >  1:667152/4‑150 (MQ=255)
                                                                                                                             aaaaGACCCGCCCAGGCGGGTCTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTcc                        >  1:1430997/4‑150 (MQ=255)
                                                                                                                              aTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAg                                                                                           <  1:1576082/82‑1 (MQ=255)
                                                                                                                              aTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAg                                                                                           >  2:1576082/1‑82 (MQ=255)
                                                                                                                                    ccgcccaggCGGGTCTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTTCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGa                 >  2:1442971/8‑150 (MQ=255)
                                                                                                                                    ccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGa                 >  2:2293812/1‑150 (MQ=255)
                                                                                                                                      gcccaggCGGGTCTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCCATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGaaa                                                            >  1:1920037/6‑105 (MQ=39)
                                                                                                                                      gcccaggCGGGTCTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGaaa                                                            <  2:1920037/100‑1 (MQ=255)
                                                                                                                                        cTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGGAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCTCTTTCGCGGGTAAACCCGCGAAAGGGCCGTCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCAc             <  1:2081820/150‑1 (MQ=255)
                                                                                                                                        cTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCAc             >  2:751963/1‑150 (MQ=255)
                                                                                                                                               ggTCTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAACGCATCTCCCAGCCGAGCGGTAGCCTTCGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGcc      <  2:1430997/150‑1 (MQ=255)
                                                                                                                                                  cTTCTGCTTCTATTTGTGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCTGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTc   <  2:1588904/150‑1 (MQ=255)
                                                                                                                                                   ttttGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTACGCCAGGGCCGGCCCTTTCGCGCGTAAACCCGCCAAAGGGCCGGCAAATCCTCCTCAGGGCTCAAACCCTTACACAATGATCACCTCAGCCTTCg  >  2:1827999/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCG  >  NC_002947/4338310‑4338606

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