Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,413,411 A→T 7.6% intergenic (+18/+18) 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,4110AT 7.6% 273.6 / 11.8 92intergenic (+18/+18)PP_1234/PP_1235permease/AhpC/TSA family protein
Reads supporting (aligned to +/- strand):  ref base A (46/39);  new base T (3/4);  total (49/43)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

GCGCGGTGCTGCTGTTTGGCGGGCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGT  >  NC_002947/1413263‑1413559
                                                                                                                                                    |                                                                                                                                                    
gcgcGGTGCTGCTGTTTGGCGGGCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAg                                                                                                                                                     >  1:1516153/1‑150 (MQ=255)
     gTGCTGCTGTTTGGCGGGCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGc                                                                                                                                                >  2:328186/1‑150 (MQ=255)
        ctgctgTTTGGCGGGCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGcccc                                                                                                                                             <  1:2161937/150‑1 (MQ=255)
        ctgctgTTTGGCGGGCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGcccc                                                                                                                                             >  2:1826603/1‑150 (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                                                                                                                                          >  1:716464/1‑144 (MQ=255)
               ttGGCGGGCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAAATGGGGCTGCCTTGCAGCCCCaat                                                                                                                                          <  2:716464/146‑3 (MQ=255)
               ttGGCGGGCTATGGGGGTTTTGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAAtt                                                                                                                                     <  1:413115/150‑1 (MQ=255)
                tGGCGGGCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAAtt                                                                                                                                     <  1:1141654/150‑1 (MQ=255)
                    gggCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCAtt                                                                                                                                          <  2:2015298/141‑1 (MQ=255)
                    gggCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCAtt                                                                                                                                          >  1:2015298/1‑141 (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                                                                                                                                          <  2:973773/97‑3 (MQ=38)
                                                                gACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAAATGGGGCTGCCTTGCAGCCCCaat                                                                                                                                          >  1:973773/1‑95 (MQ=38)
                                                                     tGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCAtt                                                                                                                                          >  1:2028813/1‑92 (MQ=255)
                                                                     tGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCAtt                                                                                                                                          <  2:2028813/92‑1 (MQ=255)
                                                                     tGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCgg                                                                                >  2:1706368/1‑150 (MQ=255)
                                                                     tGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGAGTGGGCGGCCGCCAATACGGCATCCACATGCCCgg                                                                                <  1:1133441/150‑1 (MQ=255)
                                                                        tCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAAttt                                                                                                                                    >  2:1295530/1‑95 (MQ=255)
                                                                        tCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAAttt                                                                                                                                    <  1:1295530/95‑1 (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                                                                                               >  1:926521/1‑128 (MQ=255)
                                                                            ggCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACgg                                                                                               <  2:926521/128‑1 (MQ=255)
                                                                                      gATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCt                                                                                                                             <  2:1643766/88‑1 (MQ=255)
                                                                                      gATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCt                                                                                                                             >  1:1643766/1‑88 (MQ=255)
                                                                                      gATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGggcggc                                                                                                          >  1:1100881/1‑107 (MQ=255)
                                                                                      gATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGggcggc                                                                                                          <  2:1100881/107‑1 (MQ=255)
                                                                                           gTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTc                                                          <  1:703802/150‑1 (MQ=255)
                                                                                             ggCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCt                                                        >  1:714504/1‑150 (MQ=255)
                                                                                              gCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCa                                                                                             >  2:604817/1‑112 (MQ=255)
                                                                                              gCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCa                                                                                             <  1:604817/112‑1 (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)
                                                                                                                                       aaTTGGGGCTGCAAGGCAGCGCCATTCAATTTAAAGCCTTGTTCAGGACTTGGGCGGCCGGCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCAGCGCGTCCTTGTCGATCAGGAAGGTGCTGCggg              >  2:949766/1‑149 (MQ=255)
                                                                                                                                        aTTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCAACCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTc             >  2:1308755/1‑150 (MQ=255)
                                                                                                                                                tGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCa     >  1:1546056/1‑150 (MQ=255)
                                                                                                                                                 gCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCAt    <  2:903084/150‑1 (MQ=255)
                                                                                                                                                  cAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCGGGCAGGAGCTGGGGTCAACGCCCATg   >  2:1299864/1‑150 (MQ=255)
                                                                                                                                                  cAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATg   >  1:1739359/1‑150 (MQ=255)
                                                                                                                                                  cAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATg   >  1:802798/1‑150 (MQ=255)
                                                                                                                                                   aaGGCAGCCCCATTAAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGt  >  2:585118/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCGCGGTGCTGCTGTTTGGCGGGCTATGGGGGTTTTGGGGGATATTCTTCGCCATCCCGCTGGCGACGCTGATCAAGGCCGTGCTGGATGCGTGGCCCAAGCAGGAGCCTAGCGTTTCACCGATGCTTTGAATGTAATTGGGGCTGCAAGGCAGCCCCATTCAATTTCAAGCCTTGTTCAAGGCTTGGGCGGCCGCCAATACGGCATCCACATGCCCGGGCACTTTCACGCCACGCCATTCCTGGCGCAGCACCCCGTCCTTGTCGATCAGGAAGGTGCTGCGGTCAACGCCCATGT  >  NC_002947/1413263‑1413559

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: