Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,869,902 T→G 7.1% E354A (GAG→GCG) ‡ 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,9020TG 7.1% 230.2 / 7.6 86E354A (GAG→GCG) ‡PP_3419sigma‑54 dependent transcriptional regulator/response regulator
Reads supporting (aligned to +/- strand):  ref base T (39/39);  new base G (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 = 8.62e-01

CCAACGGCTCGTCAACCAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGGAAGAACAGGTCACGC  >  NC_002947/3869753‑3870048
                                                                                                                                                     |                                                                                                                                                  
ccAACGGCTAGTCAACCAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGACTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAAACGAGCAGCAGGCTTtcggcct                                                                                                                                                    >  1:879081/1‑150 (MQ=255)
     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                                                                                                                                     <  2:680279/127‑1 (MQ=255)
                                      tGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCa                                                                                                                                     >  1:680279/1‑127 (MQ=255)
                                      tGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCgg                                                                                                              >  1:819119/1‑150 (MQ=255)
                                       gATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGc                                                                                                             <  1:1610810/150‑1 (MQ=255)
                                        aTGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGaaa                                                                                                                                           <  1:100759/119‑1 (MQ=255)
                                        aTGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGaaa                                                                                                                                           >  2:100759/1‑119 (MQ=255)
                                           cgcTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTAtg                                                                                                         <  2:655232/150‑1 (MQ=255)
                                            gcTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCt                                                                                                                            <  1:1916301/130‑1 (MQ=255)
                                            gcTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCt                                                                                                                            >  2:1916301/1‑130 (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                                                                        >  2:80149/1‑143 (MQ=255)
                                                                                   tCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGAt                                                                        <  1:80149/143‑1 (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                          <  2:819119/150‑1 (MQ=255)
                                                                                                                          cTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATg                          <  1:2236521/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               <  2:2077960/142‑1 (MQ=255)
                                                                                                                                       tgcgcttttCGGCCGAGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGgaa               >  1:2077960/7‑148 (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)
                                                                                                                                                     |                                                                                                                                                  
CCAACGGCTCGTCAACCAAGGTCATGCACAGGTTCAACTGATGCGCTTGCACCACCTCGTTCGGCGCCAGCAACACCGTCTGCTCCAGCATGTTGCGCAATTCACGCACATTGCCAGGCCAGCTATAACCGAGCAGCAGGCTTTCGGCCTCGGCCGAAAAGCGCAGGTTCGGCTTGCCATAGCGCCGGCTATGGTGGGCAAGGAAATGCCGCGCCAGTAGCAGGATGTCCTGCCCACGCGAATACAGGCGCGGCACCTTCAGGGCGATAATGCGAAGGCGGAAGAACAGGTCACGC  >  NC_002947/3869753‑3870048

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