Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,103,269 G→T 6.3% A10S (GCT→TCT)  murA → UDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,103,2690GT 6.3% 201.2 / 3.1 63A10S (GCT→TCT) murAUDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base G (30/29);  new base T (2/2);  total (32/31)
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.00e-01

ATCTATGCTCACCTGAATCCGTGGATCGCCAATGGCAGCATCCATGCGGTAACCATGAAATTTTTCAGCAGCGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCACCGTGGGCAACCTGCCGCACCTGCACGACATCACCACCATGATCGAGCTGTTC  >  NC_002947/1103123‑1103412
                                                                                                                                                  |                                                                                                                                               
aTCTATGCTCACCTGAATCCGTGGATCGCCAATGGCAGCATCCATGCGGTAACCATGAAATTTTTCAGCAGCGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCt                                                                                                                                               >  5:100477/1‑149 (MQ=255)
  cTATGCTCACCTGAATCCGTGGATCGCCAATGGCAGCATCCATGCGGTAACCATGAAATTTTTCAGCAGCGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGATCTCg                                                                                                                                             >  6:243534/1‑149 (MQ=255)
     tGCTCACCTGAATCCGTGGATCGCCAATGGCAGCATCCATGCGGTAACCATGAAATTTTTCAGCAGCGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGTCAAACTGATTATTACTGGCGGCGCt                                                                                                                                               <  4:240920/144‑1 (MQ=255)
     tGCTCACCTGAATCCGTGGATCGCCAATGGCAGCATCCATGCGGTAACCATGAAATTTTTCAGCAGCGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTc                                                                                                                                              <  1:174002/145‑1 (MQ=255)
          aCCTGAATCCGTGGATCGCCAATGGCAGCATCCATGCGGTAACCATGAAATTTTTCAGCAGCGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGa                                                                                                                                       <  2:178313/147‑1 (MQ=255)
                 tCCGTGGATCGCCAATGGCATCATCCATTCTGTAACCATGAAATTTTTCCGCAGCGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGAt                                                                                                                              <  4:194798/149‑1 (MQ=255)
                 tCCGTGGATCGCCAATGGCAGCATCCATGCGGTAACCATGAAATTTTTCAGCAGCGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGc                                                                                                                                   >  1:33469/1‑144 (MQ=255)
                 tCCGTGGATCGCCAATGGCAGCATCCATGCGGTAACCATGAAATTTTTCAGCAGAGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGATCTCGCCGTCAAGGCgagcg                                                                                                                              >  3:194760/1‑147 (MQ=255)
                    gTGGATCGCCAATGGCAGCATCCATGCGGTAACCATGAAATTTTTCAGCAGCGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCg                                                                                                                           >  8:136553/1‑149 (MQ=255)
                     tGGATCGCCAATGGCAGCATCCATGCGGTAACCATGAAATTTTTCAGCAGCGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGATCTCGCCGTCAAGGCGAGcgccgc                                                                                                                          <  4:194760/149‑1 (MQ=255)
                       gATCGCCAATGGCAGCATCCATGCGGTAACCATGAAATTTTTCAGCAGCGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCa                                                                                                                         <  8:221710/148‑1 (MQ=255)
                       gATCGCCAATGGCAGCATCCATGCGGTAACCATGAAATTTTTCAGCAGCGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGATCTCGCCGTCAAGGCGAGcgccgcca                                                                                                                        <  5:243534/149‑3 (MQ=255)
                           gCCAATGGCAGCATCCATGCGGTAACCATGAAATTTTTCAGCAGCGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCg                                                                                                                    >  4:268011/1‑149 (MQ=255)
                            ccAATGGCAGCATCCATGCGGTAACCATGAAATTTTTCAGCAGCGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGATAAACTGATTATTACAGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCgg                                                                                                                   >  1:284367/1‑149 (MQ=255)
                            ccAATGGCAGCATCCATGCGGTAACCATGAAATTTTTCAGCAGCGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCgg                                                                                                                   >  6:249987/1‑149 (MQ=255)
                            ccAATGGCAGCATCCATGCGGTAACCATGAAATTTTTCAGCAGCGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCgg                                                                                                                   <  3:268011/149‑1 (MQ=255)
                                   cagcaTCCATGCGGTAACCATGAAATTTTTCAGCAGCGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGa                                                                                                             >  8:74733/1‑148 (MQ=255)
                                    agcaTCCATGCGGTAACCATGAAATTTTTCAGCAGCGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGgctgcgtt                                                                                                           >  5:175033/1‑143 (MQ=255)
                                    agcaTCCATGCGGTAACCATGAAATTTTTCAGCAGCGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGaa                                                                                                            >  2:262795/1‑148 (MQ=255)
                                                    ccATGAAATTTTTCAGCAGCGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGa                                                                                           >  8:10812/1‑149 (MQ=255)
                                                     cATGAAATTTTTCAGCAGCGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGccc                                                                                                 >  3:207878/1‑142 (MQ=255)
                                                     cATGAAATTTTTCAGCAGCGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGccc                                                                                                 <  4:207878/142‑1 (MQ=255)
                                                     cATGAAATTTTTCAGCAGCGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGccag                                                                                           >  3:30039/1‑146 (MQ=255)
                                                            tttttCAGCAGCGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTAGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCa                                                                                   <  1:262795/149‑1 (MQ=255)
                                                               ttCAGCAGCGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGc                                                                                 <  2:183324/148‑1 (MQ=255)
                                                                      gcgcagcCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCctgc                                                                         >  4:286247/1‑149 (MQ=255)
                                                                        gcagcCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCctgctg                                                                       <  6:100477/149‑1 (MQ=255)
                                                                          agcCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGc                                                                     <  4:12000/149‑1 (MQ=255)
                                                                          agcCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACAGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGc                                                                     <  2:284367/149‑1 (MQ=255)
                                                                              tGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGcc                                                                                                                                           <  3:226765/75‑1 (MQ=255)
                                                                              tGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGcc                                                                                                                                           >  4:226765/1‑75 (MQ=255)
                                                                                ggCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATg                                                                <  7:282645/148‑1 (MQ=255)
                                                                                  cTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGc                                                                                                                                            >  6:283964/1‑70 (MQ=255)
                                                                                  cTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGc                                                                                                                                            <  5:283964/70‑1 (MQ=255)
                                                                                  cTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGcc                                                             <  5:249987/149‑1 (MQ=255)
                                                                                     aGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCg                                                           <  5:94619/148‑1 (MQ=255)
                                                                                          aCCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCAcc                                                     >  1:180935/1‑149 (MQ=255)
                                                                                          aCCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCAcc                                                     >  8:15085/1‑149 (MQ=255)
                                                                                          aCCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCAcc                                                     >  2:224433/1‑149 (MQ=255)
                                                                                            cTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCACCGt                                                   <  7:15085/149‑1 (MQ=255)
                                                                                               aGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCACCGTgg                                                 <  2:180935/148‑1 (MQ=255)
                                                                                                   aaCTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCACCGTGGGCAAc                                            >  1:93360/1‑149 (MQ=255)
                                                                                                   aaCTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCACGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCACCGTGGGCAAc                                            >  2:101452/1‑149 (MQ=255)
                                                                                                    aCTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCACCGTGGGCAAcc                                           <  7:74733/149‑1 (MQ=255)
                                                                                                        ggcggcGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCACCGTGGGCAACCTGc                                        <  2:33469/148‑1 (MQ=255)
                                                                                                            gcgAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCAccc                                                                            <  2:34904/108‑1 (MQ=255)
                                                                                                            gcgAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCAccc                                                                            >  1:34904/1‑108 (MQ=255)
                                                                                                             cgAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTATGGCAGCCACCCTGCTGGCCGATGGCCCGGTCACCGTGGGCAACCTGCCGCAcc                                  >  4:224344/1‑149 (MQ=255)
                                                                                                                 aCTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGACATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCACCGTGGGCAACCTGCCGCACCTGCa                              >  3:60327/1‑149 (MQ=255)
                                                                                                                  cTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCACGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCACCGTGGGCAACCTGCCGCACCTGCAc                             <  1:101452/149‑1 (MQ=255)
                                                                                                                     cAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCACCGTggg                                                <  8:80492/127‑1 (MQ=255)
                                                                                                                     cAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCACCGTggg                                                >  7:80492/1‑127 (MQ=255)
                                                                                                                     cAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCACCGTGGGCAACCTGCCGCACCTGCACg                            <  1:302029/147‑1 (MQ=255)
                                                                                                                     cAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCACCGTGGGCAACCTGCCGCACCTGCACGa                           >  3:137261/1‑148 (MQ=255)
                                                                                                                     cAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCACCGTGGGCAACCTGCCGCACCTGCACGAc                          <  7:10812/149‑1 (MQ=255)
                                                                                                                      aaTGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCACCGTGGGCAACCTGCCGCACCTGCACGACa                         >  1:60635/1‑149 (MQ=255)
                                                                                                                            cAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCACCGTGGGCAACCTGCCGCACCTGCACGACATcacca                   <  5:111732/149‑1 (MQ=255)
                                                                                                                               aCTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCACCGTGGGCAACCTGc                                        >  1:121325/1‑125 (MQ=255)
                                                                                                                               aCTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCACCGTGGGCAACCTGc                                        <  2:121325/125‑1 (MQ=255)
                                                                                                                                   attattACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCACCGTGGGCAACCTGCCGCACCTGCACGACATCACCACCATGa             >  2:165959/1‑148 (MQ=255)
                                                                                                                                       ttaCTGGCGGCGCTCGCCTTGACGGCGACATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCACCGTGGGCAACCTGCCGCACCTGCACGACATCACCACCATGATCGa         <  4:60327/148‑1 (MQ=255)
                                                                                                                                          cTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCACCGTGGGCAACCTGCCGCACCTGCACGACATCACCACCATGATCGAGCTg     >  4:197632/1‑149 (MQ=255)
                                                                                                                                           tGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCACCGTGGGCAACCTGCCGCACCTGCACGACATCACCACCATGATCGAGCTGt    >  5:97618/1‑149 (MQ=255)
                                                                                                                                             gcggcgCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCACCGTGGGCAACCTGCCGCACCTGCACGACATCACCACCATGATCGAGCTGTTc  >  7:186163/1‑149 (MQ=255)
                                                                                                                                             gcggcgCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGACCCGGTCACCGTGGGCAACCTGCCGCACCTGCACGCCATCACCACCATGATCGAGCTGTTc  >  7:266559/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
ATCTATGCTCACCTGAATCCGTGGATCGCCAATGGCAGCATCCATGCGGTAACCATGAAATTTTTCAGCAGCGCAGCCTGGGCTGAGCGCACCTGAGCCAACTTGGCGGCGAGACTCCAATGGACAAACTGATTATTACTGGCGGCGCTCGCCTTGACGGCGAGATCCGCATTTCGGGCGCGAAAAACGCAGCCCTGCCGATTCTGGCAGCCACCCTGCTGGCCGATGGCCCGGTCACCGTGGGCAACCTGCCGCACCTGCACGACATCACCACCATGATCGAGCTGTTC  >  NC_002947/1103123‑1103412

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