Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 571,621 T→C 18.5% intergenic (+44/+47) ssb → / ← PP_0486 single‑stranded DNA‑binding protein/GntR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947571,6210TC18.5% 175.0 / 31.3 81intergenic (+44/+47)ssb/PP_0486single‑stranded DNA‑binding protein/GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (33/33);  new base C (8/7);  total (41/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

GCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCGAT  >  NC_002947/571472‑571763
                                                                                                                                                     |                                                                                                                                              
gcagGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGt                                                                                                                                                >  1:1343983/1‑150 (MQ=255)
    gCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGcc                                                                                                                                            <  1:2205013/150‑1 (MQ=255)
     cccAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGccc                                                                                                                                           <  2:2206246/150‑1 (MQ=255)
         gcgcCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTggg                                                                                                                                       >  2:944591/1‑150 (MQ=255)
         gcgcCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAACTGCCCTggg                                                                                                                                       <  2:1379229/150‑1 (MQ=255)
           gcCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGCGCTCAAACCCAGGGCAGTTGCCCTgggct                                                                                                                                     >  1:1210538/1‑148 (MQ=255)
              cGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGttttt                                                                                                                                  >  1:2237053/1‑150 (MQ=255)
               gCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGtttttt                                                                                                                                 <  1:1177166/150‑1 (MQ=255)
               gCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAACTGCCCTGGGtttctg                                                                                                                                 <  1:2193285/150‑4 (MQ=255)
                 aaCAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGccc                                                                                                                                           <  1:480532/138‑1 (MQ=255)
                 aaCAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGccc                                                                                                                                           >  2:480532/1‑138 (MQ=255)
                   cAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATg                                                                                                                             <  2:1431167/150‑1 (MQ=255)
                    aGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATgg                                                                                                                            <  1:895824/150‑1 (MQ=255)
                    aGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATgg                                                                                                                            <  2:2091088/150‑1 (MQ=255)
                        cccTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGgcaa                                                                                                                        >  1:93822/1‑149 (MQ=255)
                          cTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGgcagca                                                                                                                      >  2:2170724/1‑150 (MQ=255)
                          cTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGgcagca                                                                                                                      >  1:356151/1‑150 (MQ=255)
                             agcgcccagcgccGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTTTGGGCAGCACcg                                                                                                                   >  2:1226187/1‑150 (MQ=255)
                                     gcgccGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAg                                                                                                           <  2:6075/150‑1 (MQ=255)
                                     gcgccGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAg                                                                                                           <  1:1762417/150‑1 (MQ=255)
                                          gcAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGCTTTTTTATGACGTGCGCCGCGTATCACGCATc                                                                                                      >  1:1329336/1‑150 (MQ=255)
                                                  cTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAg                                                                                              >  1:2379418/1‑150 (MQ=255)
                                                    gcgcCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTa                                                                                            <  1:2109550/150‑1 (MQ=255)
                                                        cgcAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGaac                                                                                        >  1:2079251/1‑150 (MQ=255)
                                                              cggcggCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCgg                                                                                  <  2:2237053/150‑1 (MQ=255)
                                                                  ggcTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCt                                                                              >  1:67526/1‑150 (MQ=255)
                                                                  ggcTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCt                                                                              >  1:1964293/1‑150 (MQ=255)
                                                                     tGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAACTGCCCTGGGtttctg                                                                                                                                 <  1:76833/96‑4 (MQ=38)
                                                                     tGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAACTGCCCTGGGtttctg                                                                                                                                 >  2:76833/1‑93 (MQ=38)
                                                                          tCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGaa                                                                      <  1:944591/150‑1 (MQ=255)
                                                                               aGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCcag                                                                 <  1:592004/150‑1 (MQ=255)
                                                                               aGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCcag                                                                 <  1:100925/150‑1 (MQ=255)
                                                                                gCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCcagc                                                                <  2:2379418/150‑1 (MQ=255)
                                                                                 cTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCa                                                                                                            <  2:816507/105‑1 (MQ=255)
                                                                                 cTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCa                                                                                                            >  1:816507/1‑105 (MQ=255)
                                                                                  ttCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCcagcag                                                              >  1:2322620/1‑150 (MQ=255)
                                                                                   tCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAgg                                                             <  1:2170724/150‑1 (MQ=255)
                                                                                                 ccGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGt                                               >  2:1781260/1‑150 (MQ=255)
                                                                                                    ttCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAg                                            <  1:1781260/150‑1 (MQ=255)
                                                                                                     tCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCcag                                                                 <  2:1490071/128‑1 (MQ=255)
                                                                                                     tCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCcag                                                                 >  1:1490071/1‑128 (MQ=255)
                                                                                                       tGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGc                                         >  2:377555/1‑150 (MQ=255)
                                                                                                          tCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCaaa                                      >  1:1867914/1‑150 (MQ=255)
                                                                                                            agagCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCGGGGCAAAGc                                    >  1:1814260/1‑150 (MQ=255)
                                                                                                            agagCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGc                                    >  1:1705546/1‑150 (MQ=255)
                                                                                                             gagCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGcc                                   <  2:37972/150‑1 (MQ=255)
                                                                                                                     gggCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCtt                           >  2:2067522/1‑150 (MQ=255)
                                                                                                                      ggCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTg                          <  1:2067522/150‑1 (MQ=255)
                                                                                                                         gTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTggcg                       <  2:1964293/150‑1 (MQ=255)
                                                                                                                              gCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCcagcag                                                              <  2:2273986/106‑1 (MQ=255)
                                                                                                                              gCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCcagcag                                                              >  1:2273986/1‑106 (MQ=255)
                                                                                                                               cAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTggg                 <  1:1032433/150‑1 (MQ=255)
                                                                                                                               cAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTggg                 >  2:1032433/1‑150 (MQ=255)
                                                                                                                                aaGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGt                <  2:356151/150‑1 (MQ=255)
                                                                                                                                  gACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTg                                                 <  1:1809685/115‑1 (MQ=255)
                                                                                                                                  gACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTg                                                 >  2:1809685/1‑115 (MQ=255)
                                                                                                                                    ccAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCa            <  1:1326070/150‑1 (MQ=255)
                                                                                                                                     cAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCtt                           >  2:1691205/1‑134 (MQ=255)
                                                                                                                                     cAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCtt                           <  1:1691205/134‑1 (MQ=255)
                                                                                                                                     cAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCAt           >  1:792586/1‑150 (MQ=255)
                                                                                                                                      aGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTACGGAACGCCTCGCGTATTAGGAATCGCTATCAGTAGATCAACCTGGGCTCCGGCAAATCCAGCAGGCTTGCCAGCTTGTTCAGGGCAAAGCCCAGGTCTTGGCGGCTGGGTGCCATc          >  1:1252674/1‑150 (MQ=255)
                                                                                                                                       gAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCa         <  1:2069376/150‑1 (MQ=255)
                                                                                                                                       aaaaCCCAGGGCAACTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGc                                                                                                    <  2:195186/58‑1 (MQ=25)
                                                                                                                                       aaaaCCCAGGGCAACTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGc                                                                                                    >  1:195186/2‑59 (MQ=25)
                                                                                                                                         agtccAGGGCGGTTGCCCTGGGTTTTTTGATGGGCAGCACTGCGTATCAGGCATCGCTTTCAGTAGAGCAACCGGTGCGCCGGCGAATCCAGCAGGCTTGCCAGGGTGGTCAGGGCAAAGCCGAGTTCTTGGCGGCTGGGTGCCATCAgg       <  2:322689/147‑1 (MQ=255)
                                                                                                                                         aaCCCAGGGCAACTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCa                                                      >  2:2076206/1‑103 (MQ=255)
                                                                                                                                         aaCCCAGGGCAACTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCa                                                      <  1:2076206/103‑1 (MQ=255)
                                                                                                                                           cccAGGGCAACTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTATAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGc     >  2:768339/1‑150 (MQ=255)
                                                                                                                                            ccAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAggg      <  1:1577321/148‑1 (MQ=255)
                                                                                                                                            ccAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAggg      >  2:1577321/1‑148 (MQ=255)
                                                                                                                                            ccAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAg        >  1:1532561/1‑146 (MQ=255)
                                                                                                                                            ccAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAg        <  2:1532561/146‑1 (MQ=255)
                                                                                                                                            ccAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCg    >  2:2345953/1‑150 (MQ=255)
                                                                                                                                            ccAGGGCAACTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTgg                         >  1:1934032/1‑129 (MQ=255)
                                                                                                                                            ccAGGGCAACTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTgg                         <  2:1934032/129‑1 (MQ=255)
                                                                                                                                            ccAGGGCAACTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAggg      >  2:1542602/1‑148 (MQ=255)
                                                                                                                                            ccAGGGCAACTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAggg      <  1:1542602/148‑1 (MQ=255)
                                                                                                                                            ccAGGGCAACTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCg    >  2:1711156/1‑150 (MQ=255)
                                                                                                                                            ccAGGGCAACTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCg    >  1:1102111/1‑150 (MQ=255)
                                                                                                                                              aGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCGAt  >  2:394922/1‑150 (MQ=255)
                                                                                                                                                   aGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCa         <  1:1457991/138‑1 (MQ=255)
                                                                                                                                                   aGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCa         >  2:1457991/1‑138 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
GCAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCGAT  >  NC_002947/571472‑571763

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