Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,245,111 C→A 6.6% intergenic (‑36/+118) malF ← / ← malE maltose ABC transporter membrane subunit MalF/maltose ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,245,1110CA 6.6% 442.6 / 12.8 131intergenic (‑36/+118)malF/malEmaltose ABC transporter membrane subunit MalF/maltose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (61/61);  new base A (6/3);  total (67/64)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 131 total aligned reads displayed.

CCTTGTGCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTGCGCGTCTTTC  >  NC_000913/4244963‑4245259
                                                                                                                                                    |                                                                                                                                                    
ccTTGTGCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCa                                                                                                                                                     >  2:212996/1‑150 (MQ=255)
     tgCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCt                                                                                                                                                <  2:553920/150‑1 (MQ=255)
      gCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCtt                                                                                                                                               >  1:771331/1‑150 (MQ=255)
      gCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCtt                                                                                                                                               <  2:625690/150‑1 (MQ=255)
       cGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTa                                                                                                                                              >  2:134472/1‑150 (MQ=255)
          aCATTAAAACAATAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATcc                                                                                                                                           >  1:1160187/1‑150 (MQ=255)
             ttAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGac                                                                                                                                        >  2:979046/1‑150 (MQ=255)
               aaaaCAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCGTGAAAGCCTTATCCgtcct                                                                                                                                      <  1:1320365/150‑5 (MQ=255)
               aaaaCAACAAGCTAACCCACCAGCAGGCCGATCAGACCTACCACTGCCCATTTCAGCGCGTCGCTTTGCCACCAAAGTATCTATTTAATCACATACATCGGGTTCTTCCTCAATCCAGCACGGATAAGTCTTTCACGCCTTATCCGacaa                                                                                                                                      >  2:9371/1‑150 (MQ=255)
                 aacaacAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGacaaca                                                                                                                                    <  1:212996/150‑1 (MQ=255)
                   caacaaGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACaacaac                                                                                                                                  >  1:1097264/1‑150 (MQ=255)
                     acaaGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGAGCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGAGATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTg                                                                                                                                <  2:527160/150‑1 (MQ=255)
                           tAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGAt                                                                                                                          >  2:524203/1‑150 (MQ=255)
                           tAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGAt                                                                                                                          >  2:868556/1‑150 (MQ=255)
                               ccaccaGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGa                                                                                                                      >  2:690863/1‑150 (MQ=255)
                                caccaGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGAc                                                                                                                     <  2:503572/150‑1 (MQ=255)
                                caccaGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGAc                                                                                                                     >  2:344586/1‑150 (MQ=255)
                                caccaGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGAc                                                                                                                     <  1:703910/150‑1 (MQ=255)
                                  ccaGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGGCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGc                                                                                                                   <  2:697569/150‑1 (MQ=255)
                                    agcagGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTg                                                                                                                 >  2:247966/1‑150 (MQ=255)
                                       agGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGAcg                                                                                                              <  1:959741/150‑1 (MQ=255)
                                         gCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGAcgcg                                                                                                            >  1:214659/1‑150 (MQ=255)
                                          ccGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGt                                                                                                           >  1:967180/1‑150 (MQ=255)
                                            gAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGGTAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCt                                                                                                         <  2:863883/150‑1 (MQ=255)
                                            gAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCt                                                                                                         <  1:1134340/150‑1 (MQ=255)
                                            gAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCt                                                                                                         <  2:1401092/150‑1 (MQ=255)
                                                 gACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGAcgcg                                                                                                            <  1:576130/142‑1 (MQ=255)
                                                 gACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGAcgcg                                                                                                            >  2:576130/1‑142 (MQ=255)
                                                   ccTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAgg                                                                                                  >  1:643817/1‑150 (MQ=255)
                                                    cTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGc                                                                                                 >  1:291211/1‑150 (MQ=255)
                                                     tAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGcc                                                                                                <  2:1367988/150‑1 (MQ=255)
                                                       gCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTa                                                                                              <  1:269377/150‑1 (MQ=255)
                                                        cACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTAc                                                                                             <  2:1031369/150‑1 (MQ=255)
                                                        cACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTAc                                                                                             <  1:524203/150‑1 (MQ=255)
                                                        cACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTAc                                                                                             <  1:868556/150‑1 (MQ=255)
                                                         aCTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACa                                                                                            <  1:1422745/150‑1 (MQ=255)
                                                            gACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATaa                                                                                         >  1:2129/1‑149 (MQ=255)
                                                             aCCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGTTAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACg                                                                                        <  2:246104/150‑1 (MQ=255)
                                                              ccATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATAAGGCCTACATACGt                                                                                       >  1:360029/1‑150 (MQ=255)
                                                                 tttCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTc                                                                                    <  2:771331/150‑1 (MQ=255)
                                                                 tttCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTc                                                                                    <  1:1471715/150‑1 (MQ=255)
                                                                  ttCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCg                                                                                   >  1:1259531/1‑150 (MQ=255)
                                                                   tCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCgg                                                                                  <  2:1368161/150‑1 (MQ=255)
                                                                   tCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCgg                                                                                  <  1:134472/150‑1 (MQ=255)
                                                                    cAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGt                                                                                 >  2:1149521/1‑150 (MQ=255)
                                                                       cgcgTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGtttt                                                                              <  1:1220391/150‑1 (MQ=255)
                                                                        gcgTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTg                                                                             >  2:650101/1‑150 (MQ=255)
                                                                        gcgTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTg                                                                             <  1:1117018/150‑1 (MQ=255)
                                                                        gcgTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTg                                                                             <  1:650101/150‑1 (MQ=255)
                                                                         cgTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGt                                                                            <  1:69981/150‑1 (MQ=255)
                                                                           tCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAg                                                                          <  2:1254157/150‑1 (MQ=255)
                                                                             gCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGc                                                                        >  1:1434607/1‑150 (MQ=255)
                                                                               tttGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGcct                                                                      <  2:1118345/150‑2 (MQ=255)
                                                                               tttGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCg                                                                      <  1:690863/150‑1 (MQ=255)
                                                                                 tGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGa                                                                    >  2:676265/1‑150 (MQ=255)
                                                                                    caccaATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACaa                                                                 >  2:1349640/1‑150 (MQ=255)
                                                                                    caccaATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACaa                                                                 >  2:864537/1‑150 (MQ=255)
                                                                                     accaATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAg                                                                >  1:1365224/1‑150 (MQ=255)
                                                                                      ccaATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGAc                                                                   >  1:261269/1‑146 (MQ=255)
                                                                                      ccaATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGAc                                                                   <  2:261269/146‑1 (MQ=255)
                                                                                      ccaATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAgg                                                               <  1:562803/150‑1 (MQ=255)
                                                                                            tttCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCa                                                         >  2:1452141/1‑150 (MQ=255)
                                                                                             ttCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCAc                                                        >  2:712550/1‑150 (MQ=255)
                                                                                              tCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCAcg                                                       >  1:1426889/1‑150 (MQ=255)
                                                                                                 ttttAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCAcgccg                                                    >  2:1338816/1‑150 (MQ=255)
                                                                                                    tAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCAt                                                 >  2:1004900/1‑150 (MQ=255)
                                                                                                      aTGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATcc                                               <  1:676265/150‑1 (MQ=255)
                                                                                                       tGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATcc                                                                                                                                           <  1:729290/57‑1 (MQ=255)
                                                                                                       tGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATcc                                                                                                                                           >  2:729290/1‑57 (MQ=255)
                                                                                                       tGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGATGCGGCGTGAAcgccttgtccg                                              <  1:1076100/150‑8 (MQ=37)
                                                                                                          cATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCa                                           <  1:340209/150‑1 (MQ=255)
                                                                                                           aTCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGATGCGGCGTGAACGCCTTGTCCGGcct                                          <  2:214659/150‑3 (MQ=37)
                                                                                                           aTCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGATGCGGCGTGAACGCCTTGTCCGGcct                                          <  2:967180/150‑3 (MQ=37)
                                                                                                           aTCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGATGCGGCGTGAACGCCTTGTCCGGcct                                          <  1:1077917/150‑3 (MQ=37)
                                                                                                             ccATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCAttt                                        >  1:1247820/1‑150 (MQ=255)
                                                                                                               aTGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCATGATGCGACGCTGCCGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGCACAAGGCGTTCACGCCGCATCCGGCATTTca                                      >  2:781983/1‑150 (MQ=255)
                                                                                                                tGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTcac                                     >  2:599448/1‑150 (MQ=255)
                                                                                                                tGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTcac                                     <  1:1004900/150‑1 (MQ=255)
                                                                                                                  gggTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAg                                   >  1:983664/1‑150 (MQ=255)
                                                                                                                   ggTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGc                                  >  1:1532223/1‑150 (MQ=255)
                                                                                                                    gTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCa                                 <  1:756063/150‑1 (MQ=255)
                                                                                                                         tcCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACt                            >  2:1461535/1‑150 (MQ=255)
                                                                                                                          ccTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACtt                           >  1:1130339/1‑150 (MQ=255)
                                                                                                                           cTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTg                          >  1:221628/1‑150 (MQ=255)
                                                                                                                              aTTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTg                       >  2:383311/1‑150 (MQ=255)
                                                                                                                               ttCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGa                      <  1:797442/150‑1 (MQ=255)
                                                                                                                                tCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGAt                     <  2:221628/150‑1 (MQ=255)
                                                                                                                                tCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGACCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGAt                     <  1:1149521/150‑1 (MQ=255)
                                                                                                                                 ccAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATa                    <  1:1292855/150‑1 (MQ=255)
                                                                                                                                  cAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATAc                   <  1:344586/150‑1 (MQ=255)
                                                                                                                                  cAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGACGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATAc                   >  1:902801/1‑150 (MQ=255)
                                                                                                                                   aGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGTCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACg                  >  1:1530192/1‑150 (MQ=255)
                                                                                                                                        ggATAAGGCGTGAAAGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGt                        >  2:1463807/1‑139 (MQ=255)
                                                                                                                                         gATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTg            >  2:752469/1‑150 (MQ=255)
                                                                                                                                         gATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTg            <  2:1160187/150‑1 (MQ=255)
                                                                                                                                         gATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTg            >  2:30870/1‑150 (MQ=255)
                                                                                                                                         gATAAGGCGTGAAAGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACtt                           <  2:72409/135‑1 (MQ=255)
                                                                                                                                           tAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTAACGCCGCATCCGGCATTTAACAGCATTACTTTGTGATACGAGTCTgcg          >  1:748519/1‑150 (MQ=255)
                                                                                                                                             aGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTGcgcg        <  1:954871/150‑1 (MQ=255)
                                                                                                                                              ggCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTGCGCGt       >  1:1304957/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCTTGTGCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTGCGCGTCTTTC  >  NC_000913/4244963‑4245259

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