Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,145,446 G→A 9.3% intergenic (‑107/+34) hemO ← / ← PP_1006 heme oxygenase/heme receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,145,4460GA 9.3% 195.3 / 9.7 75intergenic (‑107/+34)hemO/PP_1006heme oxygenase/heme receptor
Reads supporting (aligned to +/- strand):  ref base G (41/27);  new base A (4/3);  total (45/30)
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.05e-01

CTGGTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAAT  >  NC_002947/1145298‑1145595
                                                                                                                                                     |                                                                                                                                                     
cTGGTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGt                                                                                                                                                      <  2:1564753/150‑1 (MQ=255)
  ggTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGc                                                                                                                                                    >  2:1885769/1‑150 (MQ=255)
            cTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCGGGCACGCCCGCGAagag                                                                                                                                          <  1:1576724/150‑1 (MQ=255)
             tGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAaga                                                                                                                                           >  2:1676828/1‑148 (MQ=255)
             tGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAaga                                                                                                                                           <  1:1676828/148‑1 (MQ=255)
                gAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGCCAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCg                                                                                                                                      >  1:596347/1‑150 (MQ=255)
                           tGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTa                                                                                                                           >  2:1603353/1‑150 (MQ=255)
                                  tcggtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTc                                                                                                                    >  1:1460828/1‑150 (MQ=255)
                                     gtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGa                                                                                                                 >  2:368061/1‑150 (MQ=255)
                                      tcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAACAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAAGGCCTCAgct                                                                                                                >  1:316785/1‑148 (MQ=255)
                                      tcATGGGGGAGTCCTTGGAAAGACAGCGTCTTATTACAGTACGAACCACCGGGACCGGTGCAGC‑ATTATCCAATTTCATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGat                                                                                                                <  1:975286/150‑1 (MQ=255)
                                       cATGGGGGAGACCTTGGAAAGAGGGGGGGTAGAAACAAGACGAAACAGCGGGAGCGAAAAAGA‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGata                                                                                                               <  2:249718/150‑1 (MQ=255)
                                        aTGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGTGATCGCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTCCAGACCTAACGCCTCAGatat                                                                                                              >  1:594001/1‑150 (MQ=255)
                                          gggggAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATcc                                                                                                            >  1:1832355/1‑150 (MQ=255)
                                           ggggAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATccc                                                                                                           <  1:994086/150‑1 (MQ=255)
                                              gAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCaaa                                                                                                        <  1:1885769/150‑1 (MQ=255)
                                              gAGTCCTGGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCaaa                                                                                                        <  1:1174744/150‑1 (MQ=255)
                                                  ccTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAAGCCGGTATTCGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACca                                                                                                     >  1:549317/1‑150 (MQ=255)
                                                   cTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACCGGCGTAACGCCTCAGATATCCCAAACcaca                                                                                                   >  1:1735216/1‑150 (MQ=255)
                                                   cTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACcaca                                                                                                   >  1:1595114/1‑150 (MQ=255)
                                                    ttGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACAt                                                                                                  <  2:1595114/150‑1 (MQ=255)
                                                    ttGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACAt                                                                                                  >  1:1303572/1‑150 (MQ=255)
                                                    ttGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGGTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACTCCTCGGTTATCCCAAACCACAt                                                                                                  >  2:1842576/1‑150 (MQ=255)
                                                       gAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCGGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGa                                                                                               >  1:114844/1‑150 (MQ=255)
                                                              ggCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTTCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGaa                                                                                        <  1:1228533/150‑1 (MQ=255)
                                                                   cTAGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTc                                                                                                                    >  1:144662/1‑117 (MQ=255)
                                                                   cTAGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTc                                                                                                                    <  2:144662/117‑1 (MQ=255)
                                                                     aGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGaa                                                                                        >  2:1417449/1‑143 (MQ=255)
                                                                     aGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGaa                                                                                        <  1:1417449/143‑1 (MQ=255)
                                                                      gATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGcc                                                                                >  2:523316/1‑150 (MQ=255)
                                                                        tACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGcccc                                                                              >  2:1592489/1‑150 (MQ=255)
                                                                          cAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTgc                                                                                   <  1:6053/143‑1 (MQ=255)
                                                                          cAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTgc                                                                                   >  2:6053/1‑143 (MQ=255)
                                                                             gACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGACGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCAt                                                                         >  1:674596/1‑150 (MQ=255)
                                                                              aCGAAACAGCAGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCaaa                                                                                                        >  2:1898911/1‑118 (MQ=255)
                                                                              aCGAAACAGCAGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCaaa                                                                                                        <  1:1898911/118‑1 (MQ=255)
                                                                                gAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGt                                                                      <  2:1832355/150‑1 (MQ=255)
                                                                                  aaCAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCa                                                                    >  1:1778512/1‑150 (MQ=255)
                                                                                      gCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGGCAGGCg                                                                >  1:1016554/1‑150 (MQ=255)
                                                                                            gCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGGACAGACGTAACGCCTCAGATATGCCAAACCACAGGGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGAt                                                          >  1:890404/1‑150 (MQ=255)
                                                                                              gagaCAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGt                                                        >  1:1407585/1‑150 (MQ=255)
                                                                                              gagaCAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGt                                                        >  2:984480/1‑150 (MQ=255)
                                                                                               agaCAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTc                                                                                                                    >  2:846647/1‑89 (MQ=255)
                                                                                               agaCAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTc                                                                                                                    <  1:846647/89‑1 (MQ=255)
                                                                                                gaCAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtg                                                      >  1:587058/1‑150 (MQ=255)
                                                                                                    gT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtggctg                                                  >  2:1043748/1‑150 (MQ=255)
                                                                                                        aaaaaaGCCGGTATTTGCGTTGCCAGTCCCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTggctggct                                               <  2:588690/150‑1 (MQ=255)
                                                                                                          aaaaGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGctggcngg                                             <  1:1043748/150‑1 (MQ=255)
                                                                                                            aaGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCTGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGGTCTTGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGt                                           <  1:1789135/150‑1 (MQ=255)
                                                                                                            aaGCCGGTATTTGCGGTGCCAGTCCCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGt                                           <  1:1710563/150‑1 (MQ=255)
                                                                                                               ccGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTcac                                        >  1:464518/1‑150 (MQ=255)
                                                                                                                cGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTcaca                                       >  2:198453/1‑150 (MQ=255)
                                                                                                                  gTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAgc                                     <  2:1559761/150‑1 (MQ=255)
                                                                                                                     tttGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAgccgc                                  <  1:1682391/150‑1 (MQ=255)
                                                                                                                          gTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCAc                             <  2:594001/150‑1 (MQ=255)
                                                                                                                           ttGCCGGTACCGGCCTCGTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCAcc                            >  2:1269724/1‑150 (MQ=255)
                                                                                                                           ttGCCAGTACCGGCCCCTCCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCAcc                            <  2:114844/150‑1 (MQ=255)
                                                                                                                              ccAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGt                                                              <  1:744443/113‑1 (MQ=255)
                                                                                                                              ccAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGt                                                              >  2:744443/1‑113 (MQ=255)
                                                                                                                              ccAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTg                         <  1:377154/150‑1 (MQ=255)
                                                                                                                                     cGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcgtcg                  >  1:1616618/1‑150 (MQ=255)
                                                                                                                                     cGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCCCAGATATCCCAAACCACCTTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcgtcg                  >  1:1917060/1‑150 (MQ=255)
                                                                                                                                       gCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGa                                                                                               >  2:1710351/1‑71 (MQ=255)
                                                                                                                                       gCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGa                                                                                               <  1:1710351/71‑1 (MQ=255)
                                                                                                                                        ccTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACcac                                                                                                    <  1:1396954/65‑1 (MQ=255)
                                                                                                                                        ccTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACcac                                                                                                    >  2:1396954/1‑65 (MQ=255)
                                                                                                                                         ctctTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAg                                                                                       >  1:1242126/1‑77 (MQ=21)
                                                                                                                                         ctctTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAg                                                                                       <  2:1242126/77‑1 (MQ=25)
                                                                                                                                         ctctTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGc              >  2:742417/1‑150 (MQ=255)
                                                                                                                                          tctTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGgctggctg                                              >  2:707371/1‑117 (MQ=38)
                                                                                                                                          tctTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGgctggctg                                              <  1:707371/117‑1 (MQ=38)
                                                                                                                                          tctTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCt             >  2:1771551/1‑150 (MQ=255)
                                                                                                                                          tctTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCTAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGATCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCt             >  2:1000337/1‑150 (MQ=255)
                                                                                                                                            ttCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCTTATATCCCATACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTgc           >  1:1281080/1‑150 (MQ=255)
                                                                                                                                            ttCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTgc           >  2:559538/1‑150 (MQ=255)
                                                                                                                                                     gTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCCAACTGGACGTCGTCGCCGCGCACCGact  <  2:1460828/150‑3 (MQ=255)
                                                                                                                                                     gTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCACCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAAt  >  1:156372/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGGTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGT‑AAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAAT  >  NC_002947/1145298‑1145595

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