Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,564,575 C→T 14.4% intergenic (‑5/+90) hfq ← / ← miaA RNA‑binding protein Hfq/tRNA (adenosine(37)‑N6)‑dimethylallyltransferase MiaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,564,5750CT14.4% 231.4 / 26.4 97intergenic (‑5/+90)hfq/miaARNA‑binding protein Hfq/tRNA (adenosine(37)‑N6)‑dimethylallyltransferase MiaA
Reads supporting (aligned to +/- strand):  ref base C (42/40);  new base T (8/6);  total (50/47)
Fisher's exact test for biased strand distribution p-value = 7.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

TTCTTCAGCAATACCACGAACTGGTCGAACGATTCGATCGAGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGACTCAGCTCAATATGGAGATGGCCCCAAGGTATTTCAAGGTGCGGGACAGATTGTCGCAGGC  >  NC_002947/5564428‑5564724
                                                                                                                                                   |                                                                                                                                                     
ttcttcAGCAATACCACGAACTGGTCGAACGATTCGATCGAGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGc                                                                                                                                                     >  1:2096367/1‑150 (MQ=255)
      aGCAATACCACGAACTGGTCGAACGATTCGATCGAGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCt                                                                                                                                               >  1:72766/1‑150 (MQ=255)
        cAATACCACGAACTGGTCGAACGATTCGATCGAGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGt                                                                                                                                             <  1:1248300/150‑1 (MQ=255)
         aaTACCACGAACTGGTCGAACGATTCGATCGAGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTa                                                                                                                                            <  1:18581/150‑1 (MQ=255)
          aTACCACGAACTGGTCGAACGATTCGATCGAGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTaa                                                                                                                                           <  2:224181/150‑1 (MQ=255)
          aTACCACGAACTGGTCGAACGATTCGATCGAGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTaa                                                                                                                                           <  1:979712/150‑1 (MQ=255)
          aTACCACGAACTGGTCGAACGATTCGATCGAGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTaa                                                                                                                                           <  1:1519493/150‑1 (MQ=255)
           tACCACGAACTGGTCGAACGATTCGATCGAGCCCTACAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAg                                                                                                                                          <  2:1690713/150‑1 (MQ=255)
              cACGAACTGGTCGAACGATTCGATCGAGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTATCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCGCTCCTGTAAGGat                                                                                                                                       >  2:1296236/1‑150 (MQ=255)
               aCGAACTGGTCGAACGATTCGATCGAGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGata                                                                                                                                      >  1:1335307/1‑150 (MQ=255)
                  aaCTGGTCGAACGATTCGATCGAGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTCCGGataaat                                                                                                                                   <  2:2313727/150‑1 (MQ=255)
                  aaCTGGTCGAACGATTCGATCGAGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGataaat                                                                                                                                   <  1:2042445/150‑1 (MQ=255)
                    cTGGTCGAACGATTCGATCGAGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGataaataa                                                                                                                                 <  2:1335307/150‑1 (MQ=255)
                     tGGTCGAACGATTCGATCGAGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAaataat                                                                                                                                >  1:1858549/1‑150 (MQ=255)
                     tGGTCGAACGATTCGATCGAGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAaataat                                                                                                                                >  2:592628/1‑150 (MQ=255)
                     tGGTCGAACGATTCGATCGAGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCATGTACGGataaggaat                                                                                                                                >  2:1995477/1‑145 (MQ=255)
                       gTCGAACGATTCGATCGAGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAg                                                                                                                              <  1:579237/150‑1 (MQ=255)
                          gAACGATTCGATCGAGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGaaa                                                                                                                           >  1:1602133/1‑150 (MQ=255)
                          gAACGATTCGATCGAGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGaaa                                                                                                                           >  2:267029/1‑150 (MQ=255)
                           aaCGATTCGATCGAGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGAGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGaaaa                                                                                                                          <  1:16184/150‑1 (MQ=255)
                           aaCGATTCGATCGAGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGAGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGaaaa                                                                                                                          >  2:16184/1‑150 (MQ=255)
                                ttCGATCGAGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTTTTTTTTTTTTTTTTCCTGTAATGATAAATAATAGAAAATCAAt                                                                                                                     >  2:697901/1‑150 (MQ=255)
                                ttCGATCGAGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCTGCACTCCTGTAAGGATAAATAATAGAAAATCAAt                                                                                                                     >  1:1851936/1‑150 (MQ=255)
                                   gatcgaGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAATGAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATaaa                                                                                                                  >  2:1045128/1‑150 (MQ=255)
                                        aGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGa                                                                                                             <  2:944519/150‑1 (MQ=255)
                                        aGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGa                                                                                                             <  2:2066319/150‑1 (MQ=255)
                                         gCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCTGCACTCCTGt                                                                                                                                             <  2:2453275/117‑1 (MQ=255)
                                         gCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCTGCACTCCTGt                                                                                                                                             >  1:2453275/1‑117 (MQ=255)
                                         gCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGAt                                                                                                            <  2:2440298/150‑1 (MQ=255)
                                          cccTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATa                                                                                                           >  2:2187040/1‑150 (MQ=255)
                                                aGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTCGTGGGAGATGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCTGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAg                                                                                                     <  2:63788/150‑1 (MQ=255)
                                                 gTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTaa                                                                                                                                           <  1:684426/111‑1 (MQ=255)
                                                 gTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTaa                                                                                                                                           >  2:684426/1‑111 (MQ=255)
                                                 gTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAgg                                                                                                    <  2:1858549/150‑1 (MQ=255)
                                                       tACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCAGTTAAATTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGccc                                                                                              <  1:1995477/150‑1 (MQ=255)
                                                          cGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCGCtt                                                                                           >  1:84681/1‑150 (MQ=255)
                                                          cGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCtt                                                                                           >  1:742110/1‑150 (MQ=255)
                                                          cGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGCCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCtt                                                                                           <  2:1319773/150‑1 (MQ=255)
                                                            ttGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCg                                                                                         <  2:2096367/150‑1 (MQ=255)
                                                                ccAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAg                                                                                                     >  2:1967280/1‑134 (MQ=255)
                                                                ccAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAg                                                                                                     <  1:1967280/134‑1 (MQ=255)
                                                                 cAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAataata                                                                                                                               >  2:2258903/1‑107 (MQ=255)
                                                                 cAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAataata                                                                                                                               <  1:2258903/107‑1 (MQ=255)
                                                                  aGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCTGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGaa                                                                                   <  1:2096627/150‑1 (MQ=255)
                                                                  aGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAAGGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGaa                                                                                   <  2:72766/150‑1 (MQ=255)
                                                                    gTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCCTTTGACATATGCCGCACTCATGTAAGGATAAATAATAGaaaa                                                                                                                          <  1:771253/109‑1 (MQ=255)
                                                                    gTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCCTTTGACATATGCCGCACTCATGTAAGGATAAATAATAGaaaa                                                                                                                          >  2:771253/1‑109 (MQ=255)
                                                                      aaaTGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGa                                                                               <  2:1602133/150‑1 (MQ=255)
                                                                        aTGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCt                                                                                                                                               >  2:325561/1‑84 (MQ=255)
                                                                        aTGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCt                                                                                                                                               <  1:325561/84‑1 (MQ=255)
                                                                         tGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACgg                                                                            <  1:2187040/150‑1 (MQ=255)
                                                                         tGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACgg                                                                            <  1:267029/150‑1 (MQ=255)
                                                                          ggATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGc                                                                           >  2:549100/1‑150 (MQ=255)
                                                                           gATACAGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCTGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCa                                                                          >  2:2434759/1‑150 (MQ=255)
                                                                                 gggACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGAAAGACTAGACGGAAATcagc                                                                    <  2:84681/150‑1 (MQ=255)
                                                                                  ggACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATcagca                                                                   <  1:592628/150‑1 (MQ=255)
                                                                                    aCCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAgg                                                                 <  2:229564/150‑1 (MQ=255)
                                                                                       ttttcttttcTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGAc                                                              >  1:1281270/1‑150 (MQ=255)
                                                                                           cttttctCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCTGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCCATCAGCACCGACTCAg                                                          <  2:1851936/150‑1 (MQ=255)
                                                                                           cttttctCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCaa                                                                                      >  1:2453177/1‑122 (MQ=255)
                                                                                           cttttctCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCaa                                                                                      <  2:2453177/122‑1 (MQ=255)
                                                                                            ttttctCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGACTCAGc                                                         >  2:1426970/1‑150 (MQ=255)
                                                                                                ctcAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGACTCAGCTCaa                                                     >  1:771014/1‑150 (MQ=255)
                                                                                                 tcAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCTGCACTCCTGTAAGGATAaataat                                                                                                                                >  1:887629/1‑74 (MQ=255)
                                                                                                 tcAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCTGCACTCCTGTAAGGATAaataat                                                                                                                                <  2:887629/74‑1 (MQ=255)
                                                                                                 tcAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAg                                                                                                     <  2:1804473/101‑1 (MQ=255)
                                                                                                 tcAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAg                                                                                                     >  1:1804473/1‑101 (MQ=255)
                                                                                                 tcAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGACTCa                                                           >  2:1230553/1‑143 (MQ=255)
                                                                                                 tcAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGACTCa                                                           <  1:1230553/143‑1 (MQ=255)
                                                                                                  cAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGataaa                                                                                                                                    <  2:2078306/69‑1 (MQ=255)
                                                                                                  cAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGataaa                                                                                                                                    >  1:2078306/1‑69 (MQ=255)
                                                                                                   aaGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTc                                                                                          >  1:247240/1‑110 (MQ=255)
                                                                                                   aaGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTc                                                                                          <  2:247240/110‑1 (MQ=255)
                                                                                                            aaGTAAGGGCCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGa                                                                                                             >  1:1894761/1‑82 (MQ=255)
                                                                                                            aaGTAAGGGCCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGa                                                                                                             <  2:1894761/82‑1 (MQ=255)
                                                                                                              gTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGACTCAGCTCAATATGGAGATGGccc                                       >  1:1954794/1‑150 (MQ=255)
                                                                                                                   ggTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGACTCAGCTCAATATGGAGATg                                           >  1:916841/1‑141 (MQ=255)
                                                                                                                   ggTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGACTCAGCTCAATATGGAGATg                                           <  2:916841/141‑1 (MQ=255)
                                                                                                                     tCTTGTAGCGAATGCCCTTTTGACATATGCTGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGACTCAGCTCAATATGGAGATGGCCCCAAGGTa                                >  2:1725530/1‑150 (MQ=255)
                                                                                                                     tCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAATGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATcagcag                                                                  >  1:1495789/1‑116 (MQ=255)
                                                                                                                     tCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATcagcag                                                                  <  2:1495789/116‑1 (MQ=255)
                                                                                                                         gTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGACTCAGCTCAATATGGAGATGGCCCCAAGGTATTTc                            >  2:444633/1‑150 (MQ=255)
                                                                                                                          tAGCGAATGCCCTTTTGACATATGCTGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATcagcag                                                                  >  2:2243607/1‑111 (MQ=255)
                                                                                                                          tAGCGAATGCCCTTTTGACATATGCTGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATcagcag                                                                  <  1:2243607/111‑1 (MQ=255)
                                                                                                                          tAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGACTCAGCTAAATATGGAGAAGTCCCCAAGGTATTTCa                           >  1:525916/1‑150 (MQ=255)
                                                                                                                            tcGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGACTCAGCTCAATATGGAGATGGCCCCAAGGTATTTCAAg                         >  1:2179783/2‑150 (MQ=255)
                                                                                                                            gCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGCTAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGACTCAGCTCAATATGGAGATGGCCCCAAGGTATTTCAAg                         >  1:51322/1‑150 (MQ=255)
                                                                                                                               aaTGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGACTCAGCTCaa                                                     >  1:1611691/1‑119 (MQ=255)
                                                                                                                               aaTGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGACTCAGCTCaa                                                     <  2:1611691/119‑1 (MQ=255)
                                                                                                                                 tGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGACTCAGCTCAATATGGAGATGGCCCCAAGGTATTTCAAGGTGCg                    <  2:742110/150‑1 (MQ=255)
                                                                                                                                    ccTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGACTCAGCTCAATATGGAGATGGCCCCAAGGTATTTCAAGGTGCGGGa                 >  1:1915926/1‑150 (MQ=255)
                                                                                                                                      ttttGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGACTCAGCTCAATATGGAGATGGCCCCAAGGTATTTCAAGGTGCGGGACa               <  1:1426970/150‑1 (MQ=255)
                                                                                                                                           aCATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGACTCAGCTCAATATGGAGATGGCCCCAAGGTATTTCAAGGTGCGGGACAGATTg          >  2:1953260/1‑150 (MQ=255)
                                                                                                                                           aCATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGACTCAGCTCAATATGGAGATGGCCCCAAGGTATTTCAAGGTGCGGGACAGATTg          >  1:1647111/1‑150 (MQ=255)
                                                                                                                                             atatGCTGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGACTCAGCTCAATATGGAGATGGCCCCAAGGTATTTCAAGGTGCGGGACAGATTGTc        >  2:1697547/1‑150 (MQ=255)
                                                                                                                                             atatGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGACTCAGCTCAATATGGAGATGGCCCCAAGGTATTTCAAGGTGCGGGACGGATTGTc        >  1:642422/1‑150 (MQ=255)
                                                                                                                                              tatGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGACTCAGCTCAATATGGAGATGGCCCCAAGGTATTTCAAGGTGCGGGACAGATTGTCg       >  2:1923653/1‑150 (MQ=255)
                                                                                                                                                  ccGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGACTCAGCTCAATATGGAGATGGCCCCAAGGTATTTCAAGGTGCGGGACAGATTGTCGcagg   >  2:67990/1‑150 (MQ=255)
                                                                                                                                                   cGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGACTCAGCTCAATATGGAGATGGCCCCAAGGTATTTCAAGGTGCGGGACAGATTGTCGcaggc  <  2:771014/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTCTTCAGCAATACCACGAACTGGTCGAACGATTCGATCGAGCCCTGCAGTTTGATACCGTTGACCAGGTAAATGGATACCGGGACCTTTTCTTTTCTCAAGGTGTTCAAGTAAGGGTCTTGTAGCGAATGCCCTTTTGACATATGCCGCACTCCTGTAAGGATAAATAATAGAAAATCAATAAAATCGATAAATTAGGCCGCCCCTTCGCAAGAATAGACGGCAATCAGCAGGGACTCAGCTCAATATGGAGATGGCCCCAAGGTATTTCAAGGTGCGGGACAGATTGTCGCAGGC  >  NC_002947/5564428‑5564724

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