Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,245,108 T→G 6.2% intergenic (‑33/+121) 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,1080TG 6.2% 409.6 / 15.7 130intergenic (‑33/+121)malF/malEmaltose ABC transporter membrane subunit MalF/maltose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (61/61);  new base G (5/3);  total (66/64)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01

Only 100 of 132 total aligned reads displayed.

TTCCCCTTGTGCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTGCGCGT  >  NC_000913/4244959‑4245254
                                                                                                                                                     |                                                                                                                                                  
 tCCCCTTGTGCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCtt                                                                                                                                                   <  1:1045793/150‑1 (MQ=255)
    ccTTGTGCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCa                                                                                                                                                >  2:212996/1‑150 (MQ=255)
     cTTGTGCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCAc                                                                                                                                               >  2:703910/1‑150 (MQ=255)
         tgCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCt                                                                                                                                           <  2:553920/150‑1 (MQ=255)
          gCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCtt                                                                                                                                          <  2:625690/150‑1 (MQ=255)
          gCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCtt                                                                                                                                          >  1:771331/1‑150 (MQ=255)
           cGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTa                                                                                                                                         >  2:134472/1‑150 (MQ=255)
              aCATTAAAACAATAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATcc                                                                                                                                      >  1:1160187/1‑150 (MQ=255)
               cATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCg                                                                                                                                     >  2:562803/1‑150 (MQ=255)
               cATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCg                                                                                                                                     >  2:1292855/1‑150 (MQ=255)
               cATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCg                                                                                                                                     >  2:1078621/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)
                               tAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGAt                                                                                                                     >  2:868556/1‑150 (MQ=255)
                               tAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGAt                                                                                                                     >  2:524203/1‑150 (MQ=255)
                                   ccaccaGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGa                                                                                                                 >  2:690863/1‑150 (MQ=255)
                                    caccaGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGAc                                                                                                                <  1:703910/150‑1 (MQ=255)
                                    caccaGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGAc                                                                                                                <  2:503572/150‑1 (MQ=255)
                                      ccaGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGc                                                                                                              <  2:239347/150‑1 (MQ=255)
                                      ccaGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGc                                                                                                              >  1:935806/1‑150 (MQ=255)
                                      ccaGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGGCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGc                                                                                                              <  2:697569/150‑1 (MQ=255)
                                           agGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGAcg                                                                                                         <  1:959741/150‑1 (MQ=255)
                                            ggCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGAcgc                                                                                                        >  1:1295671/1‑150 (MQ=255)
                                             gCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGAcgcg                                                                                                       >  1:214659/1‑150 (MQ=255)
                                              ccGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGt                                                                                                      >  1:967180/1‑150 (MQ=255)
                                                gAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGGTAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCt                                                                                                    <  2:863883/150‑1 (MQ=255)
                                                gAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCt                                                                                                    <  2:1401092/150‑1 (MQ=255)
                                                  gCAGACCTAGCACTGACCATTTCAGCGCGGCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTa                                                                                                  <  1:422174/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)
                                                       ccTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTAACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAgg                                                                                             >  1:1138671/1‑150 (MQ=255)
                                                        cTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGc                                                                                            >  1:291211/1‑150 (MQ=255)
                                                         tAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGcc                                                                                           <  1:1112937/150‑1 (MQ=255)
                                                           gCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTa                                                                                         >  1:1368161/1‑150 (MQ=255)
                                                            cACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTAc                                                                                        <  1:524203/150‑1 (MQ=255)
                                                            cACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTAc                                                                                        <  1:868556/150‑1 (MQ=255)
                                                            cACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTAc                                                                                        <  2:1031369/150‑1 (MQ=255)
                                                             aCTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACa                                                                                       <  1:1422745/150‑1 (MQ=255)
                                                             aCTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACa                                                                                       <  2:733774/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                                                                               <  1:1471715/150‑1 (MQ=255)
                                                                     tttCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTc                                                                               >  2:317575/1‑150 (MQ=255)
                                                                     tttCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTc                                                                               <  2:771331/150‑1 (MQ=255)
                                                                      ttCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCg                                                                              >  1:1259531/1‑150 (MQ=255)
                                                                       tCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCgg                                                                             <  2:1368161/150‑1 (MQ=255)
                                                                        cAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGt                                                                            >  2:796824/1‑150 (MQ=255)
                                                                            gcgTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTg                                                                        >  2:650101/1‑150 (MQ=255)
                                                                            gcgTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTg                                                                        <  1:1117018/150‑1 (MQ=255)
                                                                             cgTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGt                                                                       <  2:643817/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)
                                                                                        caccaATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACaa                                                            >  2:864537/1‑150 (MQ=255)
                                                                                        caccaATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACaa                                                            >  2:1349640/1‑150 (MQ=255)
                                                                                         accaATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAg                                                           >  1:1365224/1‑150 (MQ=255)
                                                                                          ccaATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGAc                                                              <  2:261269/146‑1 (MQ=255)
                                                                                          ccaATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGAc                                                              >  1:261269/1‑146 (MQ=255)
                                                                                                tttCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCa                                                    <  1:796824/150‑1 (MQ=255)
                                                                                                tttCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCa                                                    >  2:1452141/1‑150 (MQ=255)
                                                                                                     ttttAATGACATCTATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCAcgccg                                               <  1:680562/150‑1 (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)
                                                                                                             aCATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCAACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGATGCGGCGTGAACGCCTTGTCCGGc                                       <  1:220968/150‑1 (MQ=37)
                                                                                                              cATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCa                                      <  1:340209/150‑1 (MQ=255)
                                                                                                               aTCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGATGCGGCGTGAACGCCTTGTCCGGcct                                     <  1:1077917/150‑3 (MQ=37)
                                                                                                               aTCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGATGCGGCGTGAACGCCTTGTCCGGcct                                     <  2:967180/150‑3 (MQ=37)
                                                                                                                 ccATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCAttt                                   >  1:1247820/1‑150 (MQ=255)
                                                                                                                   aTGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCATGATGCGACGCTGCCGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGCACAAGGCGTTCACGCCGCATCCGGCATTTca                                 >  2:781983/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)
                                                                                                                         ttcttcCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCAt                           <  1:979046/150‑1 (MQ=255)
                                                                                                                              ccTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACtt                      >  1:1130339/1‑150 (MQ=255)
                                                                                                                               cTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTg                     >  1:221628/1‑150 (MQ=255)
                                                                                                                                 cATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGt                   <  1:317575/150‑1 (MQ=255)
                                                                                                                                  aTTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTg                  >  2:383311/1‑150 (MQ=255)
                                                                                                                                  aTTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTg                  <  2:179519/150‑1 (MQ=255)
                                                                                                                                   ttCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGa                 <  2:1365224/150‑1 (MQ=255)
                                                                                                                                    tCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGAt                <  2:221628/150‑1 (MQ=255)
                                                                                                                                    tCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGACCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGAt                <  1:1149521/150‑1 (MQ=255)
                                                                                                                                      cAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGACGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATAc              >  1:902801/1‑150 (MQ=255)
                                                                                                                                            ggATAAGGCGTGAAAGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGt                   >  2:1463807/1‑139 (MQ=255)
                                                                                                                                            ggATAAGGCGTGAAAGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGt                   <  1:1463807/139‑1 (MQ=255)
                                                                                                                                            ggATAAGGCGTGAAAGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCt        >  1:1535867/1‑150 (MQ=255)
                                                                                                                                            ggATAAGGCGTGAAAGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCt        >  1:1227687/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)
                                                                                                                                             gATAAGGCGTGAAAGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACtt                      >  1:72409/1‑135 (MQ=255)
                                                                                                                                               tAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTgcg     <  2:508834/150‑1 (MQ=255)
                                                                                                                                                  ggCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTGCGCGt  >  1:1304957/1‑150 (MQ=255)
                                                                                                                                                  ggCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTGCGCGt  >  1:607142/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TTCCCCTTGTGCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTGCGCGT  >  NC_000913/4244959‑4245254

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: