Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 997,736 A→G 5.3% intergenic (‑138/+94) 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,7360AG 5.3% 341.9 / 5.5 113intergenic (‑138/+94)PP_0860/PP_0861sulfite reductase flavoprotein subunit/outer membrane ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base A (60/45);  new base G (3/3);  total (65/48)
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.73e-01

Only 100 of 114 total aligned reads displayed.

TCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCGCCTGGTTGGCG  >  NC_002947/997589‑997884
                                                                                                                                                     |                                                                                                                                                    
tcttctTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑aaa                                                                                                                                                    <  1:584777/150‑1 (MQ=255)
  ttcttcACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGc                                                                                                                                                  <  2:291496/150‑1 (MQ=255)
     ttcACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGCCTGTGCCGGCCCTTTCGCGGCT‑‑GAAGCCGc                                                                                                                                               >  2:2195199/1‑150 (MQ=255)
       caccacGTTTCTCACCTTCAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGCGGGCTGCTTTTCGTGTTGTCTGTTCGGGCCCCTTCGCGGCG‑‑AAcgccgcca                                                                                                                                             >  1:1070171/1‑148 (MQ=255)
          cacGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTg                                                                                                                                          >  2:466678/1‑150 (MQ=255)
           acGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGc                                                                                                                                         >  2:1519894/1‑150 (MQ=255)
             gTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCAc                                                                                                                                       >  1:839806/1‑150 (MQ=255)
              tttCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGCATGAACGGGGCCGGGAGGGATGCTTTTGGTGTTGTTTGTGCGGGCCTCGTCGCGGCT‑‑AAAGCCGCTCCTGcccc                                                                                                                                      >  1:1904256/1‑150 (MQ=255)
                    aCCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAAc                                                                                                                                >  2:1322954/1‑150 (MQ=255)
                     ccTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACg                                                                                                                               >  1:1037516/1‑150 (MQ=255)
                      cTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑CAAGCCGCTCCTGCACCGGCAACGc                                                                                                                              >  1:2195185/1‑150 (MQ=255)
                       ttGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCa                                                                                                                             <  1:1732928/150‑1 (MQ=255)
                          aGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGa                                                                                                                          <  2:1268376/150‑1 (MQ=255)
                              gctgcCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCgg                                                                                                                      <  1:2148301/150‑1 (MQ=255)
                                  ccGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTc                                                                                                                  <  1:466678/150‑1 (MQ=255)
                                   cgcgACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCa                                                                                                                 >  1:1137157/1‑150 (MQ=255)
                                    gcgACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAg                                                                                                                >  2:1690723/1‑150 (MQ=255)
                                       aCCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCcg                                                                                                             >  2:2319363/1‑150 (MQ=255)
                                        cccGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCcgc                                                                                                            <  2:2212289/150‑1 (MQ=255)
                                         ccGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCcgcg                                                                                                           >  1:2306227/1‑150 (MQ=255)
                                            cgcgGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCCCCTGCCCCGGCAACGCAGGAGCGGCTTCAGCCGCGaaa                                                                                                        <  2:2021058/150‑1 (MQ=255)
                                              cggcgTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑GAAGCCGCTCCTGCGTTGCCGGTGCAGGAGCGGCTTTAGCCGCGaagag                                                                                                      <  1:2017046/150‑4 (MQ=21)
                                                  gTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCg                                                                                                  <  1:1951162/150‑1 (MQ=255)
                                                  gTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCCGGGCGGGCCTCTTCGCGGCTCCCGAACCGCGCCT‑‑ACCGGCATGGCAGGAGCGGCTTCAGCCGCGAAAGGGCCg                                                                                                  >  1:2270617/1‑150 (MQ=255)
                                                       ttATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGcaca                                                                                             >  1:2250742/1‑150 (MQ=255)
                                                            ccTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGgcgc                                                                                        >  2:826187/1‑150 (MQ=255)
                                                              ttGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑GAAGCCGCTc                                                                                                                                             >  1:612714/1‑95 (MQ=255)
                                                              ttGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑GAAGCCGCTc                                                                                                                                             <  2:612714/95‑1 (MQ=255)
                                                                  tACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCt                                                                                  >  2:2269431/1‑150 (MQ=255)
                                                                   acacAAAAGCCCCGCGCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCtt                                                                                 <  2:2195185/150‑1 (MQ=255)
                                                                    cacaAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTc                                                                                <  1:2214865/150‑1 (MQ=255)
                                                                     acaAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTcc                                                                               >  2:1538854/1‑149 (MQ=255)
                                                                     acaAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCa                                                                               <  2:1957538/150‑1 (MQ=255)
                                                                     acaAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCa                                                                               >  1:1709177/1‑150 (MQ=255)
                                                                     acaAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCa                                                                               >  1:743229/1‑150 (MQ=255)
                                                                      caAAAGCCCCGCTCACGGGACTGACCGGCCACGGCGGGTGTCCGGTGCGGGTTTTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCaa                                                                              <  1:2269805/150‑1 (MQ=255)
                                                                      caAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCaa                                                                              <  1:826187/150‑1 (MQ=255)
                                                                      caAAAGCCCCGATCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCaa                                                                              >  2:336746/1‑150 (MQ=255)
                                                                           gCCCCGCTCACGGGAATGAACGGGTCCGGGAGGGCTTCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAActgct                                                                         <  1:1248448/150‑5 (MQ=255)
                                                                           cccccGCTCACGGGAATGATCCGGTCTTGGGTGGTTTCTTTTCGTTTTTTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCaccgatc                                                                         <  1:1248421/149‑7 (MQ=255)
                                                                            ccccGCTCACGGGAATGAACGGGGCCGGGGGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATc                                                                        <  1:719572/150‑1 (MQ=255)
                                                                            ccccGCGCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATc                                                                        <  1:926319/150‑1 (MQ=255)
                                                                              ccGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCTGGCACAGGCGCAGCCCTTCAACAGATcaa                                                                      >  2:817387/1‑150 (MQ=255)
                                                                              ccGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATcaa                                                                      <  1:2023026/150‑1 (MQ=255)
                                                                              ccGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTGGGCTTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCAGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATcaa                                                                      >  1:1354294/1‑150 (MQ=255)
                                                                                 atCAGGGGCATAAACGAGGCCAGGAGGGCCGCCGTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATcaaccc                                                                   >  2:273557/2‑150 (MQ=255)
                                                                                  tCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATcaaccca                                                                  >  1:379660/1‑150 (MQ=255)
                                                                                    aCGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATcaacccaac                                                                >  2:588931/1‑150 (MQ=255)
                                                                                    aCGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATcaacccaac                                                                <  2:85645/150‑1 (MQ=255)
                                                                                     cGGGAGTGAACGGGGCCGGGAGGCCTGCTTGTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATcaacccaacc                                                               >  2:1256958/1‑150 (MQ=255)
                                                                                     cGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATcaacccaacc                                                               <  2:1924850/150‑1 (MQ=255)
                                                                                      gggAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCaacccaaccc                                                              >  1:623260/1‑150 (MQ=255)
                                                                                        gAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCaa                                                                                                                                 >  1:1779691/1‑81 (MQ=255)
                                                                                          aTGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCTGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGAc                                                          <  1:817387/150‑1 (MQ=255)
                                                                                                  cgCCGGACGGGTTGCTCATCGTGTGCTCTGCTCGGGCCTCTTCGCGGCG‑‑AAAGCCGCTCTGTTACCGGCAACGCAGGAGCGGCTTCAGCCGCGACAGGACCGGCACAGGCGCAGCCATTCAACAGACCAACCCCACCCTGACATCAGAAg                                                  >  2:26112/2‑150 (MQ=255)
                                                                                                   gCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGt                                                 <  2:2112564/150‑1 (MQ=255)
                                                                                                   gCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACACCAGAAGt                                                 <  2:623260/150‑1 (MQ=255)
                                                                                                    ccGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTg                                                >  2:2080379/1‑150 (MQ=255)
                                                                                                       ggAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCAc                                             >  2:618409/1‑150 (MQ=255)
                                                                                                          gggCTGTTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCt                                          >  2:358167/1‑150 (MQ=255)
                                                                                                           ggCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCa                                                                               >  1:1739159/1‑112 (MQ=255)
                                                                                                           ggCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCa                                                                               <  2:1739159/112‑1 (MQ=255)
                                                                                                             ctgctTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCtggt                                       >  1:421102/1‑150 (MQ=255)
                                                                                                             ctgctTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACAGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCtggt                                       >  2:147532/1‑150 (MQ=255)
                                                                                                              tgctTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGa                                                           >  2:607084/1‑129 (MQ=255)
                                                                                                              tgctTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGa                                                           <  1:607084/129‑1 (MQ=255)
                                                                                                              tgctTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCtggtg                                      >  1:1620362/1‑150 (MQ=255)
                                                                                                                ctTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTggtggt                                    >  1:1611201/1‑150 (MQ=255)
                                                                                                                ctTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACACAACCCTGACATCAGAAGTGCACGCTggtggt                                    <  2:743229/150‑1 (MQ=255)
                                                                                                                 ttttCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCgg                                                                                                                      <  2:2217682/67‑1 (MQ=255)
                                                                                                                 ttttCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCgg                                                                                                                      >  1:2217682/1‑67 (MQ=255)
                                                                                                                 ttttCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTc                                   <  1:1538854/150‑1 (MQ=255)
                                                                                                                     cGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACa                               <  2:403839/150‑1 (MQ=255)
                                                                                                                       tgttgtCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGt                             >  1:894040/1‑150 (MQ=255)
                                                                                                                        gttgtCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCaa                                                                              >  1:1101295/1‑100 (MQ=255)
                                                                                                                        gttgtCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCaa                                                                              <  2:1101295/100‑1 (MQ=255)
                                                                                                                         ttgtCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCg                                                                                                                       >  2:185102/1‑58 (MQ=255)
                                                                                                                         ttgtCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCg                                                                                                                       <  1:185102/58‑1 (MQ=255)
                                                                                                                            tCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGaaa                                                                                                        >  2:1730241/1‑70 (MQ=255)
                                                                                                                            tCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGaaa                                                                                                        <  1:1730241/70‑1 (MQ=255)
                                                                                                                            tCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGt                        >  2:369078/1‑150 (MQ=255)
                                                                                                                            ccTGTGCCGGCCCTTTCGCGGCT‑‑GAAGCCGCTCCTGCACCGGCAACGCAGGAGc                                                                                                                        >  2:1225384/2‑54 (MQ=1)
                                                                                                                            ccTGTGCCGGCCCTTTCGCGGCT‑‑GAAGCCGCTCCTGCACCGGCAACGCAGGAGc                                                                                                                        <  1:1225384/53‑1 (MQ=1)
                                                                                                                                   gggCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCg                 >  2:843351/1‑150 (MQ=255)
                                                                                                                                   gggCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCg                 >  1:1225790/1‑150 (MQ=255)
                                                                                                                                     gccctttCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGgc               >  1:1383758/6‑150 (MQ=255)
                                                                                                                                     gCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGgc               >  1:1482557/1‑150 (MQ=255)
                                                                                                                                       ctctTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGgccg             <  1:1391863/150‑1 (MQ=255)
                                                                                                                                       ctctTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGgccg             >  2:1571918/1‑150 (MQ=255)
                                                                                                                                        tctTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGgccgc            >  1:74686/1‑150 (MQ=255)
                                                                                                                                        tctTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGgccgc            <  2:2270617/150‑1 (MQ=255)
                                                                                                                                          ttCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTc                                                                                                                  <  2:922244/46‑1 (MQ=255)
                                                                                                                                          ttCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTc                                                                                                                  >  1:922244/1‑46 (MQ=255)
                                                                                                                                             gcgGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCt                                                                                  >  1:2260068/1‑75 (MQ=255)
                                                                                                                                             gcgGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCt                                                                                  <  2:2260068/75‑1 (MQ=255)
                                                                                                                                             gcgGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCGCCTGGt       <  1:2080379/150‑1 (MQ=255)
                                                                                                                                                 cT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCGCCTGGTTGGc   >  1:663198/1‑150 (MQ=255)
                                                                                                                                                  t‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGaa                                                   <  1:2071111/101‑1 (MQ=255)
                                                                                                                                                  t‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGaa                                                   >  2:2071111/1‑101 (MQ=255)
                                                                                                                                                  t‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCGCCTGGTTGGCg  >  2:767497/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCT‑‑AAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCGCCTGGTTGGCG  >  NC_002947/997589‑997884

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: