Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,245,110 T→G 5.2% intergenic (‑35/+119) malF ← / ← malE maltose transporter subunit/maltose transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,245,1100TG 5.2% 279.8 / 4.9 97intergenic (‑35/+119)malF/malEmaltose transporter subunit/maltose transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (42/50);  new base G (2/3);  total (44/53)
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.69e-01

CCCTTGTGCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTGCGCGTCTTTC  >  NC_000913/4244962‑4245259
                                                                                                                                                    |                                                                                                                                                     
cccTTGTGCGTACATTAAAACAACAAGGTACCCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTAAATGACATCCATGGTGTTCTTCCTCATTCCAGGACGGATAAGGCTTTc                                                                                                                                                      >  1:425074/1‑150 (MQ=255)
cccTTGTGCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTc                                                                                                                                                      >  2:406911/1‑150 (MQ=255)
 ccTTGTGCGTACAATAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCa                                                                                                                                                     <  2:1047526/150‑1 (MQ=255)
      tgCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCt                                                                                                                                                <  2:517312/150‑1 (MQ=255)
      tgCGTAAAATAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGGTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCt                                                                                                                                                <  2:73730/150‑1 (MQ=255)
       gCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCtt                                                                                                                                               >  2:693723/1‑150 (MQ=255)
        cGTACATTAAAACACCAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGCTAAGGCTTTCACGCCTTa                                                                                                                                              >  2:1110881/1‑150 (MQ=255)
         gTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTAt                                                                                                                                             >  2:728806/1‑150 (MQ=255)
          tACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATc                                                                                                                                            >  1:1235801/1‑150 (MQ=255)
          tACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATc                                                                                                                                            <  2:359800/150‑1 (MQ=255)
          tACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCGTGAAAGCCTTATc                                                                                                                                            <  1:496517/150‑1 (MQ=255)
           aCATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCACATTCCAGGACGGATAAGGCTTTTACGCCTTATcc                                                                                                                                           <  1:1210943/150‑1 (MQ=255)
            cATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTATCACTGACCATTTCCGCGCGTCGCTTTGCCACCAATGTGTCTTTTTAATGCCATCCATGGGGTTCTTCCGCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCg                                                                                                                                          >  2:978967/1‑150 (MQ=255)
            cATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCGTGAAAGCCTTATCCg                                                                                                                                          <  2:609976/150‑1 (MQ=255)
              ttAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGac                                                                                                                                        <  2:1058262/150‑1 (MQ=255)
              ttAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGac                                                                                                                                        >  1:30023/1‑150 (MQ=255)
              ttAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCGTGAAAGCCTTATCCgtc                                                                                                                                        <  1:634296/150‑3 (MQ=255)
                   acaacaAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGATTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTGCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGAcaacaa                                                                                                                                   <  2:710905/150‑1 (MQ=255)
                    caacaaGGTAACCCACCAGCAGGCCGAGAAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTGCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACaacaac                                                                                                                                  >  1:540839/1‑150 (MQ=255)
                      acaaGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTg                                                                                                                                >  2:663176/1‑150 (MQ=255)
                        aaGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCGTTCACGCCTTATCCGACAACAACTGcc                                                                                                                              <  2:850525/150‑1 (MQ=255)
                        aaGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATATCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTGCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGcc                                                                                                                              <  1:1099198/150‑1 (MQ=255)
                         aGGTACCCCACCAGCAGGCCGAGCAGACCTCGCACTGACCATTTCAGCGCGTAGCTTTGCCCGCAATCTTTCGATTTAATGACATCCCTGGGGTTCTTCCTCATTGCAGGACGGATACGGCTTTCCCGCCATAGCCGACAAAAACTGcca                                                                                                                             >  1:640662/1‑149 (MQ=255)
                          ggTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTg                                                                                                                            >  1:1145766/1‑150 (MQ=255)
                           gTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGa                                                                                                                           >  2:179864/1‑150 (MQ=255)
                              aCCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGc                                                                                                                        <  1:323729/150‑1 (MQ=255)
                                    cagcagGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCt                                                                                                                  <  2:517100/150‑1 (MQ=255)
                                       cagGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGAc                                                                                                               >  2:651326/1‑150 (MQ=255)
                                        ttgCCGAGCATACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAATGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGAcg                                                                                                              <  2:526513/148‑1 (MQ=255)
                                        agGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGAcg                                                                                                              <  1:1166100/150‑1 (MQ=255)
                                         ggCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGAcgc                                                                                                             >  2:77202/1‑150 (MQ=255)
                                          gCCGAGCAGACCTAGCACTGACCATTTCAGCGCGGCGCTTTGCCACCAAAGATTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCCGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGAcgcg                                                                                                            <  1:910458/150‑1 (MQ=255)
                                            cGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCCTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTc                                                                                                          <  1:940833/150‑1 (MQ=255)
                                            cGAGCAGACCTAGCACTGACCATTTCAGCGCGGCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATCAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTc                                                                                                          <  1:95195/150‑1 (MQ=255)
                                               gCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTa                                                                                                       <  1:693723/150‑1 (MQ=255)
                                                      tcgCACTCACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGCATAAGGCTTTCTCGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGcc                                                                                                <  2:329042/148‑1 (MQ=255)
                                                      tAGCACTGAGCATTTCAGCGCGTCGCTTTGCGACCAATGTTTGTTTTTAATGTCATCCATGGGGTTGTTCCTCTTTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGcc                                                                                                <  1:860378/150‑1 (MQ=38)
                                                      tAGCACTGACGCTTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTATAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGcc                                                                                                <  2:89845/150‑1 (MQ=255)
                                                      tAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGcc                                                                                                <  2:1039065/150‑1 (MQ=255)
                                                      tAGCACTGACCATTTCAGCGCGGCGCTTTGCCACCAATGTTTCTTTTTCATGACTTCCCTGGGGTGCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGcc                                                                                                <  1:33388/150‑1 (MQ=38)
                                                         cACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTAc                                                                                             >  1:502394/1‑150 (MQ=255)
                                                          aCTGACCATTTCAGCGCGTAGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGTGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACa                                                                                            <  1:903184/150‑1 (MQ=255)
                                                             gACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACATAc                                                                                         >  2:1032615/1‑150 (MQ=255)
                                                             gACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATAc                                                                                         <  2:771569/150‑1 (MQ=255)
                                                                cATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGtt                                                                                      >  2:717597/1‑150 (MQ=255)
                                                                   ttCAGCGCGTCGCTTTGCCACTAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATATGGATTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCg                                                                                   <  2:494671/150‑1 (MQ=255)
                                                                    tCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCgg                                                                                  >  1:1262186/1‑150 (MQ=255)
                                                                    tCAGCGCGGCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCgg                                                                                  <  2:1031812/150‑1 (MQ=255)
                                                                     aagagcgTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGt                                                                                 <  2:811570/146‑1 (MQ=255)
                                                                        cgcgTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGtttt                                                                              <  2:174586/150‑1 (MQ=255)
                                                                          cgTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAAGTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGt                                                                            <  1:7317/150‑1 (MQ=255)
                                                                          cgTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGt                                                                            >  2:547620/1‑150 (MQ=255)
                                                                          cgTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTGCTTCCTCATTCCTGTACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGt                                                                            <  1:231129/150‑1 (MQ=255)
                                                                           gTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTa                                                                           <  2:156216/150‑1 (MQ=255)
                                                                           gTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCATTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTa                                                                           <  1:820623/150‑1 (MQ=255)
                                                                            tCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACATCTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAg                                                                          <  1:28885/150‑1 (MQ=255)
                                                                            tCGCTTTGCCACAATTGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACTCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAg                                                                          <  1:1299301/150‑1 (MQ=255)
                                                                              gCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGc                                                                        <  1:556566/150‑1 (MQ=255)
                                                                                    ccaccaATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGCCGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACa                                                                  >  2:312654/1‑150 (MQ=255)
                                                                                     caccaATGTTTCTTTTTAATGACATCCATTGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTAACGCCTTATCCGACAACCACTGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAACATACGTTACGGATTTGTAGGCCGGACaa                                                                 >  1:1033262/1‑150 (MQ=255)
                                                                                        caATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGTGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGc                                                              <  1:837834/150‑1 (MQ=255)
                                                                                        caATGTTTCTTTTTAATGACATCCATGGGGGTCGTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAGCAACTGCCTGATGTGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGc                                                              <  1:567748/150‑1 (MQ=255)
                                                                                            gTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTc                                                          <  1:17318/150‑1 (MQ=255)
                                                                                            gTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTc                                                          >  2:733428/1‑150 (MQ=255)
                                                                                                gtttttAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCAcgc                                                      <  1:1298032/149‑1 (MQ=255)
                                                                                                  ttttAATGACATCCATTGGGAGCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTAATCCGACAACAACTGCCTGATGCGACGCTGACGCGGCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCAcgccg                                                    <  2:1128712/150‑1 (MQ=38)
                                                                                                  ttttAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCAcgccg                                                    <  2:1235801/150‑1 (MQ=255)
                                                                                                   tttAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCAcgccgc                                                   >  1:450866/1‑150 (MQ=255)
                                                                                                          ccATCCATGGGGTGCTTCCTCATTCCAGGACGTATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGc                                            <  1:690164/149‑1 (MQ=255)
                                                                                                           cATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACCAGGCGTACACGCCGCATCCGGCa                                           >  2:1179329/1‑150 (MQ=255)
                                                                                                           cATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCa                                           >  2:1184758/1‑150 (MQ=255)
                                                                                                           cATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCa                                           >  1:7187/1‑150 (MQ=255)
                                                                                                             tCCATTGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGATGCGGCGTGAACGCCTTGTCCGGccta                                         <  1:658139/150‑4 (MQ=25)
                                                                                                               cATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGATGCGGCGTGAACGCCTTGTCCGGcctaca                                       <  1:165688/150‑6 (MQ=37)
                                                                                                               cATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCAttta                                       >  2:1182749/1‑149 (MQ=255)
                                                                                                                  ggggTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTcaca                                    <  2:1269293/150‑1 (MQ=255)
                                                                                                                    ggTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCCCGCCGCATCCGGCATTTCACAGc                                  >  2:155469/1‑150 (MQ=255)
                                                                                                                    ggTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGc                                  >  1:771278/1‑150 (MQ=255)
                                                                                                                     gTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCa                                 <  2:517956/150‑1 (MQ=255)
                                                                                                                      ttcttcCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGAATTTCACAGCAt                                >  1:1233563/1‑150 (MQ=255)
                                                                                                                       tcttcCTCATTCCAGGACGGATAAGGCTTTCATGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGCCAAGGCGTTCACGCCGCATCCGGCATTTCACAGCAtt                               >  2:705714/1‑150 (MQ=255)
                                                                                                                         ttcCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCAGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTGACAGCATTAc                             <  2:1055760/150‑1 (MQ=255)
                                                                                                                          tcCTCATTCCAGGATGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACt                            <  1:367629/150‑1 (MQ=255)
                                                                                                                              cATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCCCGCCGCATCCGGCATTTCACAGCATTACTTGGt                        >  2:552503/1‑150 (MQ=255)
                                                                                                                               aTTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTg                       <  1:520688/150‑1 (MQ=255)
                                                                                                                                   cAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATAc                   >  2:46691/1‑150 (MQ=255)
                                                                                                                                     ggACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCATACATACGGTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCTGGCATTTCACAGCATTACTTGGTGATACGa                 >  2:211676/1‑150 (MQ=255)
                                                                                                                                      gACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAg                >  2:295168/1‑150 (MQ=255)
                                                                                                                                         ggATAAGGCGTGAAAGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGACATTTCACAGCATTACTTGGTGATACGAGTCt             >  2:1115409/1‑150 (MQ=38)
                                                                                                                                         ggATAAGGCGTGAAAGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGACATTTCACAGCATTACTTGGTGATACGAGTCt             >  2:1125430/1‑150 (MQ=38)
                                                                                                                                          gATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGCTGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTg            <  2:310893/150‑1 (MQ=255)
                                                                                                                                               ggCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCCCAGCATTACTTGGTGATACGAGTCTGCGCGt       >  2:745701/1‑150 (MQ=255)
                                                                                                                                               ggCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTGCGCGt       >  2:785484/1‑150 (MQ=255)
                                                                                                                                               ggCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTGCGCGt       >  2:881566/1‑150 (MQ=255)
                                                                                                                                               ggCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTGCGCGt       >  1:1161590/1‑150 (MQ=255)
                                                                                                                                                gCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTGCGCGTc      >  2:588863/1‑150 (MQ=255)
                                                                                                                                                gCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCAATACTTGGTGATACGAGTCTGCGCGTc      >  1:1275721/1‑150 (MQ=255)
                                                                                                                                                    tCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTGCGCGTCTTTc  <  1:1270264/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCCTTGTGCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTGCGCGTCTTTC  >  NC_000913/4244962‑4245259

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: