Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,079,126 T→A 6.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,1260TA 6.7% 282.5 / 5.7 97intergenic (‑371/+26)mgtE/PP_t56magnesium transporter/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base T (39/51);  new base A (3/4);  total (42/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGC  >  NC_002947/5078981‑5079274
                                                                                                                                                  |                                                                                                                                                         
tATCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTAtt                                                                                                                                                        >  5:259550/1‑149 (MQ=255)
 aTCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTAtt                                                                                                                                                        >  3:311975/1‑148 (MQ=255)
  tCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTAtt                                                                                                                                                        <  4:43468/147‑1 (MQ=255)
  tCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGa                                                                                                                                                      <  3:2027/149‑1 (MQ=255)
     gTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCg                                                                                                                                                   >  1:320886/1‑149 (MQ=255)
       gCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCgg                                                                                                                                                  <  2:50275/148‑1 (MQ=255)
        cgacgaTTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCggg                                                                                                                                                 <  4:311975/148‑1 (MQ=255)
          acgaTTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCtt                                                                                                                                             <  7:312994/149‑1 (MQ=255)
          acgaTTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCACAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCtt                                                                                                                                             <  6:211746/149‑1 (MQ=255)
           cgaTTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCttt                                                                                                                                            <  7:73631/149‑1 (MQ=255)
            gaTTGTTGTGGATGTGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCtttt                                                                                                                                           >  5:255823/1‑149 (MQ=255)
            gaTTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCtttt                                                                                                                                           <  1:196538/149‑1 (MQ=255)
             aTTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCttttt                                                                                                                                          >  8:219367/1‑149 (MQ=255)
             aTTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCGGGCttttt                                                                                                                                          <  6:366589/149‑1 (MQ=255)
                 ttgtGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATg                                                                                                                                      >  2:209441/1‑149 (MQ=255)
                 ttgtGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATg                                                                                                                                      <  1:147775/149‑1 (MQ=255)
                 ttgtGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATg                                                                                                                                      >  5:114088/1‑149 (MQ=255)
                    tgGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATgggg                                                                                                                                   <  5:237934/149‑1 (MQ=255)
                    tgGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATggg                                                                                                                                    <  6:114088/148‑1 (MQ=255)
                          ggAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGTCAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGc                                                                                                                             <  1:71368/149‑1 (MQ=255)
                          ggAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATgg                                                                                                                              <  6:46935/148‑1 (MQ=255)
                          ggAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATgg                                                                                                                              <  3:190632/148‑1 (MQ=255)
                          ggAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGc                                                                                                                             <  2:209517/149‑1 (MQ=255)
                           gAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCg                                                                                                                            <  5:15579/149‑1 (MQ=255)
                            aaaGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCg                                                                                                                            >  7:23732/1‑148 (MQ=255)
                            aaaGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCg                                                                                                                            <  2:286134/148‑1 (MQ=255)
                            aaaGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCg                                                                                                                            >  1:286134/1‑148 (MQ=255)
                                cacaTTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACt                                                                                                                        >  8:281216/1‑148 (MQ=255)
                                 acaTTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACt                                                                                                                        >  3:344591/1‑147 (MQ=255)
                                   aTTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTcag                                                                                                                     >  7:178947/1‑148 (MQ=255)
                                    tttGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTcagca                                                                                                                   <  5:119280/149‑1 (MQ=255)
                                    tttGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCGGGCttttt                                                                                                                                          <  8:64843/126‑1 (MQ=255)
                                    tttGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCGGGCttttt                                                                                                                                          >  7:64843/1‑126 (MQ=255)
                                    tttGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCACAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCATTTTTGATGGGGTATGGCGGACTcagc                                                                                                                    <  1:12072/149‑1 (MQ=255)
                                     ttGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTcagcag                                                                                                                  >  5:76965/1‑149 (MQ=255)
                                        aTTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGAt                                                                                                               <  4:39367/149‑1 (MQ=255)
                                         tttttGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGAtt                                                                                                             >  6:19469/1‑149 (MQ=255)
                                             tGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGa                                                                                                          <  3:335857/148‑1 (MQ=255)
                                             cgAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGa                                                                                                          <  8:181621/147‑1 (MQ=255)
                                               aGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGTTCGaatcc                                                                                                      >  6:325853/1‑146 (MQ=255)
                                                ggCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑t                                                                                                   <  6:259550/149‑1 (MQ=255)
                                                ggCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑t                                                                                                   <  2:320886/149‑1 (MQ=255)
                                                 gCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑Tg                                                                                                  >  6:333879/1‑149 (MQ=255)
                                                 gCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑Tg                                                                                                  >  2:364918/1‑149 (MQ=255)
                                                     cTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGTTCGAATCC‑‑TGctga                                                                                              <  5:325853/149‑4 (MQ=255)
                                                           ggCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGATCAAGAATTCATTGCCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTATTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCg                                                                                        >  8:229789/1‑149 (MQ=255)
                                                              tCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGAATCAGCAGTATTCGAACC‑‑TGCGCCCGCTCGGtt                                                                                     >  6:275773/1‑149 (MQ=255)
                                                                cAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGtt                                                                                     <  5:752/147‑1 (MQ=255)
                                                                 aaaCTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGt                                                                                  <  1:209441/149‑1 (MQ=255)
                                                                  aaCTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGAt                                                                                                               <  6:82806/123‑1 (MQ=255)
                                                                  aaCTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGAt                                                                                                               >  5:82806/1‑123 (MQ=255)
                                                                  aaCTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTa                                                                                 <  7:67976/149‑1 (MQ=255)
                                                                   aCTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAg                                                                                <  2:348680/149‑1 (MQ=255)
                                                                    cTACCCCCCAAGCATTCCAGCAGAGCGTAGATCCACAGGCTCAACAATTCATTGGCGCCCAAACAAAAACCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGc                                                                               <  4:322055/149‑1 (MQ=37)
                                                                       cacCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGa                                                                            <  1:285294/149‑1 (MQ=255)
                                                                          cGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGt                                                                          >  4:288723/1‑148 (MQ=255)
                                                                           gCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTAc                                                                        >  5:129919/1‑149 (MQ=255)
                                                                               tCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGc                                                                                                 >  8:204519/1‑120 (MQ=255)
                                                                               tCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGc                                                                                                 <  7:204519/120‑1 (MQ=255)
                                                                               tCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTa                                                                    <  4:90603/149‑1 (MQ=255)
                                                                                      agaagaGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCt                                                             <  7:281216/149‑1 (MQ=255)
                                                                                      agaagaGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCt                                                             >  1:16079/1‑149 (MQ=255)
                                                                                          gagCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGc                                                         <  8:178947/149‑1 (MQ=255)
                                                                                            gctgaGATCGACAGGCTCATGAATTCATTGGCTCACATACAAAAAGCGCGCTCTTTTGAGCGGGCTTTTTGTTGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCGGAGTACTCTATCCAGCTGAGCTa                                                       <  5:380896/145‑1 (MQ=25)
                                                                                                 gATCGACAGGCTCAAGAATTGATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCTTAGCCGAGTACTCTATCCAGCTGAGCTATGAGt                                                  <  8:23732/149‑1 (MQ=255)
                                                                                                    cGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCg                                               >  3:93928/1‑149 (MQ=255)
                                                                                                     gACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGt                                              >  1:22728/1‑149 (MQ=255)
                                                                                                       cAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCt                                                        >  8:106487/1‑137 (MQ=255)
                                                                                                       cAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCt                                                        <  7:106487/137‑1 (MQ=255)
                                                                                                        aGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGAt                                           <  4:93928/149‑1 (MQ=255)
                                                                                                          gCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATgg                                         <  2:16079/149‑1 (MQ=255)
                                                                                                            tCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCATCAATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATgg                                         >  1:57719/1‑148 (MQ=255)
                                                                                                            tCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTa                                                                         <  6:1853/115‑1 (MQ=255)
                                                                                                            tCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTa                                                                         >  5:1853/1‑115 (MQ=255)
                                                                                                            tCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATg                                                     <  3:356079/135‑1 (MQ=255)
                                                                                                            tCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATg                                                     >  4:356079/1‑135 (MQ=255)
                                                                                                            tCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTa                                       >  6:339666/1‑149 (MQ=255)
                                                                                                              aaGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGt                                     >  5:249648/1‑149 (MQ=255)
                                                                                                              aaGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGAACGCTAGGTTAGTAGGCGTGTAAGCTATCCATCTGAGCTATGTGTCAGTGATGGTAGt                                     >  1:236301/1‑149 (MQ=37)
                                                                                                              aaGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGAGAACGCACGCGTCGTCGCCGACTAAGCACTCCAGCCGCGATATCAGGACTTGATGGTAGt                                     >  1:230688/1‑149 (MQ=255)
                                                                                                                gAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGttt                                   <  7:219367/149‑1 (MQ=255)
                                                                                                                      aTTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGAGATTGTAGttttttaa                              >  5:384545/1‑144 (MQ=255)
                                                                                                                           cgcACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGtt                        <  6:76965/149‑1 (MQ=255)
                                                                                                                           cgcACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGtt                        <  3:288723/149‑1 (MQ=255)
                                                                                                                              acaaacaaAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGAGATTGGAAGTTTAAACCaagataac                     >  8:68979/1‑142 (MQ=18)
                                                                                                                                aaacaaaAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCAc                   >  4:283038/1‑149 (MQ=255)
                                                                                                                                   caaaAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGTGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTgg                <  6:255823/149‑1 (MQ=255)
                                                                                                                                    aataaGCCCGCTCTATTGATCGTGCTTTTTGATGGGGTATGGCTGACTCAGCAGGATTGGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCATCTGTGCTATGAGTCCGTGATGGTATTTTTAGACCAGGTTACCACTGGc               <  3:335466/146‑1 (MQ=21)
                                                                                                                                    aaaaaGCCCGGTCTATTGAGCGGGGTTTTTTATGGGGTTTGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGc               <  2:22728/149‑1 (MQ=255)
                                                                                                                                    aaaaaGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGc               <  1:364918/149‑1 (MQ=255)
                                                                                                                                     aaaaGCCCGCTCAATAGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTAGTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCt              >  3:37289/1‑149 (MQ=255)
                                                                                                                                     aaaaGCCCGCTCAATAGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTa                       <  1:157216/140‑1 (MQ=255)
                                                                                                                                       aaGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATCTAGAACCATTTATTTCCCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGc               <  1:126987/149‑1 (MQ=21)
                                                                                                                                       aaGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCtt             >  7:318578/1‑148 (MQ=255)
                                                                                                                                           ccGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTg         >  7:42444/1‑148 (MQ=255)
                                                                                                                                             gCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACt          <  2:289153/145‑1 (MQ=255)
                                                                                                                                             gCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATGCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACt          >  1:289153/1‑145 (MQ=255)
                                                                                                                                                  tATTGAGCGGGCTTTTTGATGGGGTATGGAGGACTCAGAAGGATTCGAACA‑‑TGCGACCGCTCGGTTCTTAGACGAGTACTGGATAAAGGTGAGATATGAGTAAGTGATGTTAGTTTTAGAGGAGGTTACCGGTGGCTTAGTGAAATGGc  >  8:247974/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                         
TATCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACC‑‑TGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGC  >  NC_002947/5078981‑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: