Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,161,003 A→G 6.2% V247V (GTA→GTG)  pabC → 4‑amino‑4‑deoxychorismate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,161,0030AG 6.2% 191.9 / 2.3 64V247V (GTA→GTG) pabC4‑amino‑4‑deoxychorismate lyase
Reads supporting (aligned to +/- strand):  ref base A (30/30);  new base G (2/2);  total (32/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.92e-01

GCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCCTGCTGTTGCTGGAAATGGGCTTGATCCTTGCCGGCCTGGCGTTTGGCT  >  NC_002947/2160865‑2161146
                                                                                                                                            |                                                                                                                                                
gCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTa                                                                                                                                         >  6:156457/1‑148 (MQ=255)
 cGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTaaa                                                                                                                                       <  8:13988/149‑1 (MQ=255)
     ccGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTg                                                                                                                                    <  2:28643/148‑1 (MQ=255)
       gCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTc                                                                                                                                 <  3:190644/149‑1 (MQ=255)
        cACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGATGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCg                                                                                                                                >  8:102885/1‑149 (MQ=255)
             ctgGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCggg                                                                                                                           >  5:292422/1‑149 (MQ=255)
             ctgGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCgg                                                                                                                            <  7:86125/148‑1 (MQ=255)
              gggcaccggCCCCGCTGCTGGGTATTGCGGTACAGGGCTGCGACCTGTCGTTCGAGGCGCTGGATCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGc                                                                                                                          <  6:17681/142‑1 (MQ=255)
                gAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGACGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGc                                                                                                                          <  6:12487/147‑1 (MQ=255)
                  aCAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCAGGGCCAGc                                                                                                                      >  1:199752/1‑149 (MQ=255)
                     ggCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCa                                                                                                                   >  5:39016/1‑149 (MQ=255)
                      gCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCa                                                                                                                   >  6:258701/1‑148 (MQ=255)
                      gCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCa                                                                                                                   <  5:258701/148‑1 (MQ=255)
                               ctgGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACt                                                                                                         >  3:15278/1‑149 (MQ=255)
                                  ggTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTTGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGc                                                                                                       <  5:156457/148‑1 (MQ=255)
                                  ggTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTg                                                                                                        <  7:102885/147‑1 (MQ=255)
                                  ggTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCa                                                                                                      <  8:7698/149‑1 (MQ=255)
                                      ttGCGGTACATGTCTTCGACCTGTCGTTCGAGGCGCTGGAGCAAGCTAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGATAAACTGGGCGACGGGCCCGCTCACTCGTAAACTGCAGGc                                                                                                   <  2:199752/148‑1 (MQ=255)
                                         cGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑tt                                                                                              <  3:144247/149‑1 (MQ=255)
                                                 gTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTAc                                                                                      >  7:304282/1‑149 (MQ=255)
                                                   cTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGt                                                                                    >  4:13894/1‑149 (MQ=255)
                                                    tGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGt                                                                                    >  2:216345/1‑148 (MQ=255)
                                                            gTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGaaa                                                                           >  3:197738/1‑149 (MQ=255)
                                                               gTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCCGGCCG‑TTGCCCGTACGTTACTGGAAACCt                                                                        >  4:20333/1‑149 (MQ=255)
                                                               gTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCt                                                                        >  6:94930/1‑149 (MQ=255)
                                                               gTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCt                                                                        >  4:288891/1‑149 (MQ=255)
                                                               gTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCt                                                                        >  5:333031/1‑149 (MQ=255)
                                                                  cGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGa                                                                      <  3:20333/148‑1 (MQ=255)
                                                                     ggCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTc                                                                  >  3:38074/1‑149 (MQ=255)
                                                                       cgcTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGt                                                                <  1:216345/149‑1 (MQ=255)
                                                                       cgcTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGt                                                                <  6:291609/149‑1 (MQ=255)
                                                                         cTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑GGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGt                                                                >  2:286107/1‑147 (MQ=255)
                                                                         cTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTga                                                              <  5:62182/149‑1 (MQ=255)
                                                                              gCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACg                                                          >  7:313325/1‑148 (MQ=255)
                                                                                 aGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGTGAAGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACATGAATTCGTGAGACGt                                                         >  1:101100/1‑148 (MQ=255)
                                                                                    gAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTc                                                   <  3:13894/149‑1 (MQ=255)
                                                                                       tGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCctg                                                <  8:313325/149‑1 (MQ=255)
                                                                                         aattATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCctgc                                               <  5:43507/145‑1 (MQ=255)
                                                                                           gTATTTGTTTGCAATAGCGTTTACGGCATTAGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCCTGCTGt                                            >  1:29590/1‑149 (MQ=255)
                                                                                              tttgtttgCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCCTGCTGTTGc                                         >  8:107962/1‑149 (MQ=255)
                                                                                                tgtttgCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCCTGCTGTTGCTg                                       >  6:153640/1‑149 (MQ=255)
                                                                                                 gtttgCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑GGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCGCTTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCCTGCTGTTGCt                                        <  1:286107/148‑1 (MQ=255)
                                                                                                 gtttgCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCCTGCTGTTGCTg                                       <  6:39016/148‑1 (MQ=255)
                                                                                                     gCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCCTGCTGTTGCTGGAAAt                                  >  8:4525/1‑149 (MQ=255)
                                                                                                       aaTAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCCTGCTGTTGCTGGAAATgg                                <  4:15278/149‑1 (MQ=255)
                                                                                                         tAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCCTGCTGTTGCTGGAAATGGGc                              >  6:64413/1‑149 (MQ=255)
                                                                                                           gCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCCTGCTGTTGCTGGAAATGGGCt                             <  5:94930/148‑1 (MQ=255)
                                                                                                           gCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCCTGCTGTTGCTGGAAATGGGCt                             <  8:304282/148‑1 (MQ=255)
                                                                                                              tttACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCCTGCTGTTGCTGGAAATGGGCTTGAt                         >  5:59713/1‑149 (MQ=255)
                                                                                                               ttACGGCATTTGGCCCGTGCGTGGCGT‑‑GGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCct                                                 >  6:79650/1‑124 (MQ=255)
                                                                                                               ttACGGCATTTGGCCCGTGCGTGGCGT‑‑GGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCct                                                 <  5:79650/124‑1 (MQ=255)
                                                                                                                  cGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCCTGCTGTTGCTGGAAATGGGCTTGATCCtt                     <  6:292422/149‑1 (MQ=255)
                                                                                                                     cATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCCTGCTGTTGCTGGAAATGGGCTTGATCCTTGcc                  >  4:303004/1‑149 (MQ=255)
                                                                                                                       tttGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCCTGCTGTTGc                                         <  2:289651/124‑1 (MQ=255)
                                                                                                                       tttGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTACAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCCTGCTGTTGc                                         >  1:289651/1‑124 (MQ=255)
                                                                                                                          ggCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGGCGTAAATTCCTGCTGTTGCTGGAAATGGGCTTGATCCTTGCCGGcc              <  4:197738/148‑1 (MQ=255)
                                                                                                                          ggCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCCTGCTGTTGCTGGAAATGGGCTTGATCCTTGCCGGcc              <  3:288891/148‑1 (MQ=255)
                                                                                                                             ccGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCCTGCTGTTGCTgg                                      <  1:26166/121‑1 (MQ=255)
                                                                                                                             ccGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCCTGCTGTTGCTgg                                      >  2:26166/1‑121 (MQ=255)
                                                                                                                             ccGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCATGCTGTTGCTGGAAATGGGCTTGATCCTTGCCGGCCTGGc          >  6:137143/1‑149 (MQ=255)
                                                                                                                                    tgGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCCTGCTGTTGCTGGAAATGGGCTTGATCCTTGCCGGCCTGGCGTTTGGc   <  6:333031/149‑1 (MQ=255)
                                                                                                                                     ggCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCCTGCTGTTGCTGGAAATGGGCTTGATCCTTGCCGGCCTGGCGTTTGGCt  >  5:319510/1‑149 (MQ=255)
                                                                                                                                         t‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCCTGCTGTTGCTGGAAATggg                               <  7:237798/116‑1 (MQ=255)
                                                                                                                                         t‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCCTGCTGTTGCTGGAAATggg                               >  8:237798/1‑116 (MQ=255)
                                                                                                                                            |                                                                                                                                                
GCGTGCCGCACTGCTGGAACAGGCCCCGCTGCTGGGTATTGCGGTACAGGTCTGCGACCTGTCGTTCGAGGCGCTGGAGCAAGCGAATGAAGTATTTGTTTGCAATAGCGTTTACGGCATTTGGCCCGTGCGTGGCGT‑‑AGCCGCGCTAAACTGGTCGCCGGGCCCGCTCACTCGTAAACTGCAGGCCG‑TTGCCCGTACGTTACTGGAAACCTGAATTCGTGAGACGTAAATTCCTGCTGTTGCTGGAAATGGGCTTGATCCTTGCCGGCCTGGCGTTTGGCT  >  NC_002947/2160865‑2161146

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: ‑