Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,218,597 T→C 6.0% M1T (ATG→ACG) † aceK → isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,218,5970TC 6.0% 186.5 / 2.4 66M1T (ATG→ACG) †aceKisocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base T (32/30);  new base C (2/2);  total (34/32)
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

GGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGT  >  NC_000913/4218449‑4218745
                                                                                                                                                    |                                                                                                                                                    
ggCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAAt                                                                                                                                                      >  2:563095/1‑149 (MQ=255)
 gCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATg                                                                                                                                                     >  2:282361/1‑149 (MQ=255)
   gTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGcc                                                                                                                                                   <  1:388330/149‑1 (MQ=255)
       aCGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGt                                                                                                                                               <  2:605141/149‑1 (MQ=255)
       aCGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAACGCCGCGt                                                                                                                                               >  2:331129/1‑149 (MQ=255)
       aCGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAACGCCGCGt                                                                                                                                               >  2:408841/1‑149 (MQ=255)
         gccgcATCCGGCAATCTGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCGGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTgg                                                                                                                                             <  1:148802/149‑1 (MQ=255)
            gcATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGCACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATTTCTCTGGTGCTGATGAGGGCGCTAAATGCCGCGTGGCCt                                                                                                                                          >  1:83365/1‑149 (MQ=255)
                 cGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAat                                                                                                                                     <  2:632757/149‑1 (MQ=255)
                   gCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAatta                                                                                                                                   >  2:643969/1‑149 (MQ=255)
                     aaTCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAattatt                                                                                                                                 <  1:928403/149‑1 (MQ=255)
                     aaTCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAattatt                                                                                                                                 <  2:618375/149‑1 (MQ=255)
                        cGGTGCACGATGCCTGATGCGACGCTTGCGGGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTattgat                                                                                                                              <  1:622886/149‑1 (MQ=255)
                            gCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCt                                                                                                                          <  1:437920/149‑1 (MQ=255)
                                 aTGCCTGATGCGACGCTTGCGCGTCTTATCATGACTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAAc                                                                                                                     <  2:742274/149‑1 (MQ=255)
                                       gATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATTAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTTAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCAtttt                                                                                                               >  1:985227/1‑149 (MQ=255)
                                            gACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAg                                                                                                          >  1:173344/1‑149 (MQ=255)
                                             aCGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGGGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAgg                                                                                                         <  2:925504/149‑1 (MQ=255)
                                               gCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCt                                                                                                       >  2:853244/1‑149 (MQ=255)
                                                cTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCtt                                                                                                      <  1:261829/149‑1 (MQ=255)
                                                  tGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCg                                                                                                    <  2:135847/149‑1 (MQ=255)
                                                        tattATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGTGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTc                                                                                              <  1:578866/147‑1 (MQ=255)
                                                         cTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATACGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAACCCATTTTGCAAGGCTTCGATgagca                                                                                             >  1:532234/1‑145 (MQ=255)
                                                                 gCCTAAAGCCTTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATGCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATgggc                                                                                     <  2:83365/149‑3 (MQ=255)
                                                                       cgCCGTTGCCGAACGTAGGCTGGATCAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCTTTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTcc                                                                               <  2:568859/148‑1 (MQ=255)
                                                                             tGCCGATCGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAg                                                                         <  1:695845/149‑1 (MQ=255)
                                                                             tGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGatg                                                                         >  2:141730/1‑147 (MQ=255)
                                                                              gCCGAACGTAGGCTTGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGt                                                                        <  2:372318/149‑1 (MQ=255)
                                                                               ccGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTg                                                                       <  1:15810/149‑1 (MQ=255)
                                                                                cGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGa                                                                      >  1:444841/1‑149 (MQ=255)
                                                                                    cGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCAGATGAGGGAGCTAAATGCCGCGGGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGTCTTCTATGCTCAGTATGGTCGATTCCTCTAAGTGACCTc                                                                  >  1:996023/1‑149 (MQ=39)
                                                                                          cTGGATTAGGCGTTTACGCCGCATCCGGCAATTCTATGATCGTTATGAGGGCGCTAAATGCAGCGTGGCCTGGAATTTTTGATTGCGCAAACCATTTTGAAAGTCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTACGGTgc                                                            >  2:934583/1‑149 (MQ=255)
                                                                                          cTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGATCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTgc                                                            >  2:247183/1‑149 (MQ=255)
                                                                                           tGGATAAGGCGTTTACGCCGCATCCGGCAAATCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATAGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTgcg                                                           >  2:780850/1‑149 (MQ=255)
                                                                                              aTAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCgcag                                                        >  1:209673/1‑149 (MQ=255)
                                                                                                  ggCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAAACATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGCCCTCCGGTGCGCagcagc                                                    >  2:937255/1‑149 (MQ=255)
                                                                                                   gCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGCACTCCGGTGCGCAGCAGCg                                                   >  1:51171/1‑149 (MQ=255)
                                                                                                   gCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTACGGTGCGCCGCAGCg                                                   >  2:700859/1‑149 (MQ=255)
                                                                                                    cGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGt                                                  <  1:282361/149‑1 (MQ=255)
                                                                                                     gTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGtt                                                 >  1:307699/1‑149 (MQ=255)
                                                                                                     gTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGGCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGtt                                                 <  2:546137/149‑1 (MQ=255)
                                                                                                     gTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGGGGCATGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTTTGGTCGATTCCTCGAAGTGACCGACGGTGCGCAGCAGCGtt                                                 >  1:548204/1‑149 (MQ=255)
                                                                                                          cgccgcATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGaa                                            <  1:693398/149‑1 (MQ=255)
                                                                                                           gccgcATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGACGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAAc                                           >  1:338309/1‑149 (MQ=255)
                                                                                                            ccgcATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACa                                          <  1:853244/149‑1 (MQ=255)
                                                                                                              gcATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAgg                                        <  1:409893/149‑1 (MQ=255)
                                                                                                              gcATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAgg                                        <  1:563095/149‑1 (MQ=255)
                                                                                                                    ggCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCCTATTGCAAGGCTTCGAAGCTCAGTATTGTCGATTCCTCGAACTGACCTCCGGTGCGCCGCAGCGTTACGAACAGCCCGACt                                  >  2:1029380/1‑149 (MQ=255)
                                                                                                                    ggCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACt                                  <  1:643969/149‑1 (MQ=255)
                                                                                                                    ggCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACt                                  >  2:376904/1‑149 (MQ=255)
                                                                                                                     gCAATTCTCTGCTCCTGATGAGGGCGCTAAACGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTg                                 <  1:331129/149‑1 (MQ=255)
                                                                                                                     gCAATTCTCTGCTCCTGATGAGGGCGCTAAACGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTg                                 <  1:408841/149‑1 (MQ=255)
                                                                                                                       aaTTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGc                               <  1:797821/149‑1 (MQ=255)
                                                                                                                          tctcTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAAACATTTTGCAAGGCTTCGATGCTCAGTATGGTAGATTCCTCGAAGTGACCTGCGGTGCGCATCAGCGTTTCGAACAGGCCGACTGGCATg                            <  1:362478/149‑1 (MQ=255)
                                                                                                                               gCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAATCCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGCCCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTc                       >  2:355552/1‑149 (MQ=255)
                                                                                                                               gCTCCTGATGACGGCGCTACATGCCGCGTGGCCTGGAAGAATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCCGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTc                       <  2:854811/149‑1 (MQ=255)
                                                                                                                                   cTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCcagc                   >  1:556328/1‑149 (MQ=255)
                                                                                                                                   cTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCcagc                   >  2:934696/1‑149 (MQ=255)
                                                                                                                                   cTGATGAGGGCGCTAAATGCCGCGTGGCCGGGATTTATTGATTGCTCAACCCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCcagc                   >  2:455138/1‑149 (MQ=255)
                                                                                                                                        gaGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCg              >  2:300211/1‑149 (MQ=255)
                                                                                                                                        gaGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCg              >  2:296125/1‑149 (MQ=255)
                                                                                                                                         aGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGa             <  2:307699/149‑1 (MQ=255)
                                                                                                                                          gggCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCCCAGCAGCGTTTCGAACATGCCGACTGCAATGCTGTCCAGCAGGCgct            >  2:906713/1‑147 (MQ=255)
                                                                                                                                              gcTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGaaa        >  1:646599/1‑149 (MQ=255)
                                                                                                                                               cTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCTTTTTGCAAGGCTTCGATGCGCCGTATGGTCGATTCCTCGAAGTGACCTCCGTTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGaaaa       <  1:546105/149‑1 (MQ=255)
                                                                                                                                                tAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGCCTGGCATGCTGTCCAGCAGGCGATGaaaac      >  1:1032827/1‑148 (MQ=255)
                                                                                                                                                 aaaTGCCGCATGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGaaaacc     >  1:883045/1‑147 (MQ=255)
                                                                                                                                                    tGCCGCGTTGCCTGTAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGt  <  2:860405/149‑1 (MQ=255)
                                                                                                                                                    tGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCCGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGt  >  1:494661/1‑149 (MQ=255)
                                                                                                                                                    tGCCGCGTGGCCGTGAATAATTGATTGCTCAAACCATTTTGCAAGGCTTCGCTGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGt  <  2:731454/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGGGCGCTAAATGCCGCGTGGCCTGGAATTATTGATTGCTCAAACCATTTTGCAAGGCTTCGATGCTCAGTATGGTCGATTCCTCGAAGTGACCTCCGGTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGT  >  NC_000913/4218449‑4218745

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