Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,502,665 C→G 5.5% W607C (TGG→TGC fryA ← putative PTS enzyme: Hpr, enzyme I and IIA components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,502,6650CG 5.5% 332.4 / 7.0 104W607C (TGG→TGCfryAputative PTS enzyme: Hpr, enzyme I and IIA components
Reads supporting (aligned to +/- strand):  ref base C (51/47);  new base G (4/2);  total (55/49)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-01

Only 100 of 107 total aligned reads displayed.

AGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACT  >  NC_000913/2502516‑2502811
                                                                                                                                                     |                                                                                                                                                  
aGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTAcc                                                                                                                                                    >  1:333300/1‑150 (MQ=255)
aGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTAcc                                                                                                                                                    >  2:1116817/1‑150 (MQ=255)
aGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTAcc                                                                                                                                                    >  1:1450600/1‑150 (MQ=255)
     ccGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCAACCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTATCACCGCCCAGAAACCCGCAAATCCATACCCAtat                                                                                                                                               >  2:891166/1‑148 (MQ=255)
        acacGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGcc                                                                                                                                            <  2:875150/150‑1 (MQ=255)
         cacGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGGGTAGGCATTTGcgg                                                                                                                                           <  2:235191/150‑3 (MQ=255)
           cGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCg                                                                                                                                         >  2:38478/1‑150 (MQ=255)
           cGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCg                                                                                                                                         >  1:135916/1‑150 (MQ=255)
            gCCTCGCTATCCAGTTGTCGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCGACTCATAGTCTGCTCATCCAGGCCAAGCCCAAGCAGTAGCGGCAGATAACGGCTATCACCGCCAAGGTCACCGCAAAAGCCGAGCCATTAGCCCCGc                                                                                                                                        >  1:295960/1‑150 (MQ=255)
             ccTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCt                                                                                                                                       <  2:42315/150‑1 (MQ=255)
              cTCGCTATCCAGTTGAAGAATCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCTGCTGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTg                                                                                                                                      <  1:1440445/150‑1 (MQ=255)
                 gCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATg                                                                                                                                   <  2:1427923/150‑1 (MQ=255)
                 gCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATg                                                                                                                                   <  2:586336/150‑1 (MQ=255)
                     tCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCa                                                                                                                               >  1:48508/1‑150 (MQ=255)
                       cAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGt                                                                                                                             >  1:1265968/1‑150 (MQ=255)
                         gTTGACGAATCTGGCTTTTCACCGCCGGAATAGGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCATTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCAATTTGCCCCGCTGATGCGCAGTgg                                                                                                                           <  2:321248/150‑1 (MQ=255)
                         gTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCCCCGCCCAGTTCCCCGCAATTGCCTACCCAGGTGCCCCGCTGAGGCGCTGTgg                                                                                                                           >  2:196470/1‑150 (MQ=255)
                         gTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTgg                                                                                                                           <  1:1116817/150‑1 (MQ=255)
                           tGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATGCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCGTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTa                                                                                                                         <  1:832049/150‑1 (MQ=255)
                            gACGCAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCCCCGCAAATGCCTACCCATTTGCCCCGCAGATGCGCAGTGGTaa                                                                                                                        >  2:1113892/1‑150 (MQ=255)
                            gACGAAGCTGGCTTTTCCCCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTaa                                                                                                                        >  2:928890/1‑150 (MQ=255)
                              cGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACt                                                                                                                      >  1:358723/1‑150 (MQ=255)
                              cGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACt                                                                                                                      >  1:1676709/1‑150 (MQ=255)
                                  gCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTAttt                                                                                                                  <  2:64259/150‑1 (MQ=255)
                                   cTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTtg                                                                                                                 <  1:33181/150‑1 (MQ=255)
                                   cTGGCTTTTAACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCAGCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTtg                                                                                                                 <  2:1489090/150‑1 (MQ=255)
                                    tGGCTTTTCTCCGCCGGAATTCGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCATGCAGTAGCGGCAGATTACGGCTTTCACCGCCCAGTTCAGCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTtgc                                                                                                                <  1:891166/150‑1 (MQ=255)
                                     ggCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTtgct                                                                                                               >  2:1356731/1‑150 (MQ=255)
                                      gCTTTTCACCGCCGGAATAAGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTtgctg                                                                                                              <  1:254869/150‑1 (MQ=255)
                                       cTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCCCCGCAAATCCCTACCCATTTTACACGGTGATGAGCAGTTGTAACTATTtggtgc                                                                                                             >  1:256541/1‑150 (MQ=255)
                                            cACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCGACCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACAt                                                                                                        >  1:54029/1‑150 (MQ=255)
                                             aCCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATg                                                                                                       >  1:955974/1‑150 (MQ=255)
                                                 ccgGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCa                                                                                                   <  2:37640/150‑1 (MQ=255)
                                                   ggAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAgg                                                                                                 >  1:1414453/1‑150 (MQ=255)
                                                   ggAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAgg                                                                                                 <  2:407542/150‑1 (MQ=255)
                                                     aaTACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGaa                                                                                               >  2:44835/1‑150 (MQ=255)
                                                       tACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACg                                                                                             <  1:38478/150‑1 (MQ=255)
                                                        aCGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTAACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGa                                                                                            >  2:533530/1‑150 (MQ=255)
                                                            gggCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCACGAACGATGGc                                                                                        >  2:745590/1‑150 (MQ=255)
                                                             ggCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCCCCGCAACTGCCTACCCATTTGCCCCGCTGATGCGCAGTGGATCCTATTTGCTGCAACATGCGCAGGACCGATGGCg                                                                                       >  1:434560/1‑150 (MQ=255)
                                                             ggCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCg                                                                                       >  2:792848/1‑150 (MQ=255)
                                                                  ctcatactcaGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATc                                                                                  >  1:903451/1‑150 (MQ=255)
                                                                    catactcaGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGGAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCgg                                                                                >  1:335910/1‑150 (MQ=255)
                                                                     atactcaGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCggg                                                                               <  2:1239587/150‑1 (MQ=255)
                                                                     atactcaGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCggg                                                                               <  2:1345945/150‑1 (MQ=255)
                                                                        ctcaGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATTACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTta                                                                            >  1:354578/1‑150 (MQ=255)
                                                                        ctcaGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAAATATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTta                                                                            >  2:1694439/1‑150 (MQ=255)
                                                                           aGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTtatat                                                                         >  1:1492381/1‑150 (MQ=255)
                                                                            gCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTtatata                                                                        >  1:580957/1‑150 (MQ=255)
                                                                               cATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCggg                                                                               >  1:1224531/1‑140 (MQ=255)
                                                                               cATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCggg                                                                               <  2:1224531/140‑1 (MQ=255)
                                                                                  ccAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAgcggcg                                                                  <  1:299338/150‑1 (MQ=255)
                                                                                   cAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGa                                                                 >  1:580138/1‑150 (MQ=255)
                                                                                   cAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGa                                                                 >  1:381612/1‑150 (MQ=255)
                                                                                     ggCCCAGCTCAAGCAGTAGCGGCTGATTTCGGCTTTCACCGCCCAGTTCACCGCAAATTCCTAACCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATa                                                               <  2:256541/150‑1 (MQ=255)
                                                                                     ggCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAACTGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATa                                                               <  2:775229/150‑1 (MQ=255)
                                                                                     ggCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATa                                                               >  1:1325081/1‑150 (MQ=255)
                                                                                     ggCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATa                                                               <  2:768844/150‑1 (MQ=255)
                                                                                     ggCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATa                                                               <  2:358723/150‑1 (MQ=255)
                                                                                         cagcccaAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATAcgcg                                                           >  1:549603/1‑150 (MQ=255)
                                                                                            cccaAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCgg                                                        >  1:978106/1‑150 (MQ=255)
                                                                                              caAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGt                                                      <  2:1676709/150‑1 (MQ=255)
                                                                                              caAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGt                                                      <  2:938034/150‑1 (MQ=255)
                                                                                                    gTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTAc                                                <  2:586470/150‑1 (MQ=255)
                                                                                                      aGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGa                                              >  2:1007056/1‑150 (MQ=255)
                                                                                                      aGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGCTTATATAGCGGCGATACGCGCGGGTTATTACGa                                              >  1:1378265/1‑150 (MQ=255)
                                                                                                        cGGCAGATAACGGCTTTCTCCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATc                                            >  2:926940/1‑150 (MQ=255)
                                                                                                        cGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATc                                            <  2:1450600/150‑1 (MQ=255)
                                                                                                         ggCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAAAGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCCACCGGCGCAGGACCGCTGGCGTATTCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCg                                           >  2:954676/1‑150 (MQ=255)
                                                                                                          gCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCgc                                          >  2:938462/1‑149 (MQ=255)
                                                                                                            aGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGAcc                                        <  2:48508/150‑1 (MQ=255)
                                                                                                             gATAACGGCTTTCACGGCCCAGTTCACCGCAAATGCCTACCCATTTGGCCCGCTGATGCGCAGTGGTATCTATTTGCTGCAACATGCGCTGGAACGATGGCGTTTTCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCg                                       <  1:891148/150‑1 (MQ=255)
                                                                                                             gATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCg                                       >  2:554148/1‑150 (MQ=255)
                                                                                                                 aCGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATa                                   >  2:278805/1‑150 (MQ=255)
                                                                                                                   ggCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACa                                 <  2:1414453/150‑1 (MQ=255)
                                                                                                                         cACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACTATGGCGTAATCGGGTTATATAGCGGCGATACGCGCTGGTTATTACGATCGACCGCATACAGATACt                           <  2:861610/150‑1 (MQ=255)
                                                                                                                           ccGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTgg                         <  1:928890/150‑1 (MQ=255)
                                                                                                                             gCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGt                       <  2:381612/150‑1 (MQ=255)
                                                                                                                             gCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGt                       >  1:669849/1‑150 (MQ=255)
                                                                                                                             gCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGt                       <  2:580138/150‑1 (MQ=255)
                                                                                                                                  gTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCatat                  <  2:580957/150‑1 (MQ=255)
                                                                                                                                   agcggggCAAATGGGTAGGCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATc                 >  1:14383/7‑150 (MQ=255)
                                                                                                                                     cggggCAAATGGGTAGGCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGt               >  2:482350/5‑150 (MQ=255)
                                                                                                                                      ggggCAAATGGGTAGGCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTACAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGtt              >  1:1292876/4‑150 (MQ=255)
                                                                                                                                      aCCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGtt              <  1:792848/150‑1 (MQ=255)
                                                                                                                                       ccGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATTGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTg             <  2:54029/150‑1 (MQ=255)
                                                                                                                                       ccGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTg             >  1:1233121/1‑150 (MQ=255)
                                                                                                                                       ccGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTg             <  2:135916/150‑1 (MQ=255)
                                                                                                                                        ggCAAATGGGTAGGCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTgg            >  2:321806/2‑150 (MQ=255)
                                                                                                                                        cGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTgg            >  2:859373/1‑150 (MQ=255)
                                                                                                                                          cAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAg          <  2:1265968/150‑1 (MQ=255)
                                                                                                                                          cAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAg          >  1:1430660/1‑150 (MQ=255)
                                                                                                                                          cAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAg          >  2:1322363/1‑150 (MQ=255)
                                                                                                                                           aactGCCTACGCATTTGCCCCGCGGATGCGCAGTGGTAACTTTTTGCTGGAACATGCGCAGGAACGATTGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTGATTACGTTCGACCGCATTCAGATACTGGGTCATATCGTTGGAGc         <  2:959119/147‑1 (MQ=255)
                                                                                                                                            aaTGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGcc        <  2:571755/150‑1 (MQ=255)
                                                                                                                                              tGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGa      >  2:517422/1‑150 (MQ=255)
                                                                                                                                               gCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATATAGATACTGGGTCATATCGTTGGAGCCGAt     >  2:1348704/1‑150 (MQ=255)
                                                                                                                                               gCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGAt     >  2:513872/1‑150 (MQ=255)
                                                                                                                                                ccTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATa    >  2:1535475/1‑150 (MQ=255)
                                                                                                                                                  tACCCATTTGCCCCGCTGATGCGCATTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGAAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGTGTCATATCGTTGGAGCCGATACt  <  1:1067383/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACT  >  NC_000913/2502516‑2502811

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: