Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,063,465 G→C 5.2% L118L (CTG→CTC ubiC → chorismate pyruvate‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,063,4650GC 5.2% 246.1 / 3.1 77L118L (CTG→CTCubiCchorismate pyruvate‑lyase
Reads supporting (aligned to +/- strand):  ref base G (34/39);  new base C (2/2);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-01

GCATCGCTGCCGGCGCCGAAGGCTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTGGGGCCGCCGCTCGCGCTTCGAGCGCAACGGCCTGGACCTGCTG  >  NC_002947/6063320‑6063612
                                                                                                                                                 |                                                                                                                                                   
gCATCGCTGCCGGCGCCGAAGGCTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGc                                                                                                                                                  >  2:273933/1‑149 (MQ=255)
 cATCGCTGCCGGCGCCGAAGGCTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGc                                                                                                                                                  >  4:193095/1‑148 (MQ=255)
     gCTGCCGGCGCCGAAGGCTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAgc                                                                                                                                              <  1:278106/148‑1 (MQ=255)
       tGCCGGCGCCGAAGGCTGGGTGCGCGAAGTGAATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAgctgc                                                                                                                                           >  1:61850/1‑149 (MQ=255)
          cGGCGCCGAAGGCTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGt                                                                                                                                        >  3:31234/1‑149 (MQ=255)
             cgcCGAAGGCTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCt                                                                                                                                     <  1:24022/149‑1 (MQ=255)
              gcCGAAGGCTGGGTGCGCGGAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTg                                                                                                                                    >  7:326327/1‑149 (MQ=255)
              gcCGAAGGCTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTg                                                                                                                                    <  3:193095/149‑1 (MQ=255)
              gcCGAAGGCTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTg                                                                                                                                    >  4:252787/1‑149 (MQ=255)
               ccGAAGGCTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCAGCTCGCCCAGCGAGcggctgcccaag                                                                                                                                   <  5:334735/149‑7 (MQ=255)
                 gAAGGCTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGa                                                                                                                                 >  7:243952/1‑149 (MQ=255)
                 gAAGGCTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCAGCTCGCCCAGCGAGcggctgcccaa                                                                                                                                    <  7:282445/146‑6 (MQ=255)
                 gAAGGCTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCAGCTCGCCCAGCGAGcggctgcccaa                                                                                                                                    >  8:282445/1‑141 (MQ=255)
                 gAAGGCTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACATTGGGTATTCGCCCGCAGCGTGGACAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGa                                                                                                                                 >  6:59669/1‑149 (MQ=255)
                  aaGGCTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGAc                                                                                                                                >  7:293814/1‑149 (MQ=255)
                     gCTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAg                                                                                                                             <  4:324186/149‑1 (MQ=255)
                      cTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAgg                                                                                                                            >  3:317472/1‑149 (MQ=255)
                      cTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAgg                                                                                                                            <  8:137422/149‑1 (MQ=255)
                         ggggcgcGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGc                                                                                                                           <  5:268761/144‑1 (MQ=255)
                         ggTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGc                                                                                                                           <  3:282337/147‑1 (MQ=255)
                         ggTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGc                                                                                                                           <  7:184272/147‑1 (MQ=255)
                         ggTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGt                                                                                                                         <  5:83762/149‑1 (MQ=255)
                                     gtATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATc                                                                                                             >  1:89290/1‑149 (MQ=255)
                                         gtggggggCCATGGCCAACCTTGGGTATTCGGCCGCAGCGTGGCCAGCCGCAGGGCCCTGGAACTGGGTGGCGTGGACGTGGAAACCTTGGGCTCCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCAc                                                                                                          <  1:284660/142‑1 (MQ=255)
                                          tGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTc                                                                                                        >  3:145920/1‑149 (MQ=255)
                                            cgcgGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCACTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTc                                                                                                        <  3:156916/147‑1 (MQ=255)
                                             gcgGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCg                                                                                                       >  8:25671/1‑147 (MQ=255)
                                             gcgGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGaa                                                                                                     <  7:115140/149‑1 (MQ=255)
                                                      gCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCAc                                                                                             <  3:9352/148‑1 (MQ=255)
                                                      gCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACt                                                                                            >  1:326289/1‑149 (MQ=255)
                                                       ccAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACt                                                                                            >  2:231121/1‑148 (MQ=255)
                                                       ccAACCTTGGGTATTCGCCCGCAGCGTGGCAATCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCCTCATCCACTCGAAGTGTGCACt                                                                                            >  6:139707/1‑148 (MQ=255)
                                                         aaCCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCa                                                                                              >  7:89083/1‑144 (MQ=255)
                                                         aaCCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCa                                                                                              <  8:89083/144‑1 (MQ=255)
                                                            cTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCca                                                                                      >  8:81592/1‑149 (MQ=255)
                                                             ttGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCca                                                                                      >  2:24637/1‑148 (MQ=255)
                                                               gggTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCGGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCcaca                                                                                    <  6:274255/148‑1 (MQ=255)
                                                                ggTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAg                                                                                   <  2:326289/148‑1 (MQ=255)
                                                                  tATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGcc                                                                                >  7:299539/1‑149 (MQ=255)
                                                                   aTTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCt                                                                               >  2:44975/1‑149 (MQ=255)
                                                                    ttCGCACGCAGCGTGGCCAGCAGCAGCGCACTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTAGCTGGGAGAGCTGCTGTTCTGCGACCAGCCGATCATCCGTCATCCACTCGAAGTGTGAACTTATCCACAGGCCTg                                                                              >  2:168602/1‑149 (MQ=255)
                                                                     tCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTgg                                                                             >  3:39089/1‑149 (MQ=255)
                                                                       gCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGc                                                                            <  5:100688/148‑1 (MQ=255)
                                                                        cccGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTg                                                                          <  8:180415/149‑1 (MQ=255)
                                                                               ggTGGCCAGCCGCTGCGCCCTGGAACGTGGGGGGCTGGACCTGGAAAGCTTGGGCAGGCGCTCGCTGGGCGAGCTGCTGTTCTGCTTCCAGGCGATCATCCGTCATCCACTCGACGTGGGCACTTATCCACAGGCCTGGCTGCCGTCCg                                                                   <  5:305740/148‑1 (MQ=255)
                                                                               cGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCg                                                                   <  3:252787/149‑1 (MQ=255)
                                                                                gTGGCCAGCCGCAGCGCCCTGGAACGTGGGGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCg                                                                   >  1:28398/1‑148 (MQ=255)
                                                                                 tGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGaa                                                                 <  1:273933/149‑1 (MQ=255)
                                                                                   gCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCGCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGc                                                               <  8:311805/149‑1 (MQ=255)
                                                                                     cAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTg                                                             >  8:187422/1‑149 (MQ=255)
                                                                                     cAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTg                                                             <  8:326327/149‑1 (MQ=255)
                                                                                      aGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTg                                                             >  5:174300/1‑148 (MQ=255)
                                                                                      aGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGc                                                            <  2:89290/149‑1 (MQ=255)
                                                                                       gCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGca                                                           <  8:293814/149‑1 (MQ=255)
                                                                                        ccGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGcac                                                          <  8:299539/149‑1 (MQ=255)
                                                                                             gcgcCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATgcg                                                     >  6:328167/1‑149 (MQ=255)
                                                                                              cgcCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCTAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATgcgg                                                    <  3:283812/149‑1 (MQ=255)
                                                                                                     gAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCttt                                              <  7:25671/148‑1 (MQ=255)
                                                                                                      aaCGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTgg                                            <  7:276189/149‑1 (MQ=255)
                                                                                                       accgggtggCCGGGACCTGCAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGGTCTGCGACCAGGCGTTCAGCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTggg                                           <  4:39046/145‑1 (MQ=255)
                                                                                                            gtggCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTGGGgccg                                       <  1:44975/148‑1 (MQ=255)
                                                                                                               gCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGGGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTGGGGCcgccgc                                   >  2:232590/1‑149 (MQ=255)
                                                                                                                    gacctggaAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTggg                                           >  1:376291/1‑136 (MQ=255)
                                                                                                                    gacctggaAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTggg                                           <  2:376291/136‑1 (MQ=255)
                                                                                                                    gacctggaAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTGGGGCCGCCGCTcgc                               <  4:31234/148‑1 (MQ=255)
                                                                                                                     acctggaAAGCTTGGGCAGCGGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTGGGGCCGCCGCTCgcgc                             <  1:168602/149‑1 (MQ=255)
                                                                                                                         ggaAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTGGGGCCGCCGCTCgcgc                             <  1:353187/145‑1 (MQ=255)
                                                                                                                         ggaAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTGGGGCCGCCGCTCgcgc                             >  2:353187/1‑145 (MQ=255)
                                                                                                                         ggaAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTACGAAGCTGCACATGCGGCGCTTTGGGGCCGCCGCTCGCGCTTc                          >  5:262780/1‑148 (MQ=255)
                                                                                                                          gaAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTGGGGCCGCCGCTCGCGCTTc                          <  4:317472/147‑1 (MQ=255)
                                                                                                                               gcagaacagcagCTCGCCCAGCGAGCGGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTGGGGCCGCCGCTCGCGCTTCGAgcgc                    >  8:259214/7‑148 (MQ=255)
                                                                                                                                   aacAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACAGGCGGCGCTTTGGGGCCGCCGCTCGCGCTTCGAGCGCAACgg               >  5:322225/3‑149 (MQ=255)
                                                                                                                                        ccGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTGGGGCCGCCGCTCGCGCTTCGAGCGCAACGGCCTgg          <  5:59669/149‑1 (MQ=255)
                                                                                                                                         cgctcgctGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTGGGGCCGCCGCTCGCGCTTCGAGCGCAACGGCCTGGa         >  8:248599/1‑149 (MQ=255)
                                                                                                                                          gctcgctGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCTCTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACAAGCGGCGCTTTGGGGCCGCCGCTCGCGCTTCGAGCGCAACGGCCTGGAc        <  4:39089/149‑1 (MQ=255)
                                                                                                                                              gctGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTGGGGCCGCCGCTCGCGCTTCGAGCGCAACGGCCTGGACctgc    <  6:174300/149‑1 (MQ=255)
                                                                                                                                                tGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTGGGGCCGCCGCTCGCGCTTCGAGCGCAACGGCCTGGACctgctg  >  1:331091/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCATCGCTGCCGGCGCCGAAGGCTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTGGGGCCGCCGCTCGCGCTTCGAGCGCAACGGCCTGGACCTGCTG  >  NC_002947/6063320‑6063612

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: