Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,245,113 C→A 6.5% intergenic (‑38/+116) 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,1130CA 6.5% 450.4 / 14.2 150intergenic (‑38/+116)malF/malEmaltose transporter subunit/maltose transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (82/58);  new base A (6/4);  total (88/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

Only 100 of 152 total aligned reads displayed.

TTGTGCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTGCGCGTCTTTCAGG  >  NC_000913/4244965‑4245262
                                                                                                                                                    |                                                                                                                                                     
ttGTGCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACg                                                                                                                                                      <  2:102587/150‑1 (MQ=255)
  gtgCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGcc                                                                                                                                                    <  1:301670/150‑1 (MQ=255)
  gtgCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGcc                                                                                                                                                    <  1:1324080/150‑1 (MQ=255)
   tgCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCt                                                                                                                                                   >  2:12915/1‑150 (MQ=255)
   tgCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCt                                                                                                                                                   >  2:572663/1‑150 (MQ=255)
     cGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTa                                                                                                                                                 >  2:431882/1‑150 (MQ=255)
       tACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATc                                                                                                                                               <  2:452555/150‑1 (MQ=255)
       tACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATc                                                                                                                                               >  2:1473351/1‑150 (MQ=255)
       tACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACAATTTCAGCGCGTCGCTATTCCACCAATGATTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATc                                                                                                                                               >  1:1012331/1‑150 (MQ=255)
         cATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCg                                                                                                                                             >  2:1660196/1‑150 (MQ=255)
         cATAAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACGAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCGTGAAAGCCTTATCCg                                                                                                                                             <  1:115936/150‑1 (MQ=255)
          aTTAAAACAACAAGGTAACCTACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCGTGAAAGCCTTATCCgt                                                                                                                                            <  2:1051559/150‑2 (MQ=255)
          aTTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGCCGGATAAGGCTTTCACGCCTTATCCgc                                                                                                                                            >  1:1365482/1‑149 (MQ=255)
          aTTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGa                                                                                                                                            >  1:536267/1‑150 (MQ=255)
          aTTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGATTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGa                                                                                                                                            >  1:1223047/1‑150 (MQ=255)
          aTTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGAGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCGTGAAAGCCTTATCCgt                                                                                                                                            <  1:343457/150‑2 (MQ=255)
           ttAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGac                                                                                                                                           >  2:288255/1‑150 (MQ=255)
           ttAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGac                                                                                                                                           >  1:1680931/1‑150 (MQ=255)
            tAAACCAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGaca                                                                                                                                          >  1:732487/1‑150 (MQ=255)
            tAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGaca                                                                                                                                          >  1:1456945/1‑150 (MQ=255)
            tAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGaca                                                                                                                                          >  2:609608/1‑150 (MQ=255)
             aaaaCAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGacaa                                                                                                                                         >  2:565914/1‑150 (MQ=255)
             aaaaCAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGacaa                                                                                                                                         <  1:471235/150‑1 (MQ=255)
               aacaacAATGTAACCCACCAGCATGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGacaaca                                                                                                                                       <  2:1246819/150‑1 (MQ=255)
               aacaacAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGacacca                                                                                                                                       >  2:96642/1‑150 (MQ=255)
                   acaaGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCATCATTACAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTg                                                                                                                                   >  2:1386957/1‑150 (MQ=255)
                      tgttAACCCACCAGCAGGCCTAGCATACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCAACCATTGTTTATTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCt                                                                                                                                <  1:1329386/147‑1 (MQ=255)
                      aGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCt                                                                                                                                >  1:1781575/1‑150 (MQ=255)
                              caccaGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGAc                                                                                                                        >  2:636585/1‑150 (MQ=255)
                              caccaGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGATTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGAc                                                                                                                        <  2:1223047/150‑1 (MQ=255)
                                  agcagGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTg                                                                                                                    >  2:960495/1‑150 (MQ=255)
                                    cagGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGAc                                                                                                                  >  1:506519/1‑150 (MQ=255)
                                    cagGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGAc                                                                                                                  >  1:616603/1‑150 (MQ=255)
                                     agGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCTATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGAcg                                                                                                                 >  2:342986/1‑150 (MQ=255)
                                      ggCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGAcgc                                                                                                                >  1:151036/1‑150 (MQ=255)
                                       gCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGAcgcg                                                                                                               <  1:1272923/150‑1 (MQ=255)
                                          gAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCt                                                                                                            <  2:1640778/150‑1 (MQ=255)
                                                aCCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAg                                                                                                      >  2:371133/1‑150 (MQ=255)
                                                 ccTATCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAgg                                                                                                     <  1:355963/150‑1 (MQ=255)
                                                 ccTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAgg                                                                                                     >  2:844705/1‑150 (MQ=255)
                                                 ccTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAgg                                                                                                     <  2:220526/150‑1 (MQ=255)
                                                  cTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGc                                                                                                    <  2:1290229/150‑1 (MQ=255)
                                                  cTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGc                                                                                                    <  2:353501/150‑1 (MQ=255)
                                                   tAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGcc                                                                                                   <  2:990265/150‑1 (MQ=255)
                                                        cTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAt                                                                                              >  1:760774/1‑150 (MQ=255)
                                                             cATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGtt                                                                                         <  2:151036/150‑1 (MQ=255)
                                                              aTTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACAAACGttt                                                                                        >  1:1694961/1‑150 (MQ=255)
                                                                ttCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCg                                                                                      <  2:1825243/150‑1 (MQ=255)
                                                                        gTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTa                                                                              >  1:501170/1‑150 (MQ=255)
                                                                           gCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGc                                                                           >  2:167165/1‑150 (MQ=255)
                                                                             tttGCCAGCTATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATATGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCg                                                                         <  1:1518055/150‑1 (MQ=255)
                                                                                  caccaATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACaa                                                                    >  1:980651/1‑150 (MQ=255)
                                                                                  caccaATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACaa                                                                    >  1:57517/1‑150 (MQ=255)
                                                                                       aTGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCAATCCAGGACGGATAATGCTTTCACGCCTTTTCCGACAACAACTGCCTGATGCGTCGCTGACGAGTCTTATCAGGCGTACATTCTTTGCGGTTTTGTAGGCCGGACAAGGCGt                                                               <  1:1620332/150‑1 (MQ=255)
                                                                                             cTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTc                                                             <  1:136694/146‑1 (MQ=255)
                                                                                             cTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTc                                                             >  2:136694/1‑146 (MQ=255)
                                                                                                tttAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCAcgccgc                                                      >  1:946923/1‑150 (MQ=255)
                                                                                                 ttAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCa                                                     >  1:833681/1‑150 (MQ=255)
                                                                                                   aaTGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATc                                                   <  1:1288818/150‑1 (MQ=255)
                                                                                                     tGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCg                                                 >  1:98697/1‑150 (MQ=255)
                                                                                                       tcATCCATGGGGTTCTTCGTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGc                                               <  2:239152/149‑1 (MQ=255)
                                                                                                       aCATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGc                                               >  2:1010787/1‑150 (MQ=255)
                                                                                                        cATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCa                                              <  1:1747522/150‑1 (MQ=255)
                                                                                                         aTCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGATGCGGCGTGAACGCCTTGTCCGGcct                                             <  2:1680931/150‑3 (MQ=37)
                                                                                                          tCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCAtt                                            <  1:371133/150‑1 (MQ=255)
                                                                                                          tCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCAtt                                            <  1:1706507/150‑1 (MQ=255)
                                                                                                            cATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTc                                          <  2:455335/150‑1 (MQ=255)
                                                                                                              gggggTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCGGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTcac                                        <  1:1660196/149‑1 (MQ=255)
                                                                                                                gggTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCTGACAAGGCGTTCACGCCGCATCCGGCATTTCACAg                                      >  1:1600758/1‑150 (MQ=255)
                                                                                                                gggTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAg                                      <  1:288255/150‑1 (MQ=255)
                                                                                                                gggTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAg                                      <  1:12915/150‑1 (MQ=255)
                                                                                                                  gTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCa                                    >  1:23368/1‑150 (MQ=255)
                                                                                                                    tcttcCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCAtt                                  <  2:1456945/150‑1 (MQ=255)
                                                                                                                        ccTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACtt                              >  2:700467/1‑150 (MQ=255)
                                                                                                                         cTCATTCCAGGAGGGATAAGGCTTTAACGCCTTATCCGACAACAACTGCCTGATGCGAAGCTGAAGTGTGTTATCAGGCCTACTTACGTTTGGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTg                             <  2:913551/150‑1 (MQ=38)
                                                                                                                          tCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTgg                            <  1:431882/150‑1 (MQ=255)
                                                                                                                           cATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGTCAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGt                           <  2:1012331/150‑1 (MQ=255)
                                                                                                                           cATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGt                           >  2:841983/1‑150 (MQ=255)
                                                                                                                           cATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGt                           <  2:1209885/150‑1 (MQ=255)
                                                                                                                           aaTTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGGCAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGt                           <  2:1112182/149‑1 (MQ=255)
                                                                                                                            aTTCCTGGACGGATAAGTCTTTCAAGCCTTATCCGACAACAACTGCCTTATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGGTTTGTAGGCCGGACAAGGCGTTCACGCCGAATCCGGCATTTCACAGCATTACTTGGTg                          <  2:732485/150‑1 (MQ=38)
                                                                                                                            aTTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATTCGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTg                          <  2:788361/150‑1 (MQ=255)
                                                                                                                            aTTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTg                          >  2:692482/1‑150 (MQ=255)
                                                                                                                            aTTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTg                          <  2:732487/150‑1 (MQ=255)
                                                                                                                                cAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATAc                      >  2:592315/1‑150 (MQ=255)
                                                                                                                                cAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATAc                      >  2:671291/1‑150 (MQ=255)
                                                                                                                                 aGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACg                     >  2:348939/1‑150 (MQ=255)
                                                                                                                                    aCGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGt                  >  2:1785375/1‑150 (MQ=255)
                                                                                                                                     cGGATAAGGCGTGAAAGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTc                 >  2:229722/1‑150 (MQ=255)
                                                                                                                                     cGGATAAGGCGTGAAAGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTc                 >  1:1198382/1‑150 (MQ=255)
                                                                                                                                     cGGATAAGGCGTGAAAGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTc                 >  2:939869/1‑150 (MQ=255)
                                                                                                                                      ggATAAGGCGTGAAAGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCt                >  1:733145/1‑150 (MQ=255)
                                                                                                                                      ggATAAGGCGTGAAAGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCt                >  1:1471664/1‑150 (MQ=255)
                                                                                                                                           aGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTGcgcg           >  2:1530389/1‑150 (MQ=255)
                                                                                                                                             gCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTGCGCGTc         >  2:1563171/1‑150 (MQ=255)
                                                                                                                                             gCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTGCGCGTc         <  1:1059424/150‑1 (MQ=255)
                                                                                                                                                   aCGCCTTATCCGACAACAACTGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTGCGCGTCTATCAg   >  1:364148/1‑150 (MQ=255)
                                                                                                                                                   aCGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTGCGCGTCTTTCAg   >  2:1502160/1‑150 (MQ=255)
                                                                                                                                                    cGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTGCGCGTCTTTCAgg  >  2:782353/1‑150 (MQ=255)
                                                                                                                                                    cGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTGCGCGTCTTTCAgg  >  1:523978/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTGTGCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTGCGCGTCTTTCAGG  >  NC_000913/4244965‑4245262

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