Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,045,979 T→G 5.2% Q262P (CAG→CCG) ‡ PP_3567 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,045,9790TG 5.2% 292.8 / 6.5 97Q262P (CAG→CCG) ‡PP_3567LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (44/48);  new base G (3/2);  total (47/50)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-01

GCAATGCGGCTTTGCGCACTGCGCAACACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGTGGCCCACTTACGTCCAGT  >  NC_002947/4045830‑4046128
                                                                                                                                                     |                                                                                                                                                     
gcAATGCGGCTTTGCGCACTGCGCAACACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCt                                                                                                                                                       >  2:90906/1‑150 (MQ=255)
gcAATGCGGCTTTGCGCACTGCGCAACACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCt                                                                                                                                                       >  2:607184/1‑150 (MQ=255)
  aaTGCGGCTTTGCGCACTGCGCAACACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGc                                                                                                                                                     <  2:45411/150‑1 (MQ=255)
       ggCTTTGCGCACTGCGCAACACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTc                                                                                                                                                <  1:2174165/150‑1 (MQ=255)
       ggCTTTGCGCACTGCGCAACACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTc                                                                                                                                                >  2:559336/1‑150 (MQ=255)
        gccttGCGGACTGCGCAACACTAGTGACCGTGGATGAGCAAGTGTTCGCAAAAGTCGAGGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTGGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCg                                                                                                                                               <  1:2040025/147‑1 (MQ=255)
          tttGCGCACTGCGCAACACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTACCAGTCggg                                                                                                                                             <  1:1460929/150‑1 (MQ=255)
                 aCTGCGCAACACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCggg                                                                                                                                             <  2:1926985/143‑1 (MQ=255)
                 aCTGCGCAACACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCggg                                                                                                                                             >  1:1926985/1‑143 (MQ=255)
                 aCTGCGCAACACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATAtt                                                                                                                                      >  2:978343/1‑150 (MQ=255)
                 aCTGCGCAACACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATAtt                                                                                                                                      >  1:1447140/1‑150 (MQ=255)
                   tGCGCAACACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCg                                                                                                                                    <  1:90906/150‑1 (MQ=255)
                    gcgcAACACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGc                                                                                                                                   <  2:1934944/150‑1 (MQ=255)
                    gcgcAACACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGc                                                                                                                                   >  2:1765284/1‑150 (MQ=255)
                     cgcAACACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCt                                                                                                                                  <  1:391144/150‑1 (MQ=255)
                      gcAACACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTg                                                                                                                                 >  2:1034323/1‑150 (MQ=255)
                       cAACACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGc                                                                                                                                <  1:2095307/150‑1 (MQ=255)
                        aaCACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTCGCCTGCCAGTCGGGGAATATTCGCTGCa                                                                                                                               <  1:89595/150‑1 (MQ=255)
                               gTGACCGTGTATTAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGccc                                                                                                                        <  1:584816/150‑1 (MQ=255)
                               gTGACCGTGTATGAGCAAATTTTCCCAAAAGTCGATGGATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCAGATACAGGGCGGTGATGATGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGccc                                                                                                                        <  1:1042933/150‑1 (MQ=255)
                                  aCCGTGTATGAGCAACTTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGc                                                                                                                     >  2:643572/1‑150 (MQ=255)
                                   ccGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGcc                                                                                                                    >  2:2302529/1‑150 (MQ=255)
                                     gtgtATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCaa                                                                                                                  <  1:1765284/150‑1 (MQ=255)
                                     gtgtATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCaa                                                                                                                  <  2:419608/150‑1 (MQ=255)
                                      tgtATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAac                                                                                                                 <  2:1384444/150‑1 (MQ=255)
                                         aTGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAAcaca                                                                                                              <  2:1546881/150‑1 (MQ=255)
                                                aaTTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTcgc                                                                                                       <  1:112779/150‑1 (MQ=255)
                                                aaTTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGGGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTcgc                                                                                                       <  1:1221804/150‑1 (MQ=255)
                                                    ttCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTTCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGcc                                                                                                                    >  2:2396716/1‑133 (MQ=255)
                                                    ttCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTTCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGcc                                                                                                                    <  1:2396716/133‑1 (MQ=255)
                                                       gCAAAAGTCGATGAATGCGTTTACCCGCCGTGGGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTAccag                                                                                                >  2:212499/1‑150 (MQ=255)
                                                       gCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCCGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACCCATCCTCGCGCAccag                                                                                                >  2:469411/1‑150 (MQ=255)
                                                        cAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTAccagc                                                                                               >  1:392168/1‑150 (MQ=255)
                                                         aaaaGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTAccagcc                                                                                              >  2:1507722/1‑150 (MQ=255)
                                                          aaaGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTAccagcca                                                                                             <  1:52113/150‑1 (MQ=255)
                                                          aaaGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTAccagcca                                                                                             <  2:500821/150‑1 (MQ=255)
                                                          aaaGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTAccagcca                                                                                             <  2:2246272/150‑1 (MQ=255)
                                                              tCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCg                                                                                         <  2:1540449/150‑1 (MQ=255)
                                                                gATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGACCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGc                                                                                       >  2:75891/1‑150 (MQ=255)
                                                                gATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATAAAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGc                                                                                       <  1:612772/150‑1 (MQ=255)
                                                                 aTGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGc                                                                                                                                   <  2:1260415/105‑1 (MQ=255)
                                                                 aTGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGc                                                                                                                                   >  1:1260415/1‑105 (MQ=255)
                                                                  tGAATGCGTTTACCCGCCGTGAGCCCGGGTGGTTGGGTATATACAGGGCGTTGATGCTGCTATTGGCAGCGCCAGGACTGGCCTGCCAGTCGGTGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCca                                                                                     <  1:469411/150‑1 (MQ=255)
                                                                   gAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCGGGCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTTGCGTACCAGCCAGTCGGCcag                                                                                    >  1:488963/1‑150 (MQ=255)
                                                                     aTGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTTCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCgagatcgg                                                                                  >  1:674249/1‑143 (MQ=255)
                                                                      tGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCcagcag                                                                                 >  2:1229316/1‑150 (MQ=255)
                                                                         gTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGc                                                                              <  1:196001/150‑1 (MQ=255)
                                                                         gTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGc                                                                              >  2:1541267/1‑150 (MQ=255)
                                                                             aCCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTACAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATa                                                                          <  2:532593/150‑1 (MQ=255)
                                                                               ccgccgTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATAcc                                                                        >  2:31189/1‑150 (MQ=255)
                                                                                cgccgTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCAGGCCAGTcct                                                                                                                                              >  1:2363729/1‑77 (MQ=255)
                                                                                cgccgTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCAGGCCAGTcct                                                                                                                                              <  2:2363729/79‑3 (MQ=255)
                                                                                  ccgTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCAct                                                                     <  1:1034323/150‑1 (MQ=255)
                                                                                        gCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCTCCCTGCCAACCCATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCActgcctgc                                                               <  2:945477/150‑1 (MQ=255)
                                                                                          cccGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCa                                                             <  2:2287173/150‑1 (MQ=255)
                                                                                          cccGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCa                                                             >  1:273660/1‑150 (MQ=255)
                                                                                          cccGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACACCATCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCa                                                             <  1:559336/150‑1 (MQ=255)
                                                                                           ccGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTCCCAg                                                            >  2:651480/1‑150 (MQ=255)
                                                                                            cGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGt                                                           >  2:2025915/1‑150 (MQ=255)
                                                                                            cGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAAGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGt                                                           <  1:978343/150‑1 (MQ=255)
                                                                                               tggtTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGcc                                                        <  1:98910/150‑1 (MQ=255)
                                                                                                      gCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCAc                                                 >  1:196348/1‑150 (MQ=255)
                                                                                                         tataCAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCa                                              >  1:1814865/1‑150 (MQ=255)
                                                                                                             cAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCActgcct                                                                 >  1:300554/1‑127 (MQ=255)
                                                                                                             cAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCActgcct                                                                 <  2:300554/127‑1 (MQ=255)
                                                                                                              aGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGGCAGGCAGTggta                                                       >  1:197210/1‑133 (MQ=255)
                                                                                                                ggCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCg                                                                                         <  1:2241868/100‑1 (MQ=255)
                                                                                                                ggCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCg                                                                                         >  2:406593/1‑100 (MQ=255)
                                                                                                                ggCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCg                                                                                         <  1:406593/100‑1 (MQ=255)
                                                                                                                ggCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCg                                                                                         >  2:2241868/1‑100 (MQ=255)
                                                                                                                   gTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGc                                                                              <  2:1285265/108‑1 (MQ=255)
                                                                                                                   gTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGc                                                                              >  1:1285265/1‑108 (MQ=255)
                                                                                                                   gTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGc                                    <  2:2014710/150‑1 (MQ=255)
                                                                                                                   gTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGc                                    >  1:2014710/1‑150 (MQ=255)
                                                                                                                      gATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCAtt                                 >  1:49869/1‑150 (MQ=255)
                                                                                                                           tgctATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTcg                                                                                                        <  1:118495/74‑1 (MQ=255)
                                                                                                                           tgctATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTcg                                                                                                        >  2:118495/1‑74 (MQ=255)
                                                                                                                           tgctATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGgcgc                            >  1:1971711/1‑150 (MQ=255)
                                                                                                                           tgctATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGgcgc                            >  2:38579/1‑150 (MQ=255)
                                                                                                                            gctATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCt                           <  1:643572/150‑1 (MQ=255)
                                                                                                                              tATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTg                                                                               >  1:950080/1‑96 (MQ=255)
                                                                                                                              tATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTg                                                                               <  2:950080/96‑1 (MQ=255)
                                                                                                                                ttGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTgcgc                       >  1:34673/1‑150 (MQ=255)
                                                                                                                                ttGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTgcgc                       >  2:495436/1‑150 (MQ=255)
                                                                                                                                 tGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCa                      <  1:607184/150‑1 (MQ=255)
                                                                                                                                 tGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCAAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCa                      <  1:2302529/150‑1 (MQ=255)
                                                                                                                                      gcgcCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGTGGc                 >  2:1026450/1‑150 (MQ=255)
                                                                                                                                        gcCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGTGGccc               <  2:392168/150‑1 (MQ=255)
                                                                                                                                          cccgACTGGCAGGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGt                    >  1:104130/4‑143 (MQ=255)
                                                                                                                                          cacgACTGGCAGGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGt                    <  2:104130/140‑1 (MQ=255)
                                                                                                                                          cAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCg                                      >  1:1413463/1‑125 (MQ=255)
                                                                                                                                          cAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCg                                      <  2:1413463/125‑1 (MQ=255)
                                                                                                                                           tcgACTGGCAGGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTCCCAGCCAGTAGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGTGGCCCACt            >  1:1552870/3‑150 (MQ=255)
                                                                                                                                             gACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGTGGCCCACTTa          >  2:1960380/1‑150 (MQ=255)
                                                                                                                                              aCTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGTGGCCCACTTAc         <  1:31189/150‑1 (MQ=255)
                                                                                                                                                  gCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTgcgc                       >  2:1584902/1‑132 (MQ=255)
                                                                                                                                                  gCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTgcgc                       <  1:1584902/132‑1 (MQ=255)
                                                                                                                                                    cTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGTGGCCCACTTACGTCCAg   >  2:1545326/1‑150 (MQ=255)
                                                                                                                                                    cTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGTGGCCCACTTACGTCCAg   >  1:622364/1‑150 (MQ=255)
                                                                                                                                                     tGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGTGGCCCACTTACGTCCAGt  >  2:48602/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAATGCGGCTTTGCGCACTGCGCAACACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGTGGCCCACTTACGTCCAGT  >  NC_002947/4045830‑4046128

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: