Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,413,400 T→A 6.9% intergenic (+7/+29) PP_1234 → / ← PP_1235 permease/AhpC/TSA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,413,4000TA 6.9% 274.0 / 9.3 89intergenic (+7/+29)PP_1234/PP_1235permease/AhpC/TSA family protein
Reads supporting (aligned to +/- strand):  ref base T (42/40);  new base A (3/3);  total (46/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

GCGATTATCTGCGCGGTGCTGCTGTTTGGCGGGCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTC  >  NC_002947/1413253‑1413548
                                                                                                                                                   |                                                                                                                                                    
gCGATTATCTGCGCGGTGCTGCTGTTTGGCGGGCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTg                                                                                                                                                    <  1:1827035/150‑1 (MQ=255)
gCGATTATCTGCGCGGTGCTGCTGTTTGGCGGGCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTg                                                                                                                                                    >  2:703802/1‑150 (MQ=255)
          gcgcGGTGCTGCTGTTTGGCGGGCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAg                                                                                                                                          >  1:1516153/1‑150 (MQ=255)
               gTGCTGCTGTTTGGCGGGCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGc                                                                                                                                     >  2:328186/1‑150 (MQ=255)
                 gctgctGTTTGGCGGGCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAAtt                                                                                                                                                     <  1:366635/132‑1 (MQ=255)
                 gctgctGTTTGGCGGGCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAAtt                                                                                                                                                     >  2:366635/1‑132 (MQ=255)
                  ctgctgTTTGGCGGGCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGcccc                                                                                                                                  >  2:1826603/1‑150 (MQ=255)
                  ctgctgTTTGGCGGGCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGcccc                                                                                                                                  <  1:2161937/150‑1 (MQ=255)
                  ctgctgTTTGGCGGGCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGcccc                                                                                                                                  <  1:1282791/150‑1 (MQ=255)
                   tgctgTTTGGCGGGCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCa                                                                                                                                 >  1:1074225/1‑150 (MQ=255)
                   tgctgTTTGGCGGGCTATGGGGGTTTTGGGGGATATTCTTAGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCa                                                                                                                                 >  1:2312342/1‑150 (MQ=255)
                     ctgTTTGGCGGGCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCCTTGCAGCCCCaat                                                                                                                               <  1:1083099/150‑3 (MQ=255)
                     ctgTTTGGCGGGCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAAATGGGGCTGCCTTGCAGCCCCaat                                                                                                                               >  1:507449/1‑148 (MQ=255)
                      tgTTTGGCGGGCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAAATGGGGCTGCCTTGCAGCCCCaatt                                                                                                                              <  2:507449/150‑4 (MQ=255)
                         ttGGCGGGCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAAATGGGGCTGCCTTGCAGCCCCaat                                                                                                                               <  2:716464/146‑3 (MQ=255)
                         ttGGCGGGCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAAATGGGGCTGCCTTGCAGCCCCaat                                                                                                                               >  1:716464/1‑144 (MQ=255)
                         ttGGCGGGCTATGGGGGTTTTGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAAtt                                                                                                                          <  1:413115/150‑1 (MQ=255)
                          tGGCGGGCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAAtt                                                                                                                          <  1:1141654/150‑1 (MQ=255)
                              gggCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCAtt                                                                                                                               >  1:2015298/1‑141 (MQ=255)
                              gggCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCAtt                                                                                                                               <  2:2015298/141‑1 (MQ=255)
                              gggCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCaa                                                                                                                      <  1:328186/150‑1 (MQ=255)
                              gggCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCaa                                                                                                                      >  1:581456/1‑150 (MQ=255)
                               ggCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAg                                                                                                                     <  1:964115/150‑1 (MQ=255)
                                 cTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGcc                                                                                                                   >  1:903084/1‑150 (MQ=255)
                                  tATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCt                                                                                                                  >  2:65332/1‑150 (MQ=255)
                                    tGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTg                                                                                                                >  2:155814/1‑150 (MQ=255)
                                      ggggTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGtt                                                                                                              <  1:1826603/150‑1 (MQ=255)
                                          ttttGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAg                                                                                                          <  2:2148098/150‑1 (MQ=255)
                                          ttttGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAg                                                                                                          <  2:1074225/150‑1 (MQ=255)
                                              gggggATTTTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCtt                                                                                                      <  1:29445/150‑1 (MQ=255)
                                                gggATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTgg                                                                                                    <  2:837653/150‑1 (MQ=255)
                                                 ggATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTCCAAGCCTTGAACAAGGCTTgaa                                                                                                   >  1:656326/1‑148 (MQ=255)
                                                    tatTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTgaaatt                                                                                                <  1:1189288/150‑6 (MQ=255)
                                                          tcGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCa                                                                                          <  1:772613/150‑1 (MQ=255)
                                                             ccATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATATAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATa                                                                                       >  1:41508/1‑150 (MQ=255)
                                                                tCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCaa                                                                                                                      >  1:516800/1‑116 (MQ=255)
                                                                tCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCaa                                                                                                                      <  2:516800/116‑1 (MQ=255)
                                                                 cccGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGc                                                                                   <  2:1516153/150‑1 (MQ=255)
                                                                   cGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGc                                                                                                                    <  2:127165/115‑1 (MQ=255)
                                                                   cGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGc                                                                                                                    >  1:127165/1‑115 (MQ=255)
                                                                   cGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTgg                                                                                                    >  1:1936466/1‑131 (MQ=255)
                                                                   cGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTgg                                                                                                    <  2:1936466/131‑1 (MQ=255)
                                                                          gACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAAATGGGGCTGCCTTGCAGCCCCaat                                                                                                                               >  1:973773/1‑95 (MQ=38)
                                                                          gACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAAATGGGGCTGCCTTGCAGCCCCaat                                                                                                                               <  2:973773/97‑3 (MQ=38)
                                                                               tGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCAtt                                                                                                                               <  2:2028813/92‑1 (MQ=255)
                                                                               tGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCAtt                                                                                                                               >  1:2028813/1‑92 (MQ=255)
                                                                               tGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCgg                                                                     >  2:1706368/1‑150 (MQ=255)
                                                                               tGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGAGTGGGCGGCCGCCAATACGGCATCCACATGCCCgg                                                                     <  1:1133441/150‑1 (MQ=255)
                                                                                  tCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAAttt                                                                                                                         <  1:1295530/95‑1 (MQ=255)
                                                                                  tCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAAttt                                                                                                                         >  2:1295530/1‑95 (MQ=255)
                                                                                  tCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCACCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCa                                                                  >  1:318470/1‑150 (MQ=255)
                                                                                    aaGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACt                                                                >  1:374536/1‑150 (MQ=255)
                                                                                     aGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATc                                                                                >  1:1451748/1‑133 (MQ=255)
                                                                                     aGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATc                                                                                <  2:1451748/133‑1 (MQ=255)
                                                                                      ggCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACgg                                                                                    <  2:926521/128‑1 (MQ=255)
                                                                                      ggCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACgg                                                                                    >  1:926521/1‑128 (MQ=255)
                                                                                                gATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCt                                                                                                                  >  1:1643766/1‑88 (MQ=255)
                                                                                                gATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCt                                                                                                                  <  2:1643766/88‑1 (MQ=255)
                                                                                                gATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGggcggc                                                                                               <  2:1100881/107‑1 (MQ=255)
                                                                                                gATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGggcggc                                                                                               >  1:1100881/1‑107 (MQ=255)
                                                                                                     gTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTc                                               <  1:703802/150‑1 (MQ=255)
                                                                                                       ggCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCt                                             >  1:714504/1‑150 (MQ=255)
                                                                                                        gCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCa                                                                                  <  1:604817/112‑1 (MQ=255)
                                                                                                        gCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCa                                                                                  >  2:604817/1‑112 (MQ=255)
                                                                                                         cccAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACTCCACGCCATTCCTgg                                           <  2:581456/150‑1 (MQ=255)
                                                                                                         cccAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTgg                                           >  1:737858/1‑150 (MQ=255)
                                                                                                          ccAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGgc                                          <  1:738966/150‑1 (MQ=255)
                                                                                                           cAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGgcg                                         >  1:501903/1‑150 (MQ=255)
                                                                                                            aaGCAGGAGCCTAGCGTTTCACCGATGCTTTGACTGTAACTAGCGCTCGAATTCCATCGGAATCAATGTTGATCCTGGGTCAAGGCTTGGCCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGgcgc                                        >  1:1212273/1‑150 (MQ=255)
                                                                                                                       tAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTggg                                                                                                   >  2:673775/1‑80 (MQ=255)
                                                                                                                       tAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTggg                                                                                                   <  1:673775/80‑1 (MQ=255)
                                                                                                                       tAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATc                                                                                <  2:1096675/99‑1 (MQ=255)
                                                                                                                       tAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATc                                                                                >  1:1096675/1‑99 (MQ=255)
                                                                                                                       tAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTc                             >  1:1053042/1‑150 (MQ=255)
                                                                                                                         gCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCCTTAAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGGACTTTCACGCCACGCCATGCCTGGCGCAGCACCGCGTCCt                           >  1:836938/1‑150 (MQ=255)
                                                                                                                         gCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGAACAAGGCTTgaa                                                                                                   <  1:1035785/78‑3 (MQ=39)
                                                                                                                         gCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGAACAAGGCTTgaa                                                                                                   >  2:1035785/1‑76 (MQ=39)
                                                                                                                            tttCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGt                        <  2:1053042/150‑1 (MQ=255)
                                                                                                                            tttCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCCCAGCACCCCGTCCTTGt                        >  2:1301922/1‑150 (MQ=255)
                                                                                                                             ttCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTc                       <  2:714504/150‑1 (MQ=255)
                                                                                                                                 ccGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATc                   <  1:1706368/150‑1 (MQ=255)
                                                                                                                                  cgctGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCAc                                   <  2:1538825/130‑1 (MQ=255)
                                                                                                                                  cGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCAc                                   >  1:1538825/1‑133 (MQ=255)
                                                                                                                                   gATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAg                 >  1:600089/1‑150 (MQ=255)
                                                                                                                                       cTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTTAAGGCTTGGGCGGCCGCCAATACGGCATCCACATTCCCGGGCACTTTCACGCCACAAGAGCCCGGGCGCAGCACCCCGTCCTTGTCGATCAGGAAg             >  1:121957/1‑150 (MQ=255)
                                                                                                                                        tttGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTATGTAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAgg            <  1:2261789/150‑1 (MQ=255)
                                                                                                                                         ttGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGcc                                                                        >  1:1315727/1‑89 (MQ=255)
                                                                                                                                         ttGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGcc                                                                        <  2:1315727/89‑1 (MQ=255)
                                                                                                                                              tGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGtgctgc      >  1:821907/1‑150 (MQ=255)
                                                                                                                                                 aaTTGGGGCTGCAAGGCAGCGCCATTCAATTTAAAGCCTTGTTCAGGACTTGGGCGGCCGGCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCAGCGCGTCCTTGTCGATCAGGAAGGTGCTGCgg  >  2:949766/1‑149 (MQ=255)
                                                                                                                                                  aTTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCAACCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTc  >  2:1308755/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCGATTATCTGCGCGGTGCTGCTGTTTGGCGGGCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTC  >  NC_002947/1413253‑1413548

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: