Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,145,440 G→A 12.1% 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,4400GA12.1% 165.7 / 11.9 66intergenic (‑107/+34)hemO/PP_1006heme oxygenase/heme receptor
Reads supporting (aligned to +/- strand):  ref base G (32/26);  new base A (4/4);  total (36/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 = 6.31e-01

GTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAA  >  NC_002947/1145295‑1145588
                                                                                                                                                  |                                                                                                                                                    
gTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGt                                                                                                                                                     <  2:165723/147‑1 (MQ=255)
   cAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGccc                                                                                                                                                 <  5:160218/148‑1 (MQ=255)
   cAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTAGCCAGTACCGGCCTCTTCGCGGGCGTGccc                                                                                                                                                 <  1:36613/148‑1 (MQ=255)
      gcgcTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCa                                                                                                                                              >  3:148427/1‑148 (MQ=255)
         cTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAaga                                                                                                                                          <  7:209547/149‑1 (MQ=255)
            ggAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCGGGCACGCCCGCGAagag                                                                                                                                         <  3:182662/147‑1 (MQ=255)
             gAGCGCAGGGCGGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGGCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGcc                                                                                                                                      <  5:344583/149‑1 (MQ=255)
              aGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGGTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCg                                                                                                                                     <  1:332687/149‑1 (MQ=255)
                ctcAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCGGGCACGCCCGCGAAGAGGCCgg                                                                                                                                    <  8:192492/146‑1 (MQ=255)
                   aGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTa                                                                                                                                  >  2:252644/1‑147 (MQ=255)
                             ggtcgttcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGc                                                                                                                      <  4:148427/149‑1 (MQ=255)
                              gtcggtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGcc                                                                                                                     >  1:217401/1‑149 (MQ=255)
                               tcggtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCt                                                                                                                    >  6:204761/1‑149 (MQ=255)
                                    cATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCGGGCACGCCCGCGAAGAGGCCg                                                                                                                                     >  4:229962/1‑127 (MQ=255)
                                    cATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCGGGCACGCCCGCGAAGAGGCCg                                                                                                                                     <  3:229962/127‑1 (MQ=255)
                                            aGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCaaa                                                                                                       >  1:159671/1‑149 (MQ=255)
                                                    gAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACAtt                                                                                                >  3:76604/1‑148 (MQ=255)
                                                         agGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTTCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCACAGAGGCCGGTACCGAACTAACGCCTCAGCTATACCAACCCACCTTGAAAc                                                                                           >  6:37633/1‑148 (MQ=255)
                                                                 gaGACACAAGACGAAACACCCGGAGCCACCCCGAAAAAAAAGCCCCTACCTCCGATGCCAGAACCGGCCTCATCGCGCGCGTGCCCCCAAAGAGGCCGGTACAGACGTAACCCCTCAGATATCCCAAACCACATTGATACCGAAGTTgc                                                                                  <  6:104540/148‑1 (MQ=255)
                                                                  aGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTgcg                                                                                 >  8:363577/1‑149 (MQ=255)
                                                                      aCAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGcccc                                                                             >  4:56902/1‑149 (MQ=255)
                                                                          gACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGc                                                                                                                                       >  6:88776/1‑87 (MQ=255)
                                                                          gACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGc                                                                                                                                       <  5:88776/87‑1 (MQ=255)
                                                                           aCGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGGGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCAt                                                                        <  4:14593/149‑1 (MQ=255)
                                                                             gacaCAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGCCGTAACGCCTCAGATATCCCAACCCACATTGCTACCGAAGTTGCGCCCCGGCATgg                                                                      >  6:87275/4‑149 (MQ=255)
                                                                              aaaCAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGt                                                                     <  5:204761/149‑1 (MQ=255)
                                                                               aaCAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGATGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTc                                                                    <  1:107396/149‑1 (MQ=255)
                                                                                aCAGCGCGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCACATCGCCGGCGGCGCCCCCAAAACGGCCCGTCACACGGAACGCCGCTGCAATACCAAAACACAAAGATACCGGAGATGGGCCCCGGCAAGGt                                                                     >  2:193633/1‑147 (MQ=255)
                                                                                 cAGCGGGAGCGAGACAGTAAAAAAAGCCGGTCTTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAg                                                                  <  5:37633/149‑1 (MQ=255)
                                                                                  aGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAgg                                                                 <  2:111300/149‑1 (MQ=255)
                                                                                    cGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCg                                                               >  3:171902/1‑149 (MQ=255)
                                                                                          cGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGAt                                                         <  1:252644/149‑1 (MQ=255)
                                                                                           gagaCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATg                                                        >  2:35433/1‑149 (MQ=255)
                                                                                            agaCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTgc                                                                                  >  6:128585/1‑122 (MQ=255)
                                                                                            agaCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTgc                                                                                  <  5:128585/122‑1 (MQ=255)
                                                                                               cAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCTCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtg                                                     <  4:76604/148‑1 (MQ=255)
                                                                                                aGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtggc                                                   >  7:233078/1‑149 (MQ=255)
                                                                                                 gTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtg                                                     <  3:327829/146‑1 (MQ=255)
                                                                                                 gTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtg                                                     >  4:327829/1‑146 (MQ=255)
                                                                                                  tAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACAt                                                                                                 >  5:144725/1‑101 (MQ=255)
                                                                                                  tAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACAt                                                                                                 <  6:144725/101‑1 (MQ=255)
                                                                                                  tAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGTTCGATGTtggctg                                                 >  5:38822/1‑149 (MQ=255)
                                                                                                    aaaaaaGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtggctggc                                               >  4:115361/1‑149 (MQ=255)
                                                                                                     aaaaaGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTggctggct                                              >  6:150464/1‑149 (MQ=255)
                                                                                                       ataGCCGGTATTTGCGGTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCGCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGctggctgg                                            <  1:193633/147‑1 (MQ=255)
                                                                                                         aGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGt                                          <  3:56902/149‑1 (MQ=255)
                                                                                                             ggTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTcaca                                      <  7:12885/149‑1 (MQ=255)
                                                                                                                aTTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGc                                                                                                                                       >  5:271253/1‑49 (MQ=255)
                                                                                                                aTTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGc                                                                                                                                       <  6:271253/49‑1 (MQ=255)
                                                                                                                   tGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAgccgc                                 <  3:188961/148‑1 (MQ=255)
                                                                                                                          ccAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCt                         >  3:231120/1‑149 (MQ=255)
                                                                                                                           cAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTg                        >  1:49461/1‑149 (MQ=255)
                                                                                                                                 cGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcgtc                  >  5:391887/1‑149 (MQ=255)
                                                                                                                                 cGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcgtc                  >  2:77173/1‑149 (MQ=255)
                                                                                                                                 cGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGACCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCAGGCagtc                  >  8:56667/1‑149 (MQ=255)
                                                                                                                                 cGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATCGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcgtc                  >  5:179267/1‑149 (MQ=255)
                                                                                                                                   gCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcgtcgt                <  5:150464/149‑1 (MQ=255)
                                                                                                                                   gCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcgtcg                 >  6:373828/1‑148 (MQ=255)
                                                                                                                                   gCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGAGACCCGGCATGGTCAGGAGGTCGATGTAGGCTGGCAGGGTCACAGCAGCCTCTCCCTGGCcgtcg                 >  6:226920/1‑148 (MQ=255)
                                                                                                                                    ccTCTTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcgtcg                 >  2:197366/1‑147 (MQ=255)
                                                                                                                                    ccTCTTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTa               <  1:197366/149‑1 (MQ=255)
                                                                                                                                    ccTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcgtcg                 >  7:92227/1‑147 (MQ=255)
                                                                                                                                    ccTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTa               <  7:363577/149‑1 (MQ=255)
                                                                                                                                     ctctTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTCATACAGAAGTTGCGCCACGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACCGCCGCCTCACCCTGTCCGTCGTac              >  6:169714/1‑148 (MQ=255)
                                                                                                                                      tctTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGc             >  4:80898/1‑149 (MQ=255)
                                                                                                                                      tctACGAGGCCGTCCCCGCAAAGCGGCCGGTACACACGTCACGCCTCAGATATCCCAAACCACAAACATAACCAAGTTGCGCCACGGCATCCACACCAGCTCCATCATGGCAGGCTCCCTCACAGCCACCTCACACTGGCCGACGTac              >  4:113622/1‑147 (MQ=255)
                                                                                                                                                  gTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGaa  >  4:7605/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAA  >  NC_002947/1145295‑1145588

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