Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 264,132 A→C 12.1% I295L (ATC→CTC)  gabD‑I → succinate‑semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947264,1320AC12.1% 556.1 / 54.4 195I295L (ATC→CTC) gabD‑Isuccinate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (84/87);  new base C (12/12);  total (96/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 201 total aligned reads displayed.

CGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATC  >  NC_002947/263984‑264281
                                                                                                                                                    |                                                                                                                                                     
cgcCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTAt                                                                                                                                                      <  1:5061718/150‑1 (MQ=255)
  cgccgaccATCAAGAAGCTTGGGCTTGAGCTTGGCGGCAACGCGCCTTTCATCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACGGCCAAACCTGCGTCTGCGCCAACCGTctct                                                                                                                                                    <  1:2777033/143‑1 (MQ=255)
  cgccgaccATCAAGAAGCTTGGGCTTGAGCTTGGCGGCAACGCGCCTTTCATCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACGGCCAAACCTGCGTCTGCGCCAACCGTctct                                                                                                                                                    <  1:2777218/143‑1 (MQ=255)
   cAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTa                                                                                                                                                   >  1:4786771/1‑150 (MQ=255)
   cAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTa                                                                                                                                                   >  1:1296172/1‑150 (MQ=255)
       gACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTg                                                                                                                                               <  1:432161/150‑1 (MQ=255)
        aCATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGc                                                                                                                                              >  2:186261/1‑150 (MQ=255)
         cATCAAGAAGCTTGGGCTTGAGCTTGGCGGCAACGCGCCTTTCATCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACa                                                                                                                                             >  1:1110172/1‑150 (MQ=255)
         cATCAAGAAGCTTGGGCTTGAGCTTGGCGGCAACGCGCCTTTCATCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACa                                                                                                                                             >  1:2092476/1‑150 (MQ=255)
         cATCAAGAAGCTTGGGCTTGAGCTTGGCGGCAACGCGCCTTTCATCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACAGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACa                                                                                                                                             <  2:713263/150‑1 (MQ=255)
           tCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAgg                                                                                                                                           >  1:671489/1‑150 (MQ=255)
             aagaagGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGAc                                                                                                                                         >  1:835213/1‑150 (MQ=255)
             aagaagGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGAc                                                                                                                                         >  1:3197777/1‑150 (MQ=255)
              agaagGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACg                                                                                                                                        >  2:1145581/1‑150 (MQ=255)
               gaagGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACgg                                                                                                                                       >  1:2267641/1‑150 (MQ=255)
                aagGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGc                                                                                                                                      <  2:1436971/150‑1 (MQ=255)
                        ccTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTacga                                                                                                                              <  1:5262636/150‑1 (MQ=255)
                        ccTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTacga                                                                                                                              >  1:3823558/1‑150 (MQ=255)
                        ccTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTacga                                                                                                                              <  1:4495837/150‑1 (MQ=255)
                        ccTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTacga                                                                                                                              >  2:4617820/1‑150 (MQ=255)
                            gAGCTTGGCGGCAACGCGCCTTTCATCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACGGCGTCTATGAtgcc                                                                                                                          <  2:886688/150‑2 (MQ=255)
                               cttggcggCAACGCGCCTTTCATCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACGGCGTCTATGATGCCttt                                                                                                                       <  1:5409097/144‑2 (MQ=255)
                                ttggcggCAACGCGCCTTTCATCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACGGCGTCTATGATGCCtttg                                                                                                                      >  1:804242/6‑148 (MQ=255)
                                tGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCg                                                                                                                      >  2:5158694/1‑150 (MQ=255)
                                 gggTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGTCGCGTTCGc                                                                                                                     >  2:3196163/1‑150 (MQ=255)
                                 gggTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGc                                                                                                                     <  1:2024866/150‑1 (MQ=255)
                                  ggcggCATCGCTCCTTTCATCGTGTTCGACTACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACCACGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACGGCGTCTATGATGCCTTTGcc                                                                                                                    <  2:2693764/147‑1 (MQ=255)
                                  ggTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTc                                                                                                                       <  1:1256096/147‑1 (MQ=255)
                                      gCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGTGTTAGCCgaga                                                                                                                >  1:3481564/1‑150 (MQ=255)
                                      gCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCgaga                                                                                                                >  2:2349931/1‑150 (MQ=255)
                                      gCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCgaga                                                                                                                >  1:2286302/1‑150 (MQ=255)
                                       cAACGCGCCTTTCATCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACGGCGTCTATGATGCCTTTGCCGAGaa                                                                                                               >  1:1152410/1‑150 (MQ=255)
                                       cAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGaa                                                                                                               <  2:2914906/150‑1 (MQ=255)
                                           gCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTg                                                                                                           <  2:4913864/150‑1 (MQ=255)
                                            cgcctttCATCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACGGCGTCTATGATGCCTTTGCCGAGAAGCtct                                                                                                          <  1:3508143/145‑3 (MQ=255)
                                             cccGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTa                                                                                                                                                       <  2:1068363/104‑1 (MQ=255)
                                              ccGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGcc                                                                                                        >  2:3521307/1‑150 (MQ=255)
                                                    aTCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACg                                                                                                                                        <  2:3557660/112‑1 (MQ=255)
                                                     tCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAg                                                                                                 >  2:1370126/1‑150 (MQ=255)
                                                         gtTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGc                                                                                             <  1:5398953/150‑1 (MQ=255)
                                                          ttCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGc                                                                                                             >  2:4216746/1‑133 (MQ=255)
                                                          ttCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGc                                                                                                             <  1:4216746/133‑1 (MQ=255)
                                                           tCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCa                                                                                           <  2:4850149/150‑1 (MQ=255)
                                                            cgacgaCGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCaa                                                                                          <  2:2261091/150‑1 (MQ=255)
                                                             gacgacGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAg                                                                                         >  2:5414162/1‑150 (MQ=255)
                                                              acgacgCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGc                                                                                        >  1:4475376/1‑150 (MQ=255)
                                                              acgacgCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGc                                                                                        >  1:1731880/1‑150 (MQ=255)
                                                                 acgCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGa                                                                                     <  2:4772501/150‑1 (MQ=255)
                                                                  cgCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGaa                                                                                    >  2:5195074/1‑150 (MQ=255)
                                                                    ccGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGa                                                                                  >  2:4714243/1‑150 (MQ=255)
                                                                     cGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGAt                                                                                 <  1:4871576/150‑1 (MQ=255)
                                                                      gACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATc                                                                                <  1:4742049/150‑1 (MQ=255)
                                                                         cTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGt                                                                             >  2:2000463/1‑150 (MQ=255)
                                                                               aaGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGcc                                                                      >  2:1638463/1‑150 (MQ=255)
                                                                                     gTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAATATCGGTAACGGCCTGgaa                                                                 >  1:1912917/1‑150 (MQ=255)
                                                                                     gTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGgaa                                                                 <  2:4347433/150‑1 (MQ=255)
                                                                                     gTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCGACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGGTGCCAAGCTGAAGATCGGTAACGGCCTGgaa                                                                 <  2:1090264/150‑1 (MQ=255)
                                                                                        gAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGgaagaa                                                              <  1:4089191/150‑1 (MQ=255)
                                                                                        gAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGgaagaa                                                              <  2:398575/150‑1 (MQ=255)
                                                                                         aGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGaagaag                                                             <  1:1037872/150‑1 (MQ=255)
                                                                                              gcgATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGcacc                                                        <  1:5275930/150‑1 (MQ=255)
                                                                                                gATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGc                                                                                                                                      >  2:3734213/1‑70 (MQ=255)
                                                                                                gATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGc                                                                                                                                      <  1:3734213/70‑1 (MQ=255)
                                                                                                gATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGcaccac                                                      >  1:893669/1‑150 (MQ=255)
                                                                                                    atcTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACt                                                  >  1:3382470/1‑150 (MQ=255)
                                                                                                     tctcCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCg                                                                               <  1:4936304/120‑1 (MQ=255)
                                                                                                     tctcCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTg                                                 <  1:185871/150‑1 (MQ=255)
                                                                                                     tctcCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTg                                                 <  2:3427251/150‑1 (MQ=255)
                                                                                                     tctcCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTg                                                 <  1:2977812/150‑1 (MQ=255)
                                                                                                      ctcCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTgg                                                <  2:1826817/150‑1 (MQ=255)
                                                                                                       tcCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGc                                               <  2:1279991/150‑1 (MQ=255)
                                                                                                        ccAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGcc                                              <  1:2349931/150‑1 (MQ=255)
                                                                                                         cAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGccc                                             >  1:1337453/1‑150 (MQ=255)
                                                                                                         cAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGccc                                             >  2:3884643/1‑150 (MQ=255)
                                                                                                               ccGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGAt                                       >  2:2806268/1‑150 (MQ=255)
                                                                                                               ccGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGAt                                       >  1:3258999/1‑150 (MQ=255)
                                                                                                                 gTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCg                                     >  1:1073607/1‑150 (MQ=255)
                                                                                                                      aacGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGc                                >  2:3618814/1‑150 (MQ=255)
                                                                                                                        cGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCaa                              <  1:5158694/150‑1 (MQ=255)
                                                                                                                        cGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACGGCGTCTATGATGCCTTTGCCGAGAAGCTCTCCAAGGCGGTCGGCTCGCTCAAGACCGGCAACGGCTTCGACGAGGGCATCAATCTCGGCCCGCTGATCGATgaggc                              <  2:1421100/150‑5 (MQ=255)
                                                                                                                         ggCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAg                             <  1:4392832/150‑1 (MQ=255)
                                                                                                                          gCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAgg                            <  2:1337453/150‑1 (MQ=255)
                                                                                                                          gCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAgg                            <  1:38943/150‑1 (MQ=255)
                                                                                                                           ccAAACCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACGGCGTCTATGATGCCTTTGCCGAGAAGCTCTCCAAGGCGGTCGGCTCGCTCAAGACCGGCAACGGCTTCGACGAGGGCATCAATCTCGGCCCGCTGATCGATgaggcggc                           >  2:5344036/1‑143 (MQ=255)
                                                                                                                             aaaCCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACGGCGTCTATGATGCCTTTGCCGAGAAGCTCTCCAAGGCGGTCGGCTCGCTCAAGACCGGCAACGGCTTCGACGAGGGCATCAATCTCGGCCCGCTGATCGAt                                   >  1:1728743/3‑140 (MQ=255)
                                                                                                                             aaaCCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACGGCGTCTATGATGCCTTTGCCGAGAAGCTCTCCAAGGCGGTCGGCTCGCTCAAGACCGGCAACGGCTTCGACGAGGGCATCAATCTCGGCCCGCTGATCGAt                                   <  2:1728743/138‑1 (MQ=255)
                                                                                                                               aCCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTc                       <  1:100638/150‑1 (MQ=255)
                                                                                                                                    cGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCaa                  <  2:2052466/150‑1 (MQ=255)
                                                                                                                                    cGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCaa                  <  1:5052712/150‑1 (MQ=255)
                                                                                                                                     gTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAg                 <  1:1390125/150‑1 (MQ=255)
                                                                                                                                       cTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGt               <  1:754439/150‑1 (MQ=255)
                                                                                                                                          cgcCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCa            <  1:3581861/150‑1 (MQ=255)
                                                                                                                                           gcCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAg           >  2:4588060/1‑150 (MQ=255)
                                                                                                                                            ccAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATg                                  <  1:2536111/126‑1 (MQ=255)
                                                                                                                                            ccAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATg                                  >  2:2536111/1‑126 (MQ=255)
                                                                                                                                            ccAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAgg          >  1:4616075/1‑150 (MQ=255)
                                                                                                                                              aaCCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGcc                                                                                                        >  1:4281649/1‑54 (MQ=255)
                                                                                                                                              aaCCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGcc                                                                                                        <  2:4281649/54‑1 (MQ=255)
                                                                                                                                                ccGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAaca      >  1:2705261/1‑150 (MQ=255)
                                                                                                                                                ccGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAaca      >  1:355824/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATC  >  NC_002947/263984‑264281

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: