Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 997,751 G→C 5.5% intergenic (‑156/+76) PP_0860 ← / ← PP_0861 sulfite reductase flavoprotein subunit/outer membrane ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947997,7510GC 5.5% 239.0 / 3.8 73intergenic (‑156/+76)PP_0860/PP_0861sulfite reductase flavoprotein subunit/outer membrane ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (32/37);  new base C (2/2);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

CCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCGCCTGGTTGGCGAAGTGGGTAGAGAAG  >  NC_002947/997607‑997896
                                                                                                                                                |                                                                                                                                                 
ccTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAAc                                                                                                                                               >  4:78901/1‑149 (MQ=255)
 cTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACg                                                                                                                                              >  1:93905/1‑149 (MQ=255)
    gAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCa                                                                                                                                            <  2:93905/148‑1 (MQ=255)
      gTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAg                                                                                                                                           >  1:249660/1‑147 (MQ=255)
       ttGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGa                                                                                                                                         <  4:365067/148‑1 (MQ=255)
            gcCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCGTTGCCGGTGCAGGAGCGGCt                                                                                                                                   >  4:316303/1‑149 (MQ=38)
              cgcgACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTc                                                                                                                                 <  6:299222/149‑1 (MQ=255)
               gcgACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCGTTGCCGGTGCAGGAGCGGCttta                                                                                                                                <  3:316303/149‑3 (MQ=37)
                  aCCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCGCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGcc                                                                                                                             <  7:235185/149‑1 (MQ=255)
                    ccGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCcgc                                                                                                                           <  2:249660/149‑1 (MQ=255)
                        gcggcgTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGaaa                                                                                                                       <  5:218808/149‑1 (MQ=255)
                            cgTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGc                                                                                                                   >  6:168584/1‑149 (MQ=255)
                              tGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCg                                                                                                                 <  4:270940/149‑1 (MQ=255)
                                cATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGc                                                                                                               >  4:24068/1‑149 (MQ=255)
                                 aTTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGca                                                                                                              >  8:261971/1‑149 (MQ=255)
                                     tGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAga                                                                                                          <  6:108314/149‑2 (MQ=255)
                                     tGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAg                                                                                                           <  3:24068/148‑1 (MQ=255)
                                       ccTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGgca                                                                                                        >  3:40326/1‑148 (MQ=255)
                                        cTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGgcg                                                                                                        >  2:86780/1‑148 (MQ=255)
                                        cTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGgcg                                                                                                        >  6:265039/1‑148 (MQ=255)
                                           gATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCGGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAg                                                                                                     <  7:94111/148‑1 (MQ=255)
                                             tACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGccc                                                                                                  <  7:116890/149‑1 (MQ=255)
                                              acacAAAAGCCCCGCTCACGGCAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCt                                                                                                 <  3:312862/149‑1 (MQ=255)
                                                acaAAAGCCCCGCTCTGGGGTATGATCGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCGTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTc                                                                                               <  6:312214/149‑1 (MQ=255)
                                                acaAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCa                                                                                                                                <  3:126556/116‑1 (MQ=255)
                                                acaAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCa                                                                                                                                >  4:126556/1‑116 (MQ=255)
                                                    aaGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGgcg                                                                                                        <  1:233961/136‑1 (MQ=255)
                                                    aaGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGgcg                                                                                                        >  2:233961/1‑136 (MQ=255)
                                                      gCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGa                                                                                         <  3:78901/149‑1 (MQ=255)
                                                      gCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCACAGCCCTTCAACAGa                                                                                         <  4:40326/149‑1 (MQ=255)
                                                        cccGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATc                                                                                       <  7:159236/149‑1 (MQ=255)
                                                         ccGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATc                                                                                       >  3:174228/1‑148 (MQ=255)
                                                          cGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCCCCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATcaa                                                                                     <  5:265039/149‑1 (MQ=255)
                                                               aCGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCACCGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATcaacccaa                                                                                >  8:161446/1‑149 (MQ=255)
                                                                 gggAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACCGGCGCAGCCCTTCAACAGATcaacccaacc                                                                              >  7:163611/1‑149 (MQ=255)
                                                                 gggAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATcaacccaacc                                                                              >  4:362437/1‑149 (MQ=255)
                                                                  ggAATGAAGGGGGCCGGGAGGGCTGCTTTTTGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCTACGCAGTAGCTGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCaacccaaccc                                                                             <  2:88053/149‑1 (MQ=255)
                                                                      tGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGa                                                                          >  5:245845/1‑148 (MQ=255)
                                                                                 gggAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGACGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTg                                                               >  3:91984/1‑148 (MQ=255)
                                                                                    aGGGCTGCTTTCCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCGAACCCTGACCTCAGAAGTGCACg                                                           >  8:238617/1‑149 (MQ=255)
                                                                                          gctTTTCGTGTTGGCTGTGCGGGCCTCTGCCCGGCTACAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCtggtg                                                     <  8:350258/149‑1 (MQ=255)
                                                                                           ctTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCtggtgg                                                    <  7:261971/149‑1 (MQ=255)
                                                                                              ttCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTggtggt                                                   <  3:310469/147‑1 (MQ=255)
                                                                                                   gttgtCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGt                                            >  8:236606/1‑149 (MQ=255)
                                                                                                    ttgtCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGa                                                                                                                                         <  7:2471/55‑1 (MQ=255)
                                                                                                    ttgtCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGa                                                                                                                                         >  8:2471/1‑55 (MQ=255)
                                                                                                    ttgtCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTg                                           >  2:201896/1‑149 (MQ=255)
                                                                                                      gtCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGcc                                         >  1:82588/1‑149 (MQ=255)
                                                                                                           tgCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGt                                       <  7:236606/146‑1 (MQ=255)
                                                                                                                  ctctTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGgcc                             <  3:362437/149‑1 (MQ=255)
                                                                                                                    tttCGCGGCTGAAGCCGCTCCTGCGTTGCCGGTGCAGGAGCGGCTTTAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGgccgc                           >  5:284676/2‑149 (MQ=25)
                                                                                                                       cgcgGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGATATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGgccgcc                          <  1:86780/147‑1 (MQ=255)
                                                                                                                         cgGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCGCCTGGt                      <  2:359746/149‑1 (MQ=255)
                                                                                                                         cgGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCGCCTGGt                      <  5:168584/149‑1 (MQ=255)
                                                                                                                           gCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCGCCTGGTTg                    <  2:82588/149‑1 (MQ=255)
                                                                                                                             tAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCGCCTGGTTGGc                  >  1:226533/1‑149 (MQ=255)
                                                                                                                              aaaGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCGGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCGCCTGGTTGGc                  >  2:347717/1‑148 (MQ=255)
                                                                                                                               aaGCCGCTCCTGCGTTGCCGGTGCAGGAGCGGCTTTAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCGCCTGGTTGGc                  <  6:284676/147‑1 (MQ=37)
                                                                                                                                 gCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCGCCTGGTTGGCGAAg              >  6:2117/1‑149 (MQ=255)
                                                                                                                                  ccGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCGCCTGGTTGGCGAAGt             <  5:141983/149‑1 (MQ=255)
                                                                                                                                    gCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCaa                                                <  7:106787/112‑1 (MQ=255)
                                                                                                                                    gCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCaa                                                >  8:106787/1‑112 (MQ=255)
                                                                                                                                    gCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCGCCTg                        <  8:352063/136‑1 (MQ=255)
                                                                                                                                    gCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCGCCTg                        >  7:352063/1‑136 (MQ=255)
                                                                                                                                    gCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCGCCTGGTTGGCGAAGTgg           <  6:245845/149‑1 (MQ=255)
                                                                                                                                    gCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCGCCTGGTTGGCGAAGTgg           <  5:2117/149‑1 (MQ=255)
                                                                                                                                    gCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCGCCTGGGTGGCGAAGTgg           >  7:342842/1‑149 (MQ=255)
                                                                                                                                       ccTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCGCCTGGTTGGCGAAGTggg          >  5:106960/1‑147 (MQ=255)
                                                                                                                                         tGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCGCCTGGTTGGCGAAGTGCGTa        >  6:12176/1‑147 (MQ=255)
                                                                                                                                          gCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCGCCTGGTTGGCGAAGTGGGTa        >  5:39436/1‑146 (MQ=255)
                                                                                                                                             ccGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCGCCTGGTTGGCGAAGTGGGTAGAGAAg  <  8:163611/149‑1 (MQ=255)
                                                                                                                                              cGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCtg                                                        <  6:292544/94‑1 (MQ=255)
                                                                                                                                              cGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCtg                                                        >  5:292544/1‑94 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCGCCTGGTTGGCGAAGTGGGTAGAGAAG  >  NC_002947/997607‑997896

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