Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,079,129 T→A 14.9% intergenic (‑374/+23) 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,1290TA14.9% 163.2 / 18.2 67intergenic (‑374/+23)mgtE/PP_t56magnesium transporter/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base T (22/35);  new base A (5/5);  total (27/40)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

ATCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGCTT  >  NC_002947/5078982‑5079276
                                                                                                                                                     |                                                                                                                                                   
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                                                                                                                                         >  1:73266/1‑66 (MQ=255)
                                                                                                tAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCtttt                                                                                                                                         <  2:73266/66‑1 (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                                                    >  6:19209/1‑144 (MQ=255)
                                                                                                       aCAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGt                                                    <  5:19209/144‑1 (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                            >  6:317197/1‑137 (MQ=255)
                                                                                                                                      aaaaGCCCGCTCAATAGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAgg                            <  5:317197/137‑1 (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              <  4:69037/143‑1 (MQ=255)
                                                                                                                                              gCTCAATAGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTa              >  3:69037/1‑143 (MQ=255)
                                                                                                                                                  tATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGc    >  3:311080/1‑149 (MQ=255)
                                                                                                                                                    ttGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGCtt  <  4:311080/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ATCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATT‑‑ACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGCTT  >  NC_002947/5078982‑5079276

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: