Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,864,792 T→C 10.2% A28A (GCT→GCC hda → DnaA regulatory inactivator Hda

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,864,7920TC10.2% 296.3 / 18.2 109A28A (GCT→GCChdaDnaA regulatory inactivator Hda
Reads supporting (aligned to +/- strand):  ref base T (48/49);  new base C (6/5);  total (55/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

Only 100 of 114 total aligned reads displayed.

ATAGACTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGCGAGCCCGCTGT  >  NC_002947/1864650‑1864941
                                                                                                                                              |                                                                                                                                                     
ataGACTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCAt                                                                                                                                                >  2:118227/1‑150 (MQ=255)
  aGACTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTg                                                                                                                                              <  1:800322/150‑1 (MQ=255)
  aGACTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTg                                                                                                                                              >  1:2246420/1‑150 (MQ=255)
   gACTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTgg                                                                                                                                             <  2:2246420/150‑1 (MQ=255)
   gACTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGAGCGCTGGAGTGCGCCTGCGCTATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTgg                                                                                                                                             >  1:1861326/1‑150 (MQ=255)
    aCTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTggg                                                                                                                                            <  2:440480/150‑1 (MQ=255)
    aCTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTggg                                                                                                                                            <  1:1833795/150‑1 (MQ=255)
     cTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGt                                                                                                                                           <  1:79769/150‑1 (MQ=255)
                     ccAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTAt                                                                                                                           >  2:2430961/1‑150 (MQ=255)
                     ccAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTAt                                                                                                                           >  1:599009/1‑150 (MQ=255)
                     ccAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTAt                                                                                                                           >  2:310000/1‑150 (MQ=255)
                      cAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCCGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCCGCAGCATTggc                                                                                                                                            <  1:758026/132‑2 (MQ=255)
                      cAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATg                                                                                                                          >  1:97997/1‑150 (MQ=255)
                      cAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCCGCAGCATTggc                                                                                                                                            >  2:758026/1‑131 (MQ=255)
                       aaaGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGc                                                                                                                         <  2:130241/150‑1 (MQ=255)
                          ggCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGaa                                                                                                                      <  2:463926/150‑1 (MQ=255)
                           gCCCCCAGCGGTCCCGTCGACCGCCCGCGTTCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAg                                                                                                                     >  1:1639387/1‑150 (MQ=255)
                              cccAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGc                                                                                                                         >  1:1998268/1‑143 (MQ=255)
                              cccAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGc                                                                                                                         <  2:1998268/143‑1 (MQ=255)
                               ccAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGa                                                                                                                 >  1:93381/1‑150 (MQ=255)
                               ccAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGa                                                                                                                 >  1:1631308/1‑150 (MQ=255)
                               ccAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGa                                                                                                                 >  2:681435/1‑150 (MQ=255)
                                cAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGAc                                                                                                                >  1:815557/1‑150 (MQ=255)
                                 aGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACg                                                                                                               >  1:2042979/1‑150 (MQ=255)
                                      tCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGc                                                                                                          <  2:647473/150‑1 (MQ=255)
                                          gTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGa                                                                                                      <  1:681435/150‑1 (MQ=255)
                                                cgcccgcGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCgaga                                                                                                >  1:495577/1‑150 (MQ=255)
                                                  cccgcGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGc                                                                                              <  2:1523770/150‑1 (MQ=255)
                                                   ccgcGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGcc                                                                                             >  1:225170/1‑150 (MQ=255)
                                                     gcgAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCtt                                                                                           >  1:2288090/1‑150 (MQ=255)
                                                        aGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATc                                                                                        <  1:118227/150‑1 (MQ=255)
                                                        aGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATc                                                                                        >  2:2286798/1‑150 (MQ=255)
                                                           aTGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTAc                                                                                     <  2:1861326/150‑1 (MQ=255)
                                                           aTGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTAc                                                                                     >  1:987327/1‑150 (MQ=255)
                                                                accaccGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCCGCAGCATTggc                                                                                                                                            >  2:1446102/1‑89 (MQ=255)
                                                                accaccGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCCGCAGCATTggc                                                                                                                                            <  1:1446102/90‑2 (MQ=255)
                                                                  caccGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTggg                                                                              <  1:1837593/150‑1 (MQ=255)
                                                                  caccGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTggg                                                                              <  1:2315871/150‑1 (MQ=255)
                                                                   accGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTGTCTACCTGTgggg                                                                             >  2:1296034/1‑150 (MQ=255)
                                                                     cGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAAtgccgca                                                                                                                                                    <  1:1492750/77‑2 (MQ=255)
                                                                       atgcAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAg                                                                         >  2:1676486/4‑150 (MQ=255)
                                                                         ccAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGACGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCa                                                                       <  1:653614/150‑1 (MQ=255)
                                                                           aGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGAc                                                                                                     <  1:143180/118‑1 (MQ=255)
                                                                           aGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGAc                                                                                                     >  2:143180/1‑118 (MQ=255)
                                                                           aGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAgg                                                                     <  2:495577/150‑1 (MQ=255)
                                                                           aGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAgg                                                                     <  2:941123/150‑1 (MQ=255)
                                                                             ttGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGgt                                                                   <  2:2199236/150‑1 (MQ=255)
                                                                                ccccTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAAtgctgc                                                                                                                                                     <  2:1670187/65‑1 (MQ=255)
                                                                                ccccTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAAtgctgc                                                                                                                                                     >  1:1670187/1‑65 (MQ=255)
                                                                                 cccTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGCGGCAAGCAGGGTGTTg                                                               >  2:1220734/1‑150 (MQ=255)
                                                                                  ccTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGtt                                                                                                                                          >  2:1040301/1‑74 (MQ=255)
                                                                                  ccTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGtt                                                                                                                                          <  1:1040301/74‑1 (MQ=255)
                                                                                  ccTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTgg                                                              <  2:2042979/150‑1 (MQ=255)
                                                                                  ccTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTgg                                                              >  2:1694479/1‑150 (MQ=255)
                                                                                   cTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGc                                                             >  2:1822523/1‑150 (MQ=255)
                                                                                   cTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGc                                                             >  2:59363/1‑150 (MQ=255)
                                                                                   cTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGCGTTGGc                                                             >  2:59050/1‑150 (MQ=255)
                                                                                   cTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGCAAGCAGGGTGTTGGcc                                                            <  1:2286798/150‑1 (MQ=255)
                                                                                           gcgcCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCa                                                     <  2:647860/150‑1 (MQ=255)
                                                                                           gcgcCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCa                                                     >  1:352089/1‑150 (MQ=255)
                                                                                             gcCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGcc                                                                                             <  1:1460602/108‑1 (MQ=255)
                                                                                                   gcgATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCtgctgc                                             >  1:2392756/1‑150 (MQ=255)
                                                                                                    cgATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCa                                            <  2:1631308/150‑1 (MQ=255)
                                                                                                       tGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGgc                                         <  2:1639387/150‑1 (MQ=255)
                                                                                                           gacACCTGCATCAACTACTATCCGGGCGCCAATGCAGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGgccgcc                                     >  2:1274770/3‑150 (MQ=255)
                                                                                                            ccACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCa                                    >  2:182738/1‑150 (MQ=255)
                                                                                                            ccACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGAATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCa                                    <  1:2229440/150‑1 (MQ=255)
                                                                                                              aCCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTgg                                                              >  2:905761/1‑122 (MQ=255)
                                                                                                              aCCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTgg                                                              <  1:905761/122‑1 (MQ=255)
                                                                                                                cTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGa                                                                                                                 >  2:419836/1‑69 (MQ=255)
                                                                                                                cTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTg                                                               >  2:2057351/1‑119 (MQ=255)
                                                                                                                cTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTg                                                               <  1:2057351/119‑1 (MQ=255)
                                                                                                                 ttCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTCTCTACCTGTGGGGCAAGCACGGCGTGGGCCGTTCCCACCTGCTGCAGGCCGCCACCCAt                               <  2:744417/150‑1 (MQ=255)
                                                                                                                 ttCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGGGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCAt                               >  2:1031767/1‑150 (MQ=255)
                                                                                                                  tcatcaACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGa                                                                                                                                   >  1:1476507/1‑49 (MQ=255)
                                                                                                                  tcatcaACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGa                                                                                                                                   <  2:1476507/49‑1 (MQ=255)
                                                                                                                     tcaACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCt                           >  1:489407/1‑150 (MQ=255)
                                                                                                                     tcaACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCt                           >  2:916240/1‑150 (MQ=255)
                                                                                                                      caACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCtt                          >  2:670410/1‑150 (MQ=255)
                                                                                                                         ctactaTCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCtt                                                                                           >  2:957891/1‑82 (MQ=255)
                                                                                                                         ctactaTCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCtt                                                                                           <  1:957891/82‑1 (MQ=255)
                                                                                                                         ctactaTCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCa                       <  1:310000/150‑1 (MQ=255)
                                                                                                                          tactaTCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAg                      >  2:843645/1‑150 (MQ=255)
                                                                                                                            ctaTCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGtt                                                         <  1:839501/113‑1 (MQ=255)
                                                                                                                            ctaTCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGtt                                                         >  2:839501/1‑113 (MQ=255)
                                                                                                                              aTCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACctgc                                                >  2:468376/1‑120 (MQ=255)
                                                                                                                              aTCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACctgc                                                <  1:468376/120‑1 (MQ=255)
                                                                                                                                 cGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATc                              <  2:621976/135‑1 (MQ=255)
                                                                                                                                 cGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATc                              >  1:621976/1‑135 (MQ=255)
                                                                                                                                 cGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTg               <  1:182738/150‑1 (MQ=255)
                                                                                                                                   ggCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGc             <  1:916240/150‑1 (MQ=255)
                                                                                                                                     cgcCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGCGa           <  1:2430961/150‑1 (MQ=255)
                                                                                                                                      cccAATGCCGCAGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTgg                                                                               >  2:2240851/2‑81 (MQ=255)
                                                                                                                                      cccAATGCCGCAGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTgg                                                                               <  1:2240851/80‑1 (MQ=255)
                                                                                                                                      cccAATGCCGCAGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCCGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGCGAg          >  2:221633/2‑150 (MQ=255)
                                                                                                                                      cccAATGCCGCAGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGCGAg          >  1:598245/2‑150 (MQ=255)
                                                                                                                                      cccAATGCCGCAGCATTGGGTTACGTCGAGCGGATATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGgcc                                        <  2:892819/119‑1 (MQ=255)
                                                                                                                                         aaTGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGCGAGccc       >  1:1523666/1‑150 (MQ=255)
                                                                                                                                            gctgcGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGCGAGCCCGCt    >  2:18794/1‑150 (MQ=255)
                                                                                                                                              tgcGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGCGAGCCCGCTGt  >  1:600433/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
ATAGACTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGCGAGCCCGCTGT  >  NC_002947/1864650‑1864941

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: