Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,503,974 C→G 8.6% intergenic (+22/‑78) frlD → / → frlR fructoselysine 6‑kinase/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,503,9740CG 8.6% 264.0 / 14.4 81intergenic (+22/‑78)frlD/frlRfructoselysine 6‑kinase/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (33/40);  new base G (3/4);  total (37/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GGGTGCCGGAGATTCGTTCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCAACTCCTCTACGCTACCGTCCGCCAGCGACTGCTGGATGATAT  >  NC_000913/3503826‑3504122
                                                                                                                                                    |                                                                                                                                                    
gggTGCCGGAGATTCGTTCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATc                                                                                                                                                      >  2:19181/1‑149 (MQ=255)
 ggTGCCGGAGATTCGTTCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATc                                                                                                                                                      <  3:449926/148‑1 (MQ=255)
  gTGCCGGAGATTCGTTCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCtt                                                                                                                                                    >  4:280483/1‑149 (MQ=255)
    gCCGGAGATTCGTTCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCa                                                                                                                                                  >  8:105638/1‑149 (MQ=255)
     ccGGAGATTCGTTCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATCGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCAc                                                                                                                                                 <  3:187352/149‑1 (MQ=255)
      cGGAGATTCGTTCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCAc                                                                                                                                                 <  8:66386/148‑1 (MQ=255)
      cGGAGATTCGTTCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCAc                                                                                                                                                 >  7:66386/1‑148 (MQ=255)
          gaTTCGTTCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAAGATGCTCACGCCaa                                                                                                                                            <  6:213476/149‑1 (MQ=255)
          gaTTCGTTCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAAGATGCTCACGCCaa                                                                                                                                            <  1:280309/149‑1 (MQ=255)
          gaTTCGGTCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAAGATGCTCACGCCaa                                                                                                                                            <  8:463516/149‑1 (MQ=255)
           aTTCGTTCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAAc                                                                                                                                           >  3:387566/1‑149 (MQ=255)
              cGTTCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGt                                                                                                                                         <  6:401321/148‑1 (MQ=255)
                ttCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCt                                                                                                                                      <  3:153081/149‑1 (MQ=255)
                   aTTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGt                                                                                                                                    <  1:134068/148‑1 (MQ=255)
                       ccGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACtt                                                                                                                               <  2:87072/149‑1 (MQ=255)
                         ggATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGc                                                                                                                             >  3:54614/1‑149 (MQ=255)
                          gATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGcc                                                                                                                            >  7:257090/1‑149 (MQ=255)
                                tttGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGaaaa                                                                                                                      >  2:68381/1‑149 (MQ=255)
                                  tGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGaaaa                                                                                                                      <  5:89939/147‑1 (MQ=255)
                                  tGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGaaaa                                                                                                                      <  6:206098/147‑1 (MQ=255)
                                   gCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGa                                                                                                                   >  4:36273/1‑149 (MQ=255)
                                    cGGCTGGTCTGCGGGGATGCCATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGa                                                                                                                   <  6:324439/148‑1 (MQ=255)
                                     ggCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGAcc                                                                                                                 <  3:36273/149‑1 (MQ=255)
                                     ggCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGAcc                                                                                                                 <  3:286182/149‑1 (MQ=255)
                                       cTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGc                                                                                                                                       <  2:66779/125‑1 (MQ=255)
                                       cTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGc                                                                                                                                       >  1:66779/1‑125 (MQ=255)
                                          gTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCtat                                                                                                             >  8:253645/1‑148 (MQ=255)
                                          gTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCta                                                                                                              <  1:68381/147‑1 (MQ=255)
                                            cTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAAt                                                                                                          <  1:124622/149‑1 (MQ=255)
                                              gCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGGAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATcc                                                                                                        >  1:228631/1‑149 (MQ=255)
                                                ggggATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGa                                                                                                      <  3:280483/149‑1 (MQ=255)
                                                  ggATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAg                                                                                                     >  8:25947/1‑148 (MQ=255)
                                                   gATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGt                                                                                                    >  6:279507/1‑148 (MQ=255)
                                                    aTGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGa                                                                                                      <  3:281118/145‑1 (MQ=255)
                                                    aTGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGa                                                                                                      >  4:281118/1‑145 (MQ=255)
                                                    aTGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTaa                                                                                                  >  3:24889/1‑149 (MQ=255)
                                                         aTTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGCCCCTATAATCCGAGTAattcat                                                                                             >  6:457316/1‑149 (MQ=255)
                                                            aCCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCATCGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAattcattc                                                                                           <  5:457316/148‑1 (MQ=255)
                                                             ccGCAGGCGATAGCGCAGGGAACGGCGCGGGCGGCGAAAACCCTTCAGTACCCCGCTGCCTGGTCGGTATAACGTTGGCGTTAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCtt                                                                                         <  3:245277/149‑1 (MQ=255)
                                                              cGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCttt                                                                                        <  4:387566/149‑1 (MQ=255)
                                                              cGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCtt                                                                                         >  1:64928/1‑148 (MQ=255)
                                                               gCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTa                                                                                       >  5:346098/1‑149 (MQ=255)
                                                               gCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTa                                                                                       <  7:105638/149‑1 (MQ=255)
                                                                    cGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCa                                                                                  <  4:24889/149‑1 (MQ=255)
                                                                          cgcAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTc                                                                             >  2:364757/1‑148 (MQ=255)
                                                                             aGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGAtt                                                                         <  7:253645/149‑1 (MQ=255)
                                                                              gggAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTa                                                                        <  5:279507/149‑1 (MQ=255)
                                                                               ggAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGCCCCTATAAACAGAGTAATTCATTCTTTATTTCAGGGTCGATTAt                                                                       >  6:344351/1‑149 (MQ=255)
                                                                                 aaCGGCGCGCGCGGCGAAATCCATTCACTACCACGGTGCCTGGTAGGTATAACGTTGGAGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGt                                                                     <  6:74581/149‑1 (MQ=255)
                                                                                  aCGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGt                                                                     >  7:370761/1‑148 (MQ=255)
                                                                                  aCGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTc                                                                    <  3:320678/149‑1 (MQ=255)
                                                                                   cGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCa                                                                   <  2:107490/149‑1 (MQ=255)
                                                                                    ggCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGAt                                                                          <  4:323711/141‑1 (MQ=255)
                                                                                    ggCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGAt                                                                          >  3:323711/1‑141 (MQ=255)
                                                                                     gCGTGCGCGGCGAAAACCATTCAGTACCACGGAGCCTGCTAGGTATAACGTTGGCGTGTGCATCTTCACGCCAACGTGCTGCTACATGCCGGAAAACGACCCTATATTCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGc                                                                 <  2:357272/149‑1 (MQ=255)
                                                                                                 aaaaCCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTa                                                     >  7:369021/1‑149 (MQ=255)
                                                                                                   aaCCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAattcattc                                                                                           <  7:83723/109‑1 (MQ=255)
                                                                                                   aaCCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAattcattc                                                                                           >  8:83723/1‑109 (MQ=255)
                                                                                                   aaCCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTAc                                                    <  2:228631/148‑1 (MQ=255)
                                                                                                    tccATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTActc                                                  <  3:386994/148‑1 (MQ=255)
                                                                                                      cATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCCACGTGCTGTTACTTGACGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTActct                                                 >  6:437228/1‑148 (MQ=255)
                                                                                                      cATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCta                                                                                                              <  1:196095/87‑1 (MQ=255)
                                                                                                      cATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCta                                                                                                              >  2:196095/1‑87 (MQ=255)
                                                                                                      cATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGa                                                                                                      >  3:245906/1‑95 (MQ=255)
                                                                                                      cATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGa                                                                                                      <  4:245906/95‑1 (MQ=255)
                                                                                                          cAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCtcatc                                             <  1:19181/148‑1 (MQ=255)
                                                                                                           aGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCaa                                           >  7:145544/1‑149 (MQ=255)
                                                                                                           aGTACCAAGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCaa                                           <  5:344351/149‑1 (MQ=255)
                                                                                                            gTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCAAc                                          <  4:334949/149‑1 (MQ=255)
                                                                                                             tACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATATTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAAGGTCGATTATGTCAGCTACGGACCTCTACTCTCATCAAct                                         >  8:214120/1‑149 (MQ=255)
                                                                                                               ccACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCAActc                                        <  8:369021/148‑1 (MQ=255)
                                                                                                                          ggtaggtaTAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCAACTCCTCTACGCTAcc                            >  6:297444/1‑149 (MQ=255)
                                                                                                                               acagcaCGTTGGCGTGAAGATGCTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCAACTCCTCTACGCTACCGTcc                        >  1:263571/6‑148 (MQ=39)
                                                                                                                               acagcaCGTTGGCGTGAAGATGCTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCAACTCCTCTACGCTACCGTcc                        <  2:263571/143‑1 (MQ=39)
                                                                                                                                tataACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCAACTCCTCTACGCTACCGTCCGc                      <  1:364757/149‑1 (MQ=255)
                                                                                                                                tataACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCAACTCCTCTACGCTACCGTCCGc                      >  2:101412/1‑149 (MQ=255)
                                                                                                                                  taACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTGATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCAACTCCTCTACGCTACCGTCCGc                      >  2:41971/1‑147 (MQ=255)
                                                                                                                                  gcaCGTTGGCGCGAAGATGCTCACGCCAACGTGCTGTTACTTGCCGGAAAACGCCCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCAACTCCTCTACGCTACCGTCCGc                      >  3:58624/3‑147 (MQ=255)
                                                                                                                                      gTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCAACTCCTCTACGCTACCGTCCGCCAGCGa                <  5:297444/149‑1 (MQ=255)
                                                                                                                                      gTTGGCGTGAAGATGCTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCACCAACTCCTCTACGCTACCGTCCGCCAGCGa                >  7:216358/1‑149 (MQ=255)
                                                                                                                                         ggCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCAACTCCTCTACGCTACCGTCCGCCAGCGAct              <  8:370761/148‑1 (MQ=255)
                                                                                                                                            gTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCAACTCCTCTACGCTACCGTCCGCCAGCGActgctg          >  5:412876/1‑149 (MQ=255)
                                                                                                                                             tGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCAACTCCTCTACGCTACCGTCCGCCAGCGActgctg          >  7:157989/1‑148 (MQ=255)
                                                                                                                                                    cTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCAACTCCTCTACGCTACCGTCCGCCAGCGACTGCTGGATGatat  >  1:20893/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGGTGCCGGAGATTCGTTCATTGCCGGATTCCTTTGCGGCTGGTCTGCGGGGATGACATTACCGCAGGCGATAGCGCAGGGAACGGCGTGCGCGGCGAAAACCATTCAGTACCACGGTGCCTGGTAGGTATAACGTTGGCGTGAGCATCTTCACGCCAACGTGCTGTTACTTGCCGGAAAACGACCCTATAATCCGAGTAATTCATTCTTTATTTCAGGGTCGATTATGTCAGCTACGGACCGCTACTCTCATCAACTCCTCTACGCTACCGTCCGCCAGCGACTGCTGGATGATAT  >  NC_000913/3503826‑3504122

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