Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,263,486 T→G 5.3% L79R (CTG→CGG)  trpF → N‑(5'‑phosphoribosyl)anthranilate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,263,4860TG 5.3% 218.0 / 4.1 75L79R (CTG→CGG) trpFN‑(5'‑phosphoribosyl)anthranilate isomerase
Reads supporting (aligned to +/- strand):  ref base T (32/39);  new base G (2/2);  total (34/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 = 4.83e-01

GTCTTCTATGCCAAGAGCCCGCGGGCGGTGGACGTGCGCCAGGCGCGGGCGATCATTGCCGAGTTGCCGCCCTTCGTGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTCGGCTTTACGCGGGTGCCCGCGGCATTTTGCTC  >  NC_002947/2263341‑2263631
                                                                                                                                                 |                                                                                                                                                 
gTCTTCTATGCCAAGAGCCCGCGGGCGGTGGACGTGCGCCAGGCGCGGGCGATCATTGCCGAGTTGCCGCCCTTCGTGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGAc                                                                                                                                               <  2:1163752/150‑1 (MQ=255)
 tcttctATGCCAAGAGCCCGCGGGCGGTGGACGTGCGCCAGGCGCGGGCGATCATTGCCGAGTTGCCGCCCTTCGTGACAACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGAcc                                                                                                                                              >  2:164499/1‑150 (MQ=255)
    tctATGCCAAGAGCCCGCGGGCGGTGGACGTGCGCCAGGCGCGGGCGATCATTGCCGAGTTGCCGCCCTTCGTGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGc                                                                                                                                           <  2:50150/150‑1 (MQ=255)
          ccAAGAGCCCGCGGGCGGTGGACGTGCGCCAGGCGCGGGCGATCATTGCCGAGTTGCCGCCCTTCGTGACCACGGTCGGGCTGTTCGTAACTGCATCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCACGCTGGACCTGCTACAGt                                                                                                                                     >  2:465686/1‑150 (MQ=255)
           cAAGAGCCCGCGGGCGGTGGACGTGCGCCAGGCGCGGGCGATCATTGCCGAGTTGCCGCCCTTCGTGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGtt                                                                                                                                    <  1:497133/150‑1 (MQ=255)
             agagCCCGCGGGCGGTGGACGTGCGCCAGGCGCGGGCGATCATTGCCGAGTTGCCGCCCTTCGTGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTc                                                                                                                                  <  2:33162/150‑1 (MQ=255)
              gagCCCGCGGGCGGTGGACGTGCGCCAGGCGCGGGCGATCATTGCCGAGTTGCCGCCCTTCGTGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCa                                                                                                                                 <  2:896154/150‑1 (MQ=255)
              gagCCCGCGGGCGGTGGACGTGCGCCAGGCGCGGGCGATCATTGCCGAGTTGCCGCCCTTCGTGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCCGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCa                                                                                                                                 <  2:984053/150‑1 (MQ=255)
                  ccGCGGGCGGTGGACGTGCGCCAGGCGCGGGCGATCATTGCCGAGTTGCCGCCCTTCGTGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGc                                                                                                                             >  2:533885/1‑150 (MQ=255)
                  ccGCGGGCGGTGGACGTGCGCCAGGCGCGGGCGATCATTGCCGAGTTGCCGCCCTTCGTGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGc                                                                                                                             <  2:321184/150‑1 (MQ=255)
                     cgggcggTGGACGTGCGCCAGGCGCGGGCGATCATTGCCGAGTTGCCGCCCTTCGTGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGAt                                                                                                                          <  2:845219/150‑1 (MQ=255)
                             ggACGTGCGCCAGGCGCGGGCGATCATTGCCGAGTTGCCGCCCTCCGTGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAAcccc                                                                                                                  <  1:758903/150‑1 (MQ=255)
                              gACGTGCGCCAGGCGCGGGCGCTCATTGCCGAGTTGCCGCCCTTCGTGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCg                                                                                                                 <  1:642805/150‑1 (MQ=255)
                              gACGTGCGCCAGGCGCGGGCGATCATTGCCGAGTTGCCGCCCTTCGTGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCg                                                                                                                 >  2:1412138/1‑150 (MQ=255)
                                cGTGCGCCAGGCGCGGGCGATCATTGCCGAGTTGCCGCCCTTCGTGAGCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCa                                                                                                               <  1:662326/150‑1 (MQ=255)
                                      ccAGGCGCGGGCGATCATTGCCGAGTTGCCGCCCTTCGTGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACtg                                                                                                         >  2:960263/1‑150 (MQ=255)
                                          gcgcGGGCGATCATTGCCGAGTTGCCTCTTGTCGTGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAACCCCCGCAGGACTGTGaa                                                                                                     <  1:930808/150‑1 (MQ=255)
                                             cgGGCGATCATTGCCGAGTTGCCGCCCTTCGTGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGc                                                                                                  <  1:1412138/150‑1 (MQ=255)
                                                  gATCATTGCCGAGTTGCCGCCCTTCGTGACCACGGTCGGGCTGTTCGTCAATGTCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGGTacca                                                                                             >  2:909489/1‑150 (MQ=255)
                                                          ccGAGTTGCCGCCCTTCGTGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCt                                                                                     <  1:960263/150‑1 (MQ=255)
                                                          ccGAGTTGCCGCCCTTCGTGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCt                                                                                     >  2:273053/1‑150 (MQ=255)
                                                                  ccgccCTTCGTGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAg                                                                             <  1:57726/150‑1 (MQ=255)
                                                                   cgccCTTCGTGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTgg                                                                                                                                                 <  2:360103/81‑1 (MQ=255)
                                                                   cgccCTTCGTGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTgg                                                                                                                                                 >  1:360103/1‑81 (MQ=255)
                                                                   cgccCTTCGTGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAgg                                                                            >  1:987633/1‑150 (MQ=255)
                                                                   cgccCTCCGTGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAgg                                                                            >  2:1138550/1‑150 (MQ=255)
                                                                     cccTTCGTGATTACGGTTGGGCTGTTCGTCAATGTCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGcc                                                                          <  2:181259/150‑1 (MQ=255)
                                                                       cTTCGTGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCAGCGGGAc                                                                                                                                               <  1:1340381/79‑1 (MQ=255)
                                                                       cTTCGTGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCAGCGGGAc                                                                                                                                               >  2:1340381/1‑79 (MQ=255)
                                                                           ttgcccACGGTGGGGCTGTGCCTCAATGTCTCGCGCTGCCAGCTGAACGAGATCCTCAAACTGGTCCCGCTGGACCTGCTACAGTTACCCCGCGACCGATCCCCGCAAGACCGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCgt                                                                    <  2:905073/146‑1 (MQ=255)
                                                                            tGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCgtg                                                                   <  1:328147/150‑1 (MQ=255)
                                                                                cACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTacca                                                                                             >  1:1232044/1‑120 (MQ=255)
                                                                                cACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTacca                                                                                             <  2:1232044/120‑1 (MQ=255)
                                                                                cACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTgcg                                                               >  1:493995/1‑150 (MQ=255)
                                                                                 aCGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAgg                                                                            >  1:864833/1‑136 (MQ=255)
                                                                                 aCGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAgg                                                                            <  2:864833/136‑1 (MQ=255)
                                                                                         gCTGTTCGTCAATGCCTCGCGCTGCTAGCTGGACTATATCCTCGAACTGGCCCCGCTGGACCTGAGCGGGTTGCTCGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGa                                                      <  1:386294/150‑1 (MQ=255)
                                                                                         gCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTTGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGa                                                      >  1:1397364/1‑150 (MQ=255)
                                                                                         gCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGa                                                      >  2:203338/1‑150 (MQ=255)
                                                                                            gTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGa                                                   <  1:1138550/150‑1 (MQ=255)
                                                                                             ttCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGAc                                                  >  2:1109044/1‑150 (MQ=255)
                                                                                             ttCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGAc                                                  >  2:292331/1‑150 (MQ=255)
                                                                                               cGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCt                                                <  1:203338/150‑1 (MQ=255)
                                                                                                gTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTg                                               >  1:713265/1‑150 (MQ=255)
                                                                                                  cAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCCCGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGa                                             <  2:520101/150‑1 (MQ=255)
                                                                                                     tGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAATGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAggc                                          <  2:723140/150‑1 (MQ=255)
                                                                                                        cTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGAt                                                     >  2:487436/1‑136 (MQ=255)
                                                                                                        cTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGAt                                                     <  1:487436/136‑1 (MQ=255)
                                                                                                        cTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAggcggc                                       >  1:529004/1‑150 (MQ=255)
                                                                                                         tCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGcc                                      <  1:1111886/150‑1 (MQ=255)
                                                                                                                   aGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTCGGccag                            <  1:292331/150‑4 (MQ=255)
                                                                                                                    gCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAggcgg                                        <  1:890274/137‑1 (MQ=255)
                                                                                                                    gCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAggcgg                                        >  2:890274/1‑137 (MQ=255)
                                                                                                                    gCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTCGGCTTTa                           >  1:533709/1‑150 (MQ=255)
                                                                                                                       gAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTCGGCTTTAcgc                        <  1:533885/150‑1 (MQ=255)
                                                                                                                        cccGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGc                                                                      <  1:674621/101‑1 (MQ=255)
                                                                                                                        aaCGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGc                                                                      >  2:674621/1‑103 (MQ=255)
                                                                                                                        aaCGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTCGGCTTTAcgcg                       <  2:343353/150‑1 (MQ=255)
                                                                                                                         aCGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTCGGCTTTACGCgg                      <  2:712738/150‑1 (MQ=255)
                                                                                                                             gaTCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTCGGCTTTACGCGGGTGc                  >  1:16371/1‑150 (MQ=255)
                                                                                                                             gaTCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTCGGCTTTACGCGGGTGc                  >  1:485930/1‑150 (MQ=255)
                                                                                                                                 cTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTCGGCTTTACGCGGGTGCCcgc              >  2:439653/1‑150 (MQ=255)
                                                                                                                                  tCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGGCTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTCGGCTTTACGCGGGTGCCcgcg             >  1:152980/1‑150 (MQ=255)
                                                                                                                                  tCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTCGGCTTTACGCGGGTGCCcgcg             >  1:541833/1‑150 (MQ=255)
                                                                                                                                  tCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTCGGCTTTACGCGGGTGCCcgcg             >  2:257321/1‑150 (MQ=255)
                                                                                                                                    gAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTTGAGGCGGCCTGTCGGCTTTACGCGGGTGCCCGCGGc           >  1:1008134/1‑150 (MQ=255)
                                                                                                                                     aaGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTCGGCTTTACGCGGGTGCCCGCGGCa          <  1:465686/150‑1 (MQ=255)
                                                                                                                                      gcaggTCCAGCGGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGAt                                                     <  2:1028824/103‑1 (MQ=255)
                                                                                                                                      gcaggTCCAGCGGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGAt                                                     >  1:1028824/4‑106 (MQ=255)
                                                                                                                                       gTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTCGGCTTTACGCGGGTGCCCGCGGCAtt        <  2:485930/150‑1 (MQ=255)
                                                                                                                                          gTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTCGGCTTTACGCGGGTGCCCGCGGCATTTTg     >  1:1073394/1‑150 (MQ=255)
                                                                                                                                           tCCCGCTGCACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTCGGCTTTACGCGGGTGCCCGCGGCATTTTGc    <  1:359999/150‑1 (MQ=255)
                                                                                                                                            cccGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTCGGCTTTACGCGGGTGCCCGCGGCATTTTGCt   >  1:244884/1‑150 (MQ=255)
                                                                                                                                             ccGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTCGGCTTTACGCGGGTGCCCGCGGCATTTTGCTc  <  2:725211/150‑1 (MQ=255)
                                                                                                                                             ccGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTCGGCTTTACGCGGGTGCCCGCGGCATTTTGCTc  <  1:164499/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GTCTTCTATGCCAAGAGCCCGCGGGCGGTGGACGTGCGCCAGGCGCGGGCGATCATTGCCGAGTTGCCGCCCTTCGTGACCACGGTCGGGCTGTTCGTCAATGCCTCGCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTCGGCTTTACGCGGGTGCCCGCGGCATTTTGCTC  >  NC_002947/2263341‑2263631

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