Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,869,903 C→A 6.9% E354* (GAG→TAG) ‡ PP_3419 ← sigma‑54 dependent transcriptional regulator/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,869,9030CA 6.9% 263.7 / 5.5 86E354* (GAG→TAG) ‡PP_3419sigma‑54 dependent transcriptional regulator/response regulator
Reads supporting (aligned to +/- strand):  ref base C (39/41);  new base A (4/2);  total (43/43)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

GGCTCGTCAACCAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGGAAGAACAGGTCACGC  >  NC_002947/3869758‑3870048
                                                                                                                                                 |                                                                                                                                                 
ggCTCGTCAACCAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTtcggcctcggcc                                                                                                                                               <  2:1442991/150‑1 (MQ=255)
 gCTCGTCAACCAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCg                                                                                                                                              >  1:265087/1‑150 (MQ=255)
  cTCGTCAACCAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGa                                                                                                                                             >  1:2223923/1‑150 (MQ=255)
      tCAACCAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAg                                                                                                                                         <  1:2425446/150‑1 (MQ=255)
         aCCAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAAGTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAgcgc                                                                                                                                      >  2:727605/1‑150 (MQ=255)
          ccAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCa                                                                                                                                     >  2:1610810/1‑150 (MQ=255)
           cAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAg                                                                                                                                    >  2:1113133/1‑150 (MQ=255)
               gTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTc                                                                                                                                <  1:1173458/150‑1 (MQ=255)
                     cacaGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTg                                                                                                                          >  1:1305685/1‑150 (MQ=255)
                          gTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCGGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATa                                                                                                                     >  2:1965074/1‑150 (MQ=255)
                          gTTCAACTGATGCGCTTGCACCACCTCGTGCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGGCAtc                                                                                                                     >  2:484218/1‑149 (MQ=255)
                                 tGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCa                                                                                                                                     >  1:680279/1‑127 (MQ=255)
                                 tGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCa                                                                                                                                     <  2:680279/127‑1 (MQ=255)
                                 tGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCgg                                                                                                              >  1:819119/1‑150 (MQ=255)
                                  gATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGc                                                                                                             <  1:1610810/150‑1 (MQ=255)
                                   aTGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGaaa                                                                                                                                           >  2:100759/1‑119 (MQ=255)
                                   aTGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGaaa                                                                                                                                           <  1:100759/119‑1 (MQ=255)
                                      cgcTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTAtg                                                                                                         <  2:655232/150‑1 (MQ=255)
                                       gcTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCt                                                                                                                            >  2:1916301/1‑130 (MQ=255)
                                       gcTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCt                                                                                                                            <  1:1916301/130‑1 (MQ=255)
                                         ttGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGTGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTAtggtg                                                                                                      <  2:620946/150‑1 (MQ=255)
                                         ttGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTAtggtg                                                                                                      <  2:640993/150‑1 (MQ=255)
                                                  cTCGTTCGGCGCCAGCAACACCGTCTGCTCCATCATGGTGCTCAATCCCTACCCATTTCCAGGCCAGCTATAACCGAGCAGCAGTCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGaa                                                                                             >  2:1165425/1‑150 (MQ=255)
                                                  cTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATCGCGCCGGCTATGGTGGGCaa                                                                                                 <  1:1209351/146‑1 (MQ=255)
                                                  cTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCaa                                                                                                 >  2:1209351/1‑146 (MQ=255)
                                                  cTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGaa                                                                                             >  1:2271353/1‑150 (MQ=255)
                                                  cTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGaa                                                                                             >  1:2327637/1‑150 (MQ=255)
                                                  cTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGaa                                                                                             >  1:648472/1‑150 (MQ=255)
                                                    cGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAAt                                                                                           <  1:2129181/150‑1 (MQ=255)
                                                              cAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTAtggtgg                                                                                                     >  1:2197660/1‑130 (MQ=255)
                                                              cAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTAtggtgg                                                                                                     <  2:2197660/130‑1 (MQ=255)
                                                                   acacCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCa                                                                            >  2:1597263/1‑150 (MQ=255)
                                                                         tCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGt                                                                      <  1:1113133/150‑1 (MQ=255)
                                                                          cTGCGCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTc                                                                     <  1:904848/150‑1 (MQ=255)
                                                                              tCCAGCATGTTGCGCATTTCACGCCCATTGCCAGGCGAGGTATAACCGAGCAGCAGGCTTTCGGCCCCCAACTAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGc                                                                 <  1:413496/150‑1 (MQ=255)
                                                                              tCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGc                                                                                         >  1:2375319/1‑126 (MQ=255)
                                                                              tCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGc                                                                                         <  2:2375319/126‑1 (MQ=255)
                                                                              tCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGGGCCAGTAGCAGGATGTCCTGc                                                                 >  1:2120389/1‑150 (MQ=255)
                                                                              tCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGAt                                                                        <  1:80149/143‑1 (MQ=255)
                                                                              tCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGAt                                                                        >  2:80149/1‑143 (MQ=255)
                                                                                cAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCtt                                                                                                                           >  1:1784515/1‑90 (MQ=255)
                                                                                cAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCtt                                                                                                                           <  2:1784515/90‑1 (MQ=255)
                                                                                    aTGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCAcgc                                                           <  2:653064/150‑1 (MQ=255)
                                                                                       ttGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGTCCTCGGCCGAAAAGCGCAGGTTCGGCtt                                                                                                                           <  2:1250323/83‑1 (MQ=255)
                                                                                       ttGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCtt                                                                                                                           >  1:1250323/1‑83 (MQ=255)
                                                                                         gcgcAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTACCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATa                                                      <  2:2109901/150‑1 (MQ=255)
                                                                                          cgcAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATAc                                                     >  2:1826809/1‑150 (MQ=255)
                                                                                           gcAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAgg                                                                                               >  2:2306500/1‑107 (MQ=255)
                                                                                           gcAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAgg                                                                                               <  1:2306500/107‑1 (MQ=255)
                                                                                           gcAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTAGGACTCTGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACa                                                    >  2:926693/1‑150 (MQ=255)
                                                                                            cAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAg                                                   <  2:879081/150‑1 (MQ=255)
                                                                                            cAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAg                                                   <  1:727605/150‑1 (MQ=255)
                                                                                                tCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTACGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGgcgc                                               <  2:1305685/150‑1 (MQ=255)
                                                                                                 cACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGcgcg                                              <  1:1162473/150‑1 (MQ=255)
                                                                                                  aCGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGc                                                                                                                         <  2:8988/74‑1 (MQ=255)
                                                                                                  aCGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGc                                                                                                                         >  1:8988/1‑74 (MQ=255)
                                                                                                  aCGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCgg                                             >  1:2241393/1‑150 (MQ=255)
                                                                                                  aCGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCgg                                             >  2:2230213/1‑150 (MQ=255)
                                                                                                   cGCACATTGCCAGGCCAGCTATAACCGAGCAGCCGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGGATACAGGCGCGGc                                            >  1:880494/1‑150 (MQ=255)
                                                                                                        aTTGCCAGGCCAGCTATAACCGGGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCtt                                       <  2:1850893/150‑1 (MQ=255)
                                                                                                         ttGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCCGTAGCAGGATGTCCTGCCCACGCGAATACAGgcgc                                               <  1:54321/141‑1 (MQ=255)
                                                                                                         ttGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGgcgc                                               >  2:54321/1‑141 (MQ=255)
                                                                                                            ccaggccagCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGc                                                                 >  1:367294/1‑120 (MQ=255)
                                                                                                            ccaggccagCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGc                                                                 <  2:367294/120‑1 (MQ=255)
                                                                                                             caggccagCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAggg                                  >  1:804424/1‑150 (MQ=255)
                                                                                                              aggccagCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGc                                 >  2:780230/1‑150 (MQ=255)
                                                                                                                  cagCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATa                             >  2:541689/1‑150 (MQ=255)
                                                                                                                     cTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATg                          <  1:2236521/150‑1 (MQ=255)
                                                                                                                     cTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATg                          <  2:819119/150‑1 (MQ=255)
                                                                                                                      tataACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGTTCAGGTTCGGCTTGCCATAGCGCCGGCTAAGATTGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGc                         <  1:413127/150‑1 (MQ=255)
                                                                                                                      tataACCGAGCAGCAAGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGc                         <  1:1813040/150‑1 (MQ=255)
                                                                                                                                cagcagGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGgaa               <  1:780230/150‑1 (MQ=255)
                                                                                                                                cagcagGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGgaa               <  2:2241393/150‑1 (MQ=255)
                                                                                                                                cagcagGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGgaa               <  1:1597263/150‑1 (MQ=255)
                                                                                                                                  tgcgcttttCGGCCGAGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGgaa               >  1:2077960/7‑148 (MQ=255)
                                                                                                                                  tgcgcttttCGGCCGAGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGgaa               <  2:2077960/142‑1 (MQ=255)
                                                                                                                                   cagGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGgaagaa            <  2:2271353/150‑1 (MQ=255)
                                                                                                                                    agGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGGAAGAAc           <  2:2029349/150‑1 (MQ=255)
                                                                                                                                     ggCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGGAAGAACa          >  2:2221684/1‑150 (MQ=255)
                                                                                                                                       ttttCGGCCGAGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACTGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGGAAGAACAgg        >  1:2090411/2‑150 (MQ=255)
                                                                                                                                       ttttCGGCCGAGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGGAAGAACAgg        >  1:2434845/2‑150 (MQ=255)
                                                                                                                                       cTTTCGGCCACGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTGGCAGGATGTCCTGCCCCCGCGAATACAGGCGCGGCACCTTCATGGCGATAATGCGAAGGCGGAAGAACAgg        <  1:1965074/150‑1 (MQ=255)
                                                                                                                                        tttCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCAcgc                                                           <  2:1866805/98‑1 (MQ=255)
                                                                                                                                        tttCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCAcgc                                                           >  1:1866805/1‑98 (MQ=255)
                                                                                                                                        tttCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGGAAGAACAGGt       >  2:2217726/1‑150 (MQ=255)
                                                                                                                                        tttCGGCCGAGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCAAATACAGGCGCGGCACCTTCAGGGCGATAATg                          <  1:1352061/131‑1 (MQ=255)
                                                                                                                                        tttCGGCCGAGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCAAATACAGGCGCGGCACCTTCAGGGCGATAATg                          >  2:1352061/1‑131 (MQ=255)
                                                                                                                                             gcctcggccGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGGAAGAACAGGTCACGc  <  2:648472/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GGCTCGTCAACCAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGGAAGAACAGGTCACGC  >  NC_002947/3869758‑3870048

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: