Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,245,108 T→G 5.2% intergenic (‑33/+121) 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,1080TG 5.2% 400.8 / 10.0 136intergenic (‑33/+121)malF/malEmaltose transporter subunit/maltose transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (60/68);  new base G (5/2);  total (66/70)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.76e-01

Only 100 of 140 total aligned reads displayed.

TTCCCCTTGTGCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTGCG  >  NC_000913/4244959‑4245251
                                                                                                                                                     |                                                                                                                                               
ttCCCCTTGTGCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCt                                                                                                                                                 >  1:886297/1‑150 (MQ=255)
ttCCCCTTGTGCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCt                                                                                                                                                 <  1:1424192/150‑1 (MQ=255)
 tCCCCTTGTGCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTATATAATGACATCCATGCGGTTCTTCCTCATTCCAGGACGGATAAGGCtt                                                                                                                                                >  1:265629/1‑150 (MQ=255)
  ccccTTGTGCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATACGGCttt                                                                                                                                               >  1:1159746/1‑150 (MQ=255)
      ttGTGCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACg                                                                                                                                           <  1:1177627/150‑1 (MQ=255)
      atgtgCGTACATTAAAACAACAAGGTAACCCACCATCATGCAGAGCAGACCTATCACTGAACATTTCAGCGCGTCCCTTTGCCACCAATGTTTCTTTTTAAAGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCagg                                                                                                                                           <  1:292852/149‑3 (MQ=255)
         tgCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCt                                                                                                                                        <  1:915548/150‑1 (MQ=255)
           cGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTa                                                                                                                                      >  1:125710/1‑150 (MQ=255)
            gTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTAt                                                                                                                                     >  1:626739/1‑150 (MQ=255)
             tACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCCGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATc                                                                                                                                    >  2:302926/1‑150 (MQ=255)
             tACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATc                                                                                                                                    >  2:1142103/1‑150 (MQ=255)
             tACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATc                                                                                                                                    >  1:409669/1‑150 (MQ=255)
             tACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATc                                                                                                                                    <  2:1460079/150‑1 (MQ=255)
              aCATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATcc                                                                                                                                   <  2:1688836/150‑1 (MQ=255)
                aTTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCGTGAAAGCCTTATCCgt                                                                                                                                 <  2:401409/150‑2 (MQ=255)
                 taaaaaCAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCATTGCCACAAAAGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGAAAAGGCTTTCACGACTTATCCGac                                                                                                                                >  1:1655153/3‑150 (MQ=255)
                  tAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGCCCTAGCACTGACCATTTCAGCGAGTCGCTATGCCACCAATGTTTCTTTTTAATGACATACATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACTCCTTATCCGaca                                                                                                                               >  1:336060/1‑150 (MQ=255)
                     aacaccATGGTAACCGACCAGCAGGCCGAGCAGACCTAGGACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGAGATCCATGGTGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGacaaca                                                                                                                            <  1:1051630/150‑1 (MQ=255)
                       caacaaGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACaacaac                                                                                                                          >  1:672349/1‑150 (MQ=255)
                       caacaaGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACaacaac                                                                                                                          >  2:231578/1‑150 (MQ=255)
                         acaccgaaaCCCACCAGCATGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGTATAAGGCTTTCACGCCTTATCCGACAACAACTg                                                                                                                        <  2:286771/143‑1 (MQ=255)
                         acaaGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTg                                                                                                                        >  2:330654/1‑150 (MQ=255)
                         acaaGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATCTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTg                                                                                                                        <  2:290802/150‑1 (MQ=255)
                             ggTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTg                                                                                                                    <  1:17317/150‑1 (MQ=255)
                             ggTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTg                                                                                                                    <  1:1366575/150‑1 (MQ=255)
                              gTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCACGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGa                                                                                                                   <  2:119719/150‑1 (MQ=255)
                               tAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTCTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGAt                                                                                                                  >  1:786626/1‑150 (MQ=255)
                                  cccACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGAtaag                                                                                                               >  2:1061189/1‑147 (MQ=255)
                                  cccACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCg                                                                                                               >  2:1271605/1‑150 (MQ=255)
                                    caccaGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGAc                                                                                                             >  2:516288/1‑150 (MQ=255)
                                    caccaGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGAc                                                                                                             <  2:1377485/150‑1 (MQ=255)
                                    caccaGACGGCCGATCAGACCTAGCACTGACCATTTCAGCGCTTCGCTTTGCCACAAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCATGACGGATAAGTCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGAc                                                                                                             <  2:1189672/150‑1 (MQ=255)
                                       cagcagGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCt                                                                                                          >  2:437270/1‑150 (MQ=255)
                                       cagcagGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTCTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCt                                                                                                          <  2:786626/150‑1 (MQ=255)
                                        agcagGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTg                                                                                                         <  2:1641227/150‑1 (MQ=255)
                                         gcagGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGa                                                                                                        <  1:1142103/150‑1 (MQ=255)
                                           agtcctaGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCGCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGAcg                                                                                                      <  1:23870/144‑1 (MQ=255)
                                           agGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGAcg                                                                                                      >  1:166512/1‑150 (MQ=255)
                                               cGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTc                                                                                                  >  1:32685/1‑150 (MQ=255)
                                                    aGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATc                                                                                             <  1:1475842/150‑1 (MQ=255)
                                                     gACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCGGACGCGTCTTATCa                                                                                            <  1:1001807/150‑1 (MQ=255)
                                                       ccTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAgg                                                                                          >  2:344477/1‑150 (MQ=255)
                                                       ccTAGCACTGACCATTTAAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTGCTTCCTCATTCCAGGACGGCTAAGGCTTACACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAgg                                                                                          <  2:248068/150‑1 (MQ=255)
                                                        cTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATAAGACGCGTCAGCGTCGCATCAGGc                                                                                         <  1:1061189/150‑1 (MQ=255)
                                                         tAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGcc                                                                                        <  1:367318/150‑1 (MQ=255)
                                                         tAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGAAGCGACGCTGACGCGTCTTATCAGGcc                                                                                        <  2:867109/150‑1 (MQ=255)
                                                          aGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCt                                                                                       <  2:1647317/150‑1 (MQ=255)
                                                          aGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATAAGACGCGTCAGCGTCGCATCAGGcag                                                                                       <  2:1714345/150‑3 (MQ=255)
                                                           gCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATAAGACGCGTCAGCGTCGCATCAGGcagt                                                                                      <  2:519707/150‑4 (MQ=255)
                                                            cACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTAc                                                                                     >  1:1142005/1‑150 (MQ=255)
                                                             aCTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACa                                                                                    <  2:78380/150‑1 (MQ=255)
                                                              cTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGCCATCCATGGGGTTCTTCCGCATTCCAGGACGGATAATGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAt                                                                                   <  2:975999/150‑1 (MQ=255)
                                                               tGACCATTTCATCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATCGGGTTCTTCCTCATTCCAGGACGGATAAGGCATTCACGCCTTATCCGACAACAACTCCCTGATGCGACGCTGACGCGACTTATCAGGCCTACaaa                                                                                  >  2:930891/1‑148 (MQ=38)
                                                                      ttCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCg                                                                           <  2:495290/150‑1 (MQ=255)
                                                                       tCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCgg                                                                          <  1:1124323/150‑1 (MQ=255)
                                                                        cAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGAcg                                                                                                      <  2:398971/121‑1 (MQ=255)
                                                                        cAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGAcg                                                                                                      >  1:398971/1‑121 (MQ=255)
                                                                           cgcgTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGtttt                                                                      >  1:1266310/1‑150 (MQ=255)
                                                                                cGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTAACCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAgg                                                                 >  2:185625/1‑150 (MQ=255)
                                                                                  cTTTGCCACCAATGTTTCTTTTTAATGACATACATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCAAGCCTTATCCGACAACAACTGCCTGATGCTACGCTGACGCTTCTTATCAGTCCTACATACGTTTCGGTTTTGTAggaa                                                               >  1:1589226/1‑148 (MQ=255)
                                                                                   tttGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCg                                                              <  2:1626144/150‑1 (MQ=255)
                                                                                   tttGCCACAAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTGTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCg                                                              <  2:121792/150‑1 (MQ=255)
                                                                                    ttGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCgg                                                             <  1:437270/150‑1 (MQ=255)
                                                                                       ccaccaATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCTGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACa                                                          >  2:960053/1‑150 (MQ=255)
                                                                                         accaATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAg                                                        >  1:624718/1‑150 (MQ=255)
                                                                                          ccaATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAgg                                                       <  1:302926/150‑1 (MQ=255)
                                                                                            aaTGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCg                                                     <  1:1075301/150‑1 (MQ=255)
                                                                                              tGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTATGCCGGACAAGGCGtt                                                   >  1:322220/1‑150 (MQ=255)
                                                                                              tGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGtt                                                   <  2:32685/150‑1 (MQ=255)
                                                                                               gTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTc                                                  >  2:92301/1‑150 (MQ=255)
                                                                                                    ttttttATGACTTCCATGGGGTTATTCCTCATTCCAGGACGGATAAGGCTTTCAGGCTTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGTACAAGGCGTTCAcgcc                                             <  2:1655153/150‑1 (MQ=38)
                                                                                                     ttttAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGCCAAGGCGTTCAcgccg                                            >  1:1483863/1‑150 (MQ=255)
                                                                                                       ttAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCa                                          >  1:608635/1‑150 (MQ=255)
                                                                                                            gACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCCCGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCgg                                     <  2:1031206/150‑1 (MQ=255)
                                                                                                                tCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCAtt                                 >  2:1314544/1‑150 (MQ=255)
                                                                                                                   aTGTGGTTCATCCTCATTCCAGGACGGATAAGGCTTTCCCGCCTTATCCGACAACAACTGCCAGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTCTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTca                              <  1:349509/150‑1 (MQ=38)
                                                                                                                    tGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTcac                             <  1:1271605/150‑1 (MQ=255)
                                                                                                                    tGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTcac                             >  2:1322422/1‑150 (MQ=255)
                                                                                                                    tGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTcac                             >  2:895794/1‑150 (MQ=255)
                                                                                                                        gTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCa                         >  2:124331/1‑150 (MQ=255)
                                                                                                                             tcCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACt                    >  1:150259/1‑150 (MQ=255)
                                                                                                                             tcCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACt                    >  1:300293/1‑150 (MQ=255)
                                                                                                                              ccTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACtt                   >  1:253478/1‑150 (MQ=255)
                                                                                                                              ccTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACtt                   >  2:723553/1‑150 (MQ=255)
                                                                                                                              ccTCATTCCAGGACGGATAAGGCTTTCACGCCTTATACGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCAGCATCCGGCATTTCACAGCATTACtt                   >  1:1180900/1‑150 (MQ=255)
                                                                                                                                 cATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGt                >  2:1175607/1‑150 (MQ=255)
                                                                                                                                 cATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGt                <  2:672349/150‑1 (MQ=255)
                                                                                                                                  aTTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTg               <  1:95224/150‑1 (MQ=255)
                                                                                                                                   ttCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGa              <  1:516288/150‑1 (MQ=255)
                                                                                                                                      cAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATAc           <  2:23964/150‑1 (MQ=255)
                                                                                                                                       aGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACg          >  2:433539/1‑150 (MQ=255)
                                                                                                                                         gtcGGATAAGGCGTGAAAGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACg          <  2:1709604/146‑1 (MQ=255)
                                                                                                                                         gtcGGATAAGGCGTGAAAGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACg          >  1:1709604/3‑148 (MQ=255)
                                                                                                                                         gtcGGATAAGGCGTGAAAGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAg        >  2:772391/3‑150 (MQ=255)
                                                                                                                                         gtcGGATAAGGCGTGAAAGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAg        >  1:654465/3‑150 (MQ=255)
                                                                                                                                         gACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAg        >  1:740239/1‑150 (MQ=255)
                                                                                                                                         gACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAg        <  1:1710481/150‑1 (MQ=255)
                                                                                                                                            ggATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCt     >  1:701708/1‑150 (MQ=255)
                                                                                                                                             gATAAGGCTTTAACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTAATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTg    <  2:873039/150‑1 (MQ=255)
                                                                                                                                               tAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTgcg  >  1:527853/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
TTCCCCTTGTGCGTACATTAAAACAACAAGGTAACCCACCAGCAGGCCGAGCAGACCTAGCACTGACCATTTCAGCGCGTCGCTTTGCCACCAATGTTTCTTTTTAATGACATCCATGGGGTTCTTCCTCATTCCAGGACGGATAAGGCTTTCACGCCTTATCCGACAACAACTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACATACGTTTCGGTTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCATTTCACAGCATTACTTGGTGATACGAGTCTGCG  >  NC_000913/4244959‑4245251

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