Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,655,155 T→A 5.6% S131C (AGC→TGC)  epd ← D‑erythrose 4‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,655,1550TA 5.6% 216.9 / 2.2 71S131C (AGC→TGC) epdD‑erythrose 4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (37/30);  new base A (2/2);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

ATTCGATGCCAAACGCCTGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGCCCAGCGCAGCCCAGTCGATGCCTTCGGGGGTGGCACTGC  >  NC_002947/5655008‑5655301
                                                                                                                                                   |                                                                                                                                                  
atTCGATGCCAAACGCCTGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTg                                                                                                                                                   >  4:83312/1‑149 (MQ=255)
  tCGATGCCAAACGCCTGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGc                                                                                                                                                 >  6:246075/1‑149 (MQ=255)
    gATGCCAAACGCCTGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCa                                                                                                                                               >  5:24012/1‑149 (MQ=255)
         cAAACGCCTGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGc                                                                                                                                          >  8:239976/1‑149 (MQ=255)
         cAAACGCCTGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGc                                                                                                                                          >  6:161329/1‑149 (MQ=255)
         cAAACGCCTGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCAGCCCATggc                                                                                                                                           >  8:56300/1‑147 (MQ=255)
          aaaCGCCTGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATgg                                                                                                                                            <  4:225119/146‑1 (MQ=255)
          aaaCGCCTGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATgg                                                                                                                                            >  3:225119/1‑146 (MQ=255)
          aaaCGCCTGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCt                                                                                                                                         >  2:112745/1‑149 (MQ=255)
          aaaCGCCTGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCt                                                                                                                                         >  2:203149/1‑149 (MQ=255)
             cGCCTGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGc                                                                                                                                       >  1:158096/1‑148 (MQ=255)
             cGCCTGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCg                                                                                                                                      <  7:29897/149‑1 (MQ=255)
               ccTGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCgag                                                                                                                                    <  8:178687/149‑1 (MQ=255)
                cTGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCAGCCCATGGCCAGCgagg                                                                                                                                   <  7:56300/149‑2 (MQ=255)
                 tGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGaa                                                                                                                                  >  8:292536/1‑149 (MQ=255)
                      cAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCa                                                                                                                             >  1:325386/1‑149 (MQ=255)
                      cAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCa                                                                                                                             <  3:83312/149‑1 (MQ=255)
                       aGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCAc                                                                                                                            >  5:160206/1‑149 (MQ=255)
                        gCACACGCAGCAGCGGCACGCCGCACTTGGTGGTGCACGAGGCGATGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCAATGGGCTGCGAGAACAGCAcc                                                                                                                           >  3:169594/1‑149 (MQ=255)
                         cacaCGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGACGGCCTCGCTGGCCATGGGCTGCGAGAACAGCAccc                                                                                                                          <  8:166694/149‑1 (MQ=255)
                           cacGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCtt                                                                                                                        <  5:302917/149‑1 (MQ=255)
                              gcagcaGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGgc                                                                                                                     <  7:239976/149‑1 (MQ=255)
                              gcagcaGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGgc                                                                                                                     <  2:207609/149‑1 (MQ=255)
                                   agcGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCAcc                                                                                                                <  5:246075/149‑1 (MQ=255)
                                   agcGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCGTGccgatgg                                                                                                                >  3:4159/1‑145 (MQ=255)
                                       gCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCt                                                                                                            >  6:255561/1‑149 (MQ=255)
                                           gccgcAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGa                                                                                                        <  1:211099/149‑1 (MQ=255)
                                           gccgcAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGa                                                                                                        <  7:292536/149‑1 (MQ=255)
                                            ccgcAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCa                                                                                                                             <  1:86892/127‑1 (MQ=255)
                                            ccgcAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTAGGCCTCGCTGGCCATGGGCTGCGAGAACAGCa                                                                                                                             >  2:86892/1‑127 (MQ=255)
                                                  ctggtggTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGATAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACg                                                                                                 >  6:106499/2‑149 (MQ=255)
                                                    ggtggtGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTACGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCt                                                                                               >  4:246335/1‑149 (MQ=255)
                                                        gtgCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGc                                                                                            >  6:69426/1‑148 (MQ=255)
                                                           cACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATc                                                                                        <  5:106499/149‑1 (MQ=255)
                                                                  cGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACgg                                                                                 <  1:203149/149‑1 (MQ=255)
                                                                      ggACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGgtg                                                                              <  6:306648/148‑1 (MQ=255)
                                                                       gACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGgtgt                                                                             <  6:24012/148‑1 (MQ=255)
                                                                         cacCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCga                                                                                                                                     <  8:274747/90‑1 (MQ=255)
                                                                         cacCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCga                                                                                                                                     >  7:274747/1‑90 (MQ=255)
                                                                          acCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTta                                                                           >  5:94042/1‑147 (MQ=255)
                                                                          acCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTGCGAGGAAACGCTGGCCATCGGCACGGGTGTta                                                                           >  3:34513/1‑147 (MQ=255)
                                                                            cAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAgg                                                                       >  4:5535/1‑149 (MQ=255)
                                                                                           ggTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCa                                                        <  2:189129/149‑1 (MQ=255)
                                                                                            gTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAg                                                       >  5:317160/1‑149 (MQ=255)
                                                                                            gTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAg                                                       <  5:98819/149‑1 (MQ=255)
                                                                                               aggcaggcCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCAc                                                    <  5:255561/149‑1 (MQ=255)
                                                                                                   aggcCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAgg                                                >  5:49005/1‑149 (MQ=255)
                                                                                                    ggcCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGt                                               <  2:325386/149‑1 (MQ=255)
                                                                                                      ccTGGTTGTTGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGGACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCg                                             <  4:169594/149‑1 (MQ=255)
                                                                                                       cTGGTTGATGCCGTAGACACCCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGa                                            >  2:28873/1‑149 (MQ=255)
                                                                                                          gTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGc                                                                                                >  5:337047/1‑94 (MQ=255)
                                                                                                          gTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGc                                                                                                <  6:337047/94‑1 (MQ=255)
                                                                                                          gTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGc                                         >  6:251881/1‑149 (MQ=255)
                                                                                                             gATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGcc                                        >  3:108157/1‑147 (MQ=255)
                                                                                                                gCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGCCCAgcg                                   >  4:51389/1‑149 (MQ=255)
                                                                                                                  cGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCTGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGa                                            <  8:241897/138‑1 (MQ=255)
                                                                                                                  cGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGa                                            >  7:241897/1‑138 (MQ=255)
                                                                                                                   gTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGCCCAGCGCAg                                <  2:158096/149‑1 (MQ=255)
                                                                                                                         aaCCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGCCCAGCGCAGCCCAGt                          <  1:125420/149‑1 (MQ=255)
                                                                                                                          acagTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGCCCAGCGCAGCCCAGTc                         >  2:152728/4‑149 (MQ=255)
                                                                                                                          aCCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGCCCAGCGCAGCCCAGTc                         <  6:94042/149‑1 (MQ=255)
                                                                                                                            cGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGCCCAGCGCAGCCCAGTCGa                       >  2:50699/1‑149 (MQ=255)
                                                                                                                            cGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCAACGGCTTCGAGGAAACGCTTGCCATAGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGCCCAGCGCAGCCCAGTCGa                       >  6:231398/1‑149 (MQ=255)
                                                                                                                                       cGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAg                                                                        <  4:196601/89‑1 (MQ=255)
                                                                                                                                       cGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAg                                                                        >  3:196601/1‑89 (MQ=255)
                                                                                                                                       cGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGCCCAGCGCAGCCCAGTCGACGCCTTCgggg            <  3:51389/149‑1 (MQ=255)
                                                                                                                                          cGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGCCCAGCGCAGCCCAGTCGATGCCTTCGGGGGTg         >  4:212147/1‑149 (MQ=255)
                                                                                                                                           gttcTCGCAGCCCATGGCCAGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGCCCAGCGCAGCCCAGTCGATGCCTTCggg             >  6:990/4‑144 (MQ=255)
                                                                                                                                           gttcTCGCAGCCCATGGCCAGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGCCCAGCGCAGCCCAGTCGATGCCTTCggg             <  5:990/141‑1 (MQ=255)
                                                                                                                                              cTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGCCCAGCGCAGCCCAGTCGa                       >  4:233500/1‑131 (MQ=255)
                                                                                                                                                 gCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGCCCAGCGCAGCCCAGTCGATGCCTTCGGGGGTGGCACTgc  <  5:251881/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATTCGATGCCAAACGCCTGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGCCCAGCGCAGCCCAGTCGATGCCTTCGGGGGTGGCACTGC  >  NC_002947/5655008‑5655301

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