Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,246,284 C→T 12.0% intergenic (+30/+61) ibpA → / ← PP_1983 small heat shock protein IbpA/PAS/PAC/GAF sensor‑containing diguanylate cyclase/phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,246,2840CT12.0% 258.8 / 22.7 100intergenic (+30/+61)ibpA/PP_1983small heat shock protein IbpA/PAS/PAC/GAF sensor‑containing diguanylate cyclase/phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (38/50);  new base T (6/6);  total (44/56)
Fisher's exact test for biased strand distribution p-value = 7.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

Only 100 of 104 total aligned reads displayed.

GTGAAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGATG  >  NC_002947/2246138‑2246431
                                                                                                                                                  |                                                                                                                                                     
gTGAAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCaga                                                                                                                                                    >  1:2316310/1‑150 (MQ=255)
       ccgccgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCGAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGccc                                                                                                                                             <  1:674623/150‑1 (MQ=255)
       ccgccgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGccc                                                                                                                                             <  2:1650708/150‑1 (MQ=255)
        cgccgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑t                                                                                                                                          <  1:2858031/150‑1 (MQ=255)
        cgccgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑t                                                                                                                                          <  2:2316310/150‑1 (MQ=255)
         gccgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑tt                                                                                                                                         <  2:143531/150‑1 (MQ=255)
           cgGCCTGGCGAATGGCCTGCTCAGCATCGACCTGCTGCGAATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCC‑‑tttt                                                                                                                                       <  2:723911/150‑1 (MQ=255)
             gCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGc                                                                                                                                               <  1:2090654/142‑1 (MQ=255)
             gCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGc                                                                                                                                               >  2:2090654/1‑142 (MQ=255)
             gCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTc                                                                                                                                     >  2:1146874/1‑150 (MQ=255)
               cTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCC‑‑tttttttc                                                                                                                                   <  2:3011673/150‑4 (MQ=255)
               cTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGt                                                                                                                                   <  2:2180933/150‑1 (MQ=255)
                    cAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCATGAAGAATCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTtc                                                                                                                              <  1:1568200/150‑1 (MQ=255)
                    aaaTGGCCTTTTCAGCATCGGCCTCCACGGCATCGTCTCCGAATAAGGCAATGCCAAGCCTATTCCGATCAATGGCGAAAAACCAGCACAGAATGCATCGAGCTGCAACGCAAAACAAGGGACCGCCACAGGCGCTC‑‑TTTTTCCTGCGtc                                                                                                                              <  2:1110599/149‑1 (MQ=255)
                      aTGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTtctc                                                                                                                            <  1:2811118/150‑1 (MQ=255)
                       tGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCa                                                                                                                           <  2:108191/150‑1 (MQ=255)
                        ggCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCAc                                                                                                                          <  1:1146874/150‑1 (MQ=255)
                           cTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTg                                                                                                                                  <  2:41258/139‑1 (MQ=255)
                           cTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTg                                                                                                                                  >  1:41258/1‑139 (MQ=255)
                               tCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTcc                                                                                                                   >  2:1192510/1‑150 (MQ=255)
                                      cGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAg                                                                                                            <  1:2479830/150‑1 (MQ=255)
                                        aCCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGt                                                                                                          >  2:414661/1‑150 (MQ=255)
                                        aCCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGt                                                                                                          <  2:625011/150‑1 (MQ=255)
                                         ccTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCC‑‑tttt                                                                                                                                       >  1:2937301/1‑120 (MQ=255)
                                         ccTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCC‑‑tttt                                                                                                                                       <  2:2937301/120‑1 (MQ=255)
                                         ccTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTg                                                                                                         >  2:206618/1‑150 (MQ=255)
                                            gctgcGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCa                                                                                                                                                      <  1:1624228/104‑1 (MQ=255)
                                            gctgcGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCa                                                                                                                                                      >  2:1624228/1‑104 (MQ=255)
                                             ctgcGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGt                                                                                                                                   <  1:1419322/120‑1 (MQ=255)
                                             ctgcGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGt                                                                                                                                   >  2:1419322/1‑120 (MQ=255)
                                                  cATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTg                                                                                                                                  <  2:2053791/116‑1 (MQ=255)
                                                  cATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTg                                                                                                                                  >  1:2053791/1‑116 (MQ=255)
                                                   aTCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTg                                                                                                         >  2:1696993/1‑140 (MQ=255)
                                                   aTCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTc                                                                                               <  1:1181306/150‑1 (MQ=255)
                                                   aTCGTGCCTGAAGAAGCCAAGGCCAAGAGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTg                                                                                                         <  1:1696993/140‑1 (MQ=255)
                                                      gTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATg                                                                                            <  1:206618/150‑1 (MQ=255)
                                                      gTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATg                                                                                            <  2:1432809/150‑1 (MQ=255)
                                                       tGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGc                                                                                           >  1:2951023/1‑150 (MQ=255)
                                                          cTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTc                                                                                                                                     >  1:2926805/1‑105 (MQ=255)
                                                          cTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTc                                                                                                                                     <  2:2926805/105‑1 (MQ=255)
                                                          cTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGAAGAAACGAAAAAAAGGGCACCTCCAGAGGTGCCA‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGAt                                                                                        <  2:1553805/150‑1 (MQ=255)
                                                           tGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATc                                                                                       >  1:424929/1‑150 (MQ=255)
                                                              agaagCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGt                                                                                    >  1:1517366/1‑150 (MQ=255)
                                                               gaagCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCGGAGGTGCCC‑‑ttt                                                                                                                                        <  2:199648/97‑1 (MQ=255)
                                                               gaagCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCGGAGGTGCCC‑‑ttt                                                                                                                                        >  1:199648/1‑97 (MQ=255)
                                                                aagCCAAGGCCAAGCGCATTACAATAAATGGCGAAAAACCAGCGCTGAACTGATCGAAAAAAAAAAAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTaa                                                                                  <  2:1348927/150‑1 (MQ=255)
                                                                   ccaaggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGc                                                                                           >  2:130766/1‑138 (MQ=255)
                                                                   ccaaggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGc                                                                                           <  1:130766/138‑1 (MQ=255)
                                                                   ccaaggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCt                                                                               >  1:86357/1‑150 (MQ=255)
                                                                     aaggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAgg                                                                                                                                                  >  1:1478959/1‑83 (MQ=255)
                                                                     aaggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAgg                                                                                                                                                  <  2:1478959/83‑1 (MQ=255)
                                                                     aaggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGc                                                                             <  2:86357/150‑1 (MQ=255)
                                                                     aagcccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGCTGCCTGGCCGTCATGCGATCAGTAATCTGc                                                                             <  1:1024606/150‑1 (MQ=255)
                                                                      aggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGcc                                                                            <  2:1598095/150‑1 (MQ=255)
                                                                           aagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGgcgc                                                                       >  2:890057/1‑150 (MQ=255)
                                                                                cATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCa                                                                  <  1:414661/150‑1 (MQ=255)
                                                                                    ccGATCAATGGCGAAAAACCAGCGCTGAACTGAGCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTCCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCaa                                                              <  2:61423/150‑1 (MQ=255)
                                                                                     cGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCC‑‑t                                                                                                                                          >  1:2227064/1‑73 (MQ=255)
                                                                                     cGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCC‑‑t                                                                                                                                          <  2:2227064/73‑1 (MQ=38)
                                                                                         cAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCt                                                         >  2:509724/1‑150 (MQ=255)
                                                                                          aaTGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTg                                                        >  2:177144/1‑150 (MQ=255)
                                                                                           aTGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCt                                                         <  2:916280/148‑1 (MQ=255)
                                                                                           aTGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCt                                                         >  1:916280/1‑148 (MQ=255)
                                                                                                gAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGaaa                                                  >  2:1862883/1‑150 (MQ=255)
                                                                                                     aCCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcg                                             <  1:1192510/150‑1 (MQ=255)
                                                                                                     aCCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcg                                             <  1:509724/150‑1 (MQ=255)
                                                                                                      ccAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcgg                                            >  1:2510332/1‑150 (MQ=255)
                                                                                                          cgcTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTccc                                                                                                                  <  1:2149798/76‑1 (MQ=255)
                                                                                                           gcTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCTTGCTTCTCACGCAGTCCCGCGTAGGTGCGCgg                                                                                                    <  1:2547670/89‑1 (MQ=255)
                                                                                                           gcTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCgg                                                                                                    >  2:2547670/1‑89 (MQ=255)
                                                                                                                 cTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATg                                 <  2:3071985/150‑1 (MQ=255)
                                                                                                                  tGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATgg                                >  2:2657862/1‑150 (MQ=255)
                                                                                                                  tGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TGTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATgg                                >  1:387614/1‑150 (MQ=255)
                                                                                                                      cGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTCCCGCCGCCTCCAGCAACATTTGCCGAATCTCCGGTGACGGTAATgcgcgg                            >  2:337464/1‑145 (MQ=255)
                                                                                                                         gCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCtg                         >  1:2377330/1‑150 (MQ=255)
                                                                                                                               gcacgaaaaaGGGCACCTCTGGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGa                   >  2:2680782/6‑150 (MQ=255)
                                                                                                                               aCGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGa                   <  1:180530/150‑1 (MQ=255)
                                                                                                                                cacgaaaaaGGGCACCTCTGGAGGTGCCCTTTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAg                    >  1:2001314/5‑150 (MQ=38)
                                                                                                                                    aaaaaGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggtaac                                        >  1:1085077/1‑124 (MQ=255)
                                                                                                                                    aaaaaGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggtaac                                        <  2:1085077/124‑1 (MQ=255)
                                                                                                                                      aaaGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCt            <  1:2657862/150‑1 (MQ=255)
                                                                                                                                      aaaGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCt            >  1:250182/1‑150 (MQ=255)
                                                                                                                                      aaaGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCt            <  1:177144/150‑1 (MQ=255)
                                                                                                                                      aaaGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCGGCCGGCTCCTCCAGCCACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCt            <  1:2041597/150‑1 (MQ=255)
                                                                                                                                      aaaGGGCACCGCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGGGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCt            <  2:424929/150‑1 (MQ=255)
                                                                                                                                       aaGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTg           >  2:1533066/1‑150 (MQ=255)
                                                                                                                                        aGGGCACCTCTGGAGGTGCCCTTTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCt            <  2:2001314/150‑1 (MQ=255)
                                                                                                                                         gggCACCTCTGGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGt         >  2:123530/1‑150 (MQ=255)
                                                                                                                                          ggCACCTCTGGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAg                                                                                                            >  1:435992/1‑50 (MQ=25)
                                                                                                                                          ggCACCTCTGGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAg                                                                                                            <  2:435992/50‑1 (MQ=25)
                                                                                                                                             aCCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAg     >  1:1806326/1‑150 (MQ=255)
                                                                                                                                               cTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACTTCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGAt   <  2:1517366/150‑1 (MQ=255)
                                                                                                                                               cTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCt            >  1:2894548/1‑141 (MQ=255)
                                                                                                                                               cTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCt            <  2:2894548/141‑1 (MQ=255)
                                                                                                                                               cTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGt         >  1:110389/1‑144 (MQ=255)
                                                                                                                                                tCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGATg  >  2:1914887/1‑150 (MQ=255)
                                                                                                                                                  cAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGc                               <  2:2508452/119‑1 (MQ=255)
                                                                                                                                                  cAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGc                               >  1:2508452/1‑119 (MQ=255)
                                                                                                                                                  cAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGa                   >  2:2465103/1‑131 (MQ=255)
                                                                                                                                                  cAGAGGTGCCC‑‑TTTTTCGAGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGa                   <  1:2465103/131‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GTGAAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGATG  >  NC_002947/2246138‑2246431

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: