Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,049,346 G→T 8.8% intergenic (‑32/+28) PP_1822 ← / ← folEA‑I hypothetical protein/GTP cyclohydrolase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,049,3460GT 8.8% 201.7 / 7.3 68intergenic (‑32/+28)PP_1822/folEA‑Ihypothetical protein/GTP cyclohydrolase I
Reads supporting (aligned to +/- strand):  ref base G (24/38);  new base T (4/2);  total (28/40)
Fisher's exact test for biased strand distribution p-value = 2.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

GCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGCACATGT  >  NC_002947/2049199‑2049495
                                                                                                                                                   |                                                                                                                                                     
gCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGaa                                                                                                                                                     <  2:1856697/150‑1 (MQ=255)
     cTCGATGCGGGCGTGGGGGGAGAAGTCGCGTTTGCGCTTAGAAGGGCGGCTGATCCAGTCAACCAAGACGAGGAGCCGGCAGAGTAAAACACGCAGAGCTTTGACGATCATGTGAAAAAAAAAAAAAAAAAAAGCCGACCCTGAAGGGTc                                                                                                                                                <  2:688363/150‑1 (MQ=255)
     cTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTc                                                                                                                                                <  1:1654757/150‑1 (MQ=255)
     cTCGATGCGGGCCTGGGCGGCGGCGTCGCGCTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTTCAGGGTc                                                                                                                                                >  2:1478904/1‑150 (MQ=255)
        gATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATACCGTCACCGAAGACGATGAGCTGGCCGATGTCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGc                                                                                                                                             <  2:1160248/150‑1 (MQ=255)
           gCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCg                                                                                                                                          <  2:1834845/150‑1 (MQ=255)
            cGGGCCTGGGCGGCGGCGTCGCGCTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTTCAGGGTCGGCtttt                                                                                                                                         <  1:1478904/150‑4 (MQ=255)
                 cTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGGAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCa                                                                                                                                    <  2:1090905/150‑1 (MQ=255)
                     gcggcggcgTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCaa                                                                                                                                >  1:352651/1‑150 (MQ=255)
                            cgTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCa                                                                                                                         <  2:1557552/150‑1 (MQ=255)
                            cgTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCa                                                                                                                         <  2:997541/150‑1 (MQ=255)
                               cgcgTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACtt                                                                                                                      <  1:548959/150‑1 (MQ=255)
                                   tttGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATc                                                                                                                  >  1:706044/1‑150 (MQ=255)
                                    ttGCGGGTGGCCGGGCGGCTGATGCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGTGCCAGCATGCAATCACTTGATCa                                                                                                                 <  2:274167/150‑1 (MQ=255)
                                      gcgcTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAgg                                                                                                               <  1:67224/150‑1 (MQ=255)
                                           tGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAg                                                                                                          <  1:1502089/150‑1 (MQ=255)
                                           tGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAg                                                                                                          <  2:1875946/150‑1 (MQ=255)
                                           tGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCAAGCAACCAATCACTTGATCAGGCTGAg                                                                                                          <  1:1204062/150‑1 (MQ=255)
                                              ccGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGaa                                                                                                       <  2:1272263/150‑1 (MQ=255)
                                               cGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGACTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAAc                                                                                                      >  1:1687483/1‑150 (MQ=255)
                                                  gcggcTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCg                                                                                                   <  2:1112634/150‑1 (MQ=255)
                                                                cACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTggaggc                                                                                     >  1:940350/1‑150 (MQ=255)
                                                                 aCCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGgcggcg                                                                                    >  2:1718047/1‑150 (MQ=255)
                                                                   tgatgtcGATGAGCCGGCCGAGGCCAGCCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGtt                                                                                  <  1:1834527/144‑1 (MQ=255)
                                                                   cGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGtt                                                                                  <  1:1174175/150‑1 (MQ=255)
                                                                   cGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGACTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGACCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGtt                                                                                  >  2:30896/1‑150 (MQ=255)
                                                                     aaGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGt                                                                                <  2:706044/150‑1 (MQ=255)
                                                                           aTGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGa                                                                          <  2:940350/150‑1 (MQ=255)
                                                                           aTGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGa                                                                          <  2:62569/150‑1 (MQ=255)
                                                                            tGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGaa                                                                         >  2:342705/1‑150 (MQ=255)
                                                                                 tGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGGTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCAc                                                                    <  2:1265066/150‑1 (MQ=255)
                                                                                 gggCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCAc                                                                    <  2:453694/149‑1 (MQ=255)
                                                                                  ggCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCAcc                                                                   >  2:368979/1‑150 (MQ=255)
                                                                                   gCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAg                                                                                                                <  1:183884/104‑1 (MQ=255)
                                                                                   gCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAg                                                                                                                >  2:183884/1‑104 (MQ=255)
                                                                                    ccGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCa                                                                 <  2:352651/150‑1 (MQ=255)
                                                                                      gAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGc                                                               >  2:1398846/1‑150 (MQ=255)
                                                                                           cAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAAc                                                                                                                               <  1:795577/81‑1 (MQ=255)
                                                                                           cAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAAc                                                                                                                               >  2:795577/1‑81 (MQ=255)
                                                                                           cAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCAc                                                          >  2:1216824/1‑150 (MQ=255)
                                                                                              cccGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGa                                                       <  1:342705/150‑1 (MQ=255)
                                                                                               ccGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAg                                                      >  2:1842837/1‑150 (MQ=255)
                                                                                               ccGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGCCGTTGTTGTCGCGGAACTCCCCCAGCATCACCGAg                                                      <  1:1808797/150‑1 (MQ=255)
                                                                                                       cTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGt                                                                                >  2:1584461/1‑116 (MQ=255)
                                                                                                       cTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGt                                                                                <  1:1584461/116‑1 (MQ=255)
                                                                                                          tGACGATCATGGGCCCCTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGa                                                                                                           <  1:1471076/86‑1 (MQ=255)
                                                                                                          tGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGa                                                                                                           >  2:1471076/1‑86 (MQ=255)
                                                                                                              gATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGaa                                       >  2:472862/1‑150 (MQ=255)
                                                                                                                tCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGa                                        >  1:1089970/1‑147 (MQ=255)
                                                                                                                tCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGa                                        <  2:1089970/147‑1 (MQ=255)
                                                                                                                         aCTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTtctc                            <  1:1191733/150‑1 (MQ=255)
                                                                                                                         aCTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTtctc                            <  1:1718047/150‑1 (MQ=255)
                                                                                                                         aCTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTtctc                            <  2:348976/150‑1 (MQ=255)
                                                                                                                         aCTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTtctc                            <  2:366057/150‑1 (MQ=255)
                                                                                                                                    gCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGCACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGGATCACCGAGGTGATCATGGTCGAATTCTGCTTATCGACACGGCGca                 >  2:92230/1‑150 (MQ=255)
                                                                                                                                      aaagGCCGACCCTGAAGGGCCCGCTCGGGGCCGGCAACCAATCACTTGATCAGGCTGGGGAACTCACTGCGGGCGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATCGTCGAATTCTTCTTCTCGACACCGCGcatc               >  2:1036811/1‑150 (MQ=255)
                                                                                                                                      aaaaGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatc               >  2:129088/1‑150 (MQ=255)
                                                                                                                                      aaaaGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatc               >  1:1855660/1‑150 (MQ=255)
                                                                                                                                        gaGCCGACCCTTCAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatcat             >  2:1704403/2‑150 (MQ=255)
                                                                                                                                        aaGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCa                                                                                                                 >  2:1391747/1‑50 (MQ=255)
                                                                                                                                        aaGACGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCa                                                                                                                 <  1:1391747/50‑1 (MQ=255)
                                                                                                                                           ccGACCCTTCAGGGTCGGCTTGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatcatcat          >  2:865984/1‑150 (MQ=255)
                                                                                                                                            cGACCCTTCAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGt                                           >  1:827209/1‑116 (MQ=255)
                                                                                                                                            cGACCCTTCAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGt                                           <  2:827209/116‑1 (MQ=255)
                                                                                                                                             gACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGAGACTGCGTTGTCGCGGAAACCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGCCACCGCGCATCATCATGc        >  1:1077256/1‑150 (MQ=255)
                                                                                                                                              aCCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGca       <  1:368979/150‑1 (MQ=255)
                                                                                                                                                ccTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGcaca     >  1:369011/1‑150 (MQ=255)
                                                                                                                                                 cTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCTGGTGGCGCAGTTGTCGCGTTACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGCACAt    <  2:1402487/150‑1 (MQ=255)
                                                                                                                                                   gAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGCACAtgt  <  2:1855660/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGCACATGT  >  NC_002947/2049199‑2049495

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: