Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,079,126 T→A 14.7% intergenic (‑371/+26) mgtE ← / ← PP_t56 magnesium transporter/tRNA‑Arg

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,079,1260TA14.7% 164.7 / 19.0 68intergenic (‑371/+26)mgtE/PP_t56magnesium transporter/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base T (22/36);  new base A (5/5);  total (27/41)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

ATCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGC  >  NC_002947/5078982‑5079274
                                                                                                                                                  |                                                                                                                                                    
aTCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTAtt                                                                                                                                                   >  4:19321/1‑148 (MQ=255)
 tCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTAtt                                                                                                                                                   <  5:346071/147‑1 (MQ=255)
    gTGCGACGATTGTTGTGGATTGGAAAGCACATTTTATTTTTGAGGCCGCTAATCGGCTCCAAAATACACCGCAATTAATCCACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAg                                                                                                                                                <  1:107643/149‑1 (MQ=255)
       cgacgaTTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGc                                                                                                                                           >  2:169440/1‑149 (MQ=255)
        gacgaTTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCgg                                                                                                                                             >  3:192981/1‑146 (MQ=255)
        gacgaTTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCgg                                                                                                                                             <  4:192981/146‑1 (MQ=255)
        gacgaTTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCGGGCt                                                                                                                                          <  2:245510/149‑1 (MQ=255)
           gaTTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCATCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCtttt                                                                                                                                       <  5:335349/149‑1 (MQ=255)
           gaTTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCTAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCtttt                                                                                                                                       <  3:195745/149‑1 (MQ=255)
           gaTTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCtttt                                                                                                                                       >  3:782/1‑149 (MQ=255)
           gaTTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCtttt                                                                                                                                       <  6:315228/149‑1 (MQ=255)
           gaTTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCGGGCtttt                                                                                                                                       >  2:208097/1‑149 (MQ=255)
              tgttgtGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGa                                                                                                                                    >  5:325221/1‑149 (MQ=255)
               gttgtGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATACAGCGGGCTTTTTgt                                                                                                                                    <  1:208097/148‑2 (MQ=255)
                ttgtGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGAt                                                                                                                                   <  4:188653/148‑1 (MQ=255)
                ttgtGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATg                                                                                                                                  >  6:158311/1‑149 (MQ=255)
                ttgtGGATGGGAAAGCACATTTGATTTTTGAGGCCGCGCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGAt                                                                                                                                   <  3:106209/148‑1 (MQ=255)
                   tgGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATggg                                                                                                                                <  4:782/148‑1 (MQ=255)
                     gATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTa                                                                                                                             >  5:158374/1‑149 (MQ=255)
                     gATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGACTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATggg                                                                                                                                <  2:14299/146‑1 (MQ=255)
                     gATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGACTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATggg                                                                                                                                >  1:14299/1‑146 (MQ=255)
                       tGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATg                                                                                                                           <  6:333039/149‑1 (MQ=255)
                         ggAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATgg                                                                                                                          <  7:119843/148‑1 (MQ=255)
                             aGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGAc                                                                                                                     >  4:18501/1‑149 (MQ=255)
                               cacaTTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTc                                                                                                                   <  3:221689/149‑1 (MQ=255)
                                   tttGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATgg                                                                                                                                 <  7:74790/131‑1 (MQ=255)
                                   tttGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATgg                                                                                                                                 >  8:74790/1‑131 (MQ=255)
                                   tttGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTcagca                                                                                                               <  6:303246/149‑1 (MQ=255)
                                      gATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAgg                                                                                                             <  6:325221/148‑1 (MQ=255)
                                          tttggagCCGCTGAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTc                                                                                                         >  3:289618/7‑148 (MQ=255)
                                          tttGAGTCCTCTCAGCGTCTCCAAACTACACCGCAATCATT‑‑ACATAAGAGCGTATATCGACAGGCTCAAGAATTAATTTGAGCACAAACAAAAAGCCCTCTCTATTGAGCGGGCTTTTTGATGGTGTATGGCGGACTCAGCAGGATTc                                                                                                         <  5:98955/148‑1 (MQ=255)
                                              aGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAAc                                                                                                     <  1:169440/148‑1 (MQ=255)
                                                gCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGTTCGaatcctg                                                                                                 <  8:69712/149‑6 (MQ=255)
                                                gCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGTTCGaatcctg                                                                                                 >  7:69712/1‑144 (MQ=255)
                                                   gCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGa                                                                                               <  3:19321/149‑1 (MQ=255)
                                                       aGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAATAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGAcc                                                                                             <  7:139543/147‑1 (MQ=255)
                                                             tCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACATAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTcagca                                                                                                               <  6:266932/123‑1 (MQ=255)
                                                             tCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACATAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTcagca                                                                                                               >  5:266932/1‑123 (MQ=255)
                                                             tCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATCGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGtt                                                                                     >  1:105549/1‑149 (MQ=255)
                                                                      cacCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGa                                                                            <  4:289618/149‑1 (MQ=255)
                                                                           cAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTAct                                                                       >  8:9824/1‑149 (MQ=255)
                                                                                              cGTAGATCGACAGGCTCAAGAATTCATTGGAGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATTGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCTTATCCGAGTACTCTATAAAGCTGAGCTAt                                                      >  2:155669/1‑149 (MQ=37)
                                                                                                tAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCtttt                                                                                                                                       <  2:73266/66‑1 (MQ=255)
                                                                                                tAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCtttt                                                                                                                                       >  1:73266/1‑66 (MQ=255)
                                                                                                tAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGa                                                    <  5:202739/149‑1 (MQ=255)
                                                                                                 aGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAg                                                   >  2:338622/1‑149 (MQ=255)
                                                                                                    tCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTcc                                                <  5:85210/149‑1 (MQ=255)
                                                                                                      gACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCg                                               <  3:18501/148‑1 (MQ=255)
                                                                                                       aCAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGTCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGt                                              <  3:62772/148‑1 (MQ=255)
                                                                                                       aCAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGt                                                  <  5:19209/144‑1 (MQ=255)
                                                                                                       aCAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGt                                                  >  6:19209/1‑144 (MQ=255)
                                                                                                       aCAGGCTCAAGAATTCATAGGCGCACAAACAAAAAGCCCGCTCTATTCAGCGCGCATATTGAAGGCGAATGGCGCACTCAGCAGGATTCCAACCTGCGACCGAACGGTTCGTAGCCGAGAACTCTATCAAGCTGAGCTATGAGTCCGt                                              >  4:62772/1‑148 (MQ=255)
                                                                                                        cAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTTGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGa                                            <  5:158311/149‑1 (MQ=255)
                                                                                                        cAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGa                                            <  7:9824/149‑1 (MQ=255)
                                                                                                                aGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGtt                                    <  7:80143/149‑1 (MQ=255)
                                                                                                                 gAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGttt                                   <  6:158374/149‑1 (MQ=255)
                                                                                                                    ttCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGttt                                   >  8:7919/1‑146 (MQ=255)
                                                                                                                           gcgcACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGt                         >  2:200159/1‑149 (MQ=255)
                                                                                                                                  aacaaaAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACt                  <  1:200159/149‑1 (MQ=255)
                                                                                                                                     aaaaaTCCCGCTCTATTTAGCGGGCTTTTTGTTGGGGTATGGCGGACTGAGCATGATTCGAACCTGCGACCGCTCTGTTCGTAGCCGAGTACTCTATCCAGCTGATCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGc               <  1:155669/149‑1 (MQ=21)
                                                                                                                                      aaaaGCCCGCTCAATAGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAgg                          <  5:317197/137‑1 (MQ=255)
                                                                                                                                      aaaaGCCCGCTCAATAGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAgg                          >  6:317197/1‑137 (MQ=255)
                                                                                                                                         aGCCCGCTCAATAGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTAc           >  4:112268/1‑149 (MQ=255)
                                                                                                                                             cGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTAc           >  2:160308/1‑145 (MQ=255)
                                                                                                                                             cGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTAc           <  1:160308/145‑1 (MQ=255)
                                                                                                                                              gCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGtt                        <  6:315600/131‑1 (MQ=255)
                                                                                                                                              gCTCAATAGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTa            >  3:69037/1‑143 (MQ=255)
                                                                                                                                              gCTCAATAGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTa            <  4:69037/143‑1 (MQ=255)
                                                                                                                                                  tATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGc  >  3:311080/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ATCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGC  >  NC_002947/5078982‑5079274

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