Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,045,978 C→A 5.1% Q262H (CAG→CAT) ‡ 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,9780CA 5.1% 312.8 / 4.5 97Q262H (CAG→CAT) ‡PP_3567LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (44/48);  new base A (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 = 4.15e-01

GCAATGCGGCTTTGCGCACTGCGCAACACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGTGGCCCACTTACGTCCAG  >  NC_002947/4045830‑4046127
                                                                                                                                                    |                                                                                                                                                     
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:1765284/1‑150 (MQ=255)
                    gcgcAACACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGc                                                                                                                                  <  2:1934944/150‑1 (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                                   >  1:2014710/1‑150 (MQ=255)
                                                                                                                   gTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGc                                   <  2:2014710/150‑1 (MQ=255)
                                                                                                                      gATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCAtt                                >  1:49869/1‑150 (MQ=255)
                                                                                                                           tgctATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTcg                                                                                                       >  2:118495/1‑74 (MQ=255)
                                                                                                                           tgctATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTcg                                                                                                       <  1:118495/74‑1 (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)
                                                                                                                                                    |                                                                                                                                                     
GCAATGCGGCTTTGCGCACTGCGCAACACTAGTGACCGTGTATGAGCAAATTTTCGCAAAAGTCGATGAATGCGTTTACCCGCCGTGAGCCCCGGTGGTTGGGCATATACAGGGCGTTGATGCTGCTATTGGCTGCGCCAGGACTGGCCTGCCAGTCGGGGAATATTCGCTGCAGGCGCCCTGCCAACACATCCTCGCGTACCAGCCAGTCGGCCAGCAGTGCGATACCACTGCCTGCCAGTGCCGTTTCACGCAACATGTCGGCATTGGCGCTGCGCAGTGGCCCACTTACGTCCAG  >  NC_002947/4045830‑4046127

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: