Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,537,640 T→G 5.5% intergenic (‑46/+19) nanS ← / ← nanM 9‑O‑acetyl N‑acetylneuraminic acid esterase/N‑acetylneuraminic acid mutarotase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,537,6400TG 5.5% 460.9 / 9.9 151intergenic (‑46/+19)nanS/nanM9‑O‑acetyl N‑acetylneuraminic acid esterase/N‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base T (84/57);  new base G (3/6);  total (87/64)
Fisher's exact test for biased strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 157 total aligned reads displayed.

TCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCC  >  NC_000913/4537491‑4537789
                                                                                                                                                     |                                                                                                                                                     
tCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCtt                                                                                                                                                       <  2:11511/150‑1 (MQ=255)
 cTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCttt                                                                                                                                                      >  1:557705/1‑150 (MQ=255)
 cTTCCCTGTCCGGTAAAGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTCATAGCAGGGCCCCttt                                                                                                                                                      >  2:1698027/1‑150 (MQ=255)
        gTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCgggg                                                                                                                                               >  2:1024556/1‑150 (MQ=255)
          ccGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGcc                                                                                                                                             <  2:710422/150‑1 (MQ=255)
           cGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGccc                                                                                                                                            <  2:414363/150‑1 (MQ=255)
            ggTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCt                                                                                                                                           <  1:650658/150‑1 (MQ=255)
             gTAATGGCAGTCCTTCGCCATACGCAATGGCATTAGACTGACCAGCAACGGTAAAAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTg                                                                                                                                          <  2:1118998/150‑1 (MQ=255)
              tAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCGGTAAAGGGGCCCTgc                                                                                                                                         <  1:150791/150‑2 (MQ=255)
               aaTGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCGTTACCGGGGCCCTGag                                                                                                                                        >  2:1584695/1‑149 (MQ=255)
               aaTGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCGGTAAAGGGGCCCTgct                                                                                                                                        <  1:602443/150‑3 (MQ=255)
               aaTGGCAGTCCTTCGCCATACGCAATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGAAATTGAAATGACAAATTGATAGCAGGGCCCCGGTAAAGGGGCCCTgct                                                                                                                                        <  2:1305612/150‑3 (MQ=255)
                    cAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGtt                                                                                                                                   >  2:124165/1‑150 (MQ=255)
                      gTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAg                                                                                                                                 >  2:317195/1‑150 (MQ=255)
                             gCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTa                                                                                                                          >  2:903563/1‑150 (MQ=255)
                                aTACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACtg                                                                                                                       >  2:889834/1‑150 (MQ=255)
                                aTACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTAGCGTAATTGATATGACAAATTGATAGCAGGGACCCTTTACAGGGGCCCTGAAGCGTTAGTTTTGTACtg                                                                                                                       >  1:1029147/1‑150 (MQ=255)
                                    gCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGAc                                                                                                                   >  1:1151459/1‑150 (MQ=255)
                                     ccATGGCATTAGACTGACCAGCAACGGTAAGAAAATAGTAATAATCGGGCGATATTATTGCGTTCATATATATATAGCGTAAATGAAAAGACAAATTGATAGCAGGGCACCTTTACCGGGGCCCTGAAGCGTTAGTTTTTAACTGTGACt                                                                                                                  >  1:315493/1‑150 (MQ=255)
                                       aTGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACttt                                                                                                                >  1:420354/1‑150 (MQ=255)
                                        gggCTTTAGAGTGACCAGCAACGGAAATAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCTTACTTGAAATGAAAAATTGATAGCAGGTCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTtta                                                                                                               <  2:236446/149‑1 (MQ=255)
                                           cATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGAGTAATTGAAATGACAAATTGATAGCATGTCCCCTTTACCGGGGCCCTGAAGCGTTAGATTTGTTCTGTGACTttatta                                                                                                            >  2:1342762/1‑150 (MQ=255)
                                             ttATACTGACCAGCAACTGTAAGAACATAGTAATAATCGGCCGATATTATTGCGTTCATATAAAATTCGCGTTATTGTAATGATAAATTCATAGCCGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATc                                                                                                          <  2:277529/150‑1 (MQ=255)
                                             ttAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATc                                                                                                          <  2:762873/150‑1 (MQ=255)
                                             ttAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATc                                                                                                          <  2:1622826/150‑1 (MQ=255)
                                              tAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTGTACCGGGGCCCTGAAGCGTTAGTTTTGTTCTGTGTCTTTATTATcc                                                                                                         >  1:1172505/1‑150 (MQ=255)
                                              tAGACTGACCAGCAACGGTAACAACAAAGTAATAATCGGGAGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATcc                                                                                                         <  1:136358/150‑1 (MQ=255)
                                               aGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCATGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCt                                                                                                        >  2:948827/1‑150 (MQ=255)
                                                  cTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTca                                                                                                     >  1:1472980/1‑150 (MQ=255)
                                                  cTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCACTCAAGCGTTAGTTTTGTACTGTGACTTTAATATCCTTca                                                                                                     >  1:415279/1‑150 (MQ=255)
                                                   tGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTcac                                                                                                    <  2:1137370/150‑1 (MQ=255)
                                                    gACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTcaca                                                                                                   >  1:593350/1‑150 (MQ=255)
                                                    gACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTACGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTcaca                                                                                                   >  1:979604/1‑150 (MQ=255)
                                                        aGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGATGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGa                                                                                               >  2:61924/1‑150 (MQ=255)
                                                         gCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGTCCCTGAAGCGGTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGAt                                                                                              <  1:209666/150‑1 (MQ=255)
                                                            aCGGTAAGAACATAGTAATAATCGGGAGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATAACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCaa                                                                                           <  2:1596817/150‑1 (MQ=255)
                                                                 aaGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTg                                                                                      <  2:354099/150‑1 (MQ=255)
                                                                 aaGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTg                                                                                      <  1:1004826/150‑1 (MQ=255)
                                                                   gAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAATCTGaa                                                                                    >  1:1675938/1‑150 (MQ=255)
                                                                   gAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGaa                                                                                    >  1:496767/1‑150 (MQ=255)
                                                                   gAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGaa                                                                                    >  1:739100/1‑150 (MQ=255)
                                                                     aCATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATc                                                                                  <  1:1730896/150‑1 (MQ=255)
                                                                      cATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCATTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATcc                                                                                 >  2:554735/1‑150 (MQ=255)
                                                                          gTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCa                                                                             >  1:1499919/1‑150 (MQ=255)
                                                                          gTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCa                                                                             >  1:1592748/1‑150 (MQ=255)
                                                                               aatCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCt                                                                        >  2:630344/1‑150 (MQ=255)
                                                                                  cGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTg                                                                     <  1:976893/150‑1 (MQ=255)
                                                                                  cGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTg                                                                     <  1:336442/150‑1 (MQ=255)
                                                                                  cGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTg                                                                     <  1:99885/150‑1 (MQ=255)
                                                                                       gATATTATAGCGTTCATATAAATTTAGCGTAATTGAAATGACAAATTGATAGCAGGGCACCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTgccgcc                                                                <  2:439483/150‑1 (MQ=255)
                                                                                        atatTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCt                                                               <  2:1201147/150‑1 (MQ=255)
                                                                                          attattGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGc                                                             >  1:1324426/1‑150 (MQ=255)
                                                                                            tattGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAg                                                           <  2:1522847/150‑1 (MQ=255)
                                                                                             attGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGt                                                          >  1:1067724/1‑150 (MQ=255)
                                                                                              ttGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGtt                                                         <  2:1400412/150‑1 (MQ=255)
                                                                                                gCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTc                                                       >  2:230752/1‑150 (MQ=255)
                                                                                                 cGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCa                                                      >  2:183228/1‑150 (MQ=255)
                                                                                                 cGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCa                                                      >  2:587843/1‑150 (MQ=255)
                                                                                                 cGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCa                                                      >  2:580429/1‑150 (MQ=255)
                                                                                                       tataAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGTGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCa                                                >  1:324619/1‑150 (MQ=255)
                                                                                                       tataAATTTCGCGTAATTGAAATGACAAATTGATAGCAGCGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCa                                                <  1:1097264/150‑1 (MQ=255)
                                                                                                         taAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCaat                                              <  2:420354/150‑1 (MQ=255)
                                                                                                         taAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCAGCAGTTTCACCGCCaat                                              >  2:179871/1‑150 (MQ=255)
                                                                                                            aTTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCaataat                                           <  2:1016792/150‑1 (MQ=255)
                                                                                                                cgcgTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGTGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAAt                                       <  1:1048161/150‑1 (MQ=255)
                                                                                                                cgcgTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAAt                                       <  2:1489298/150‑1 (MQ=255)
                                                                                                                  cgTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAg                                     <  1:1024556/150‑1 (MQ=255)
                                                                                                                    tAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGAc                                   >  1:460190/1‑150 (MQ=255)
                                                                                                                     aaTTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACt                                  <  2:1189528/150‑1 (MQ=255)
                                                                                                                        tGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACtatt                               >  1:1054865/1‑150 (MQ=255)
                                                                                                                        tGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACtatt                               >  1:111610/1‑150 (MQ=255)
                                                                                                                        tGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACtatt                               >  2:1144816/1‑150 (MQ=255)
                                                                                                                         gAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACtatta                              <  2:979604/150‑1 (MQ=255)
                                                                                                                         gAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACtatta                              >  2:1346975/1‑150 (MQ=255)
                                                                                                                            aTGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTc                           >  2:124930/1‑150 (MQ=255)
                                                                                                                              gACAAATTGATAGCAGGGCACCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCTAAAGTGAATCCGTCACCGCTTTACCGCCTGCAGATTCACCGCCAATAATCAATAGACTATTATTCCa                         <  2:969283/150‑1 (MQ=255)
                                                                                                                                cAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAgg                       >  2:1209677/1‑150 (MQ=255)
                                                                                                                                cAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAgg                       >  2:448135/1‑150 (MQ=255)
                                                                                                                                cAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAgg                       >  2:480606/1‑150 (MQ=255)
                                                                                                                                 aaaTTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGATTCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAggg                      >  1:147055/1‑150 (MQ=255)
                                                                                                                                  aaTTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTTACCTTCACAGTGATAAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAAAAAACAATAGACTATTATTCCAGGGc                     >  2:1563994/1‑150 (MQ=255)
                                                                                                                                      gATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATg                 >  1:769962/1‑150 (MQ=255)
                                                                                                                                        tAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGCCTTTATTATCCTTCACAGTGCTCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCCCCGCCAATAAGCAATAGACTATTCTTCCAGGGCAGTGAt               >  1:1211131/1‑150 (MQ=255)
                                                                                                                                        tAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGAt               >  1:522046/1‑150 (MQ=255)
                                                                                                                                        cttcAGGGCCCCGGTAAAGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACt                                  <  2:1593252/128‑1 (MQ=255)
                                                                                                                                        cttcAGGGCCCCGGTAAAGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACt                                  >  1:1593252/4‑131 (MQ=255)
                                                                                                                                        catcAGGGCCCCGGTAAAGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAAAAGACTATAATTCCAGGGCAATGAt               >  2:760399/4‑150 (MQ=255)
                                                                                                                                          gCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTTGTTTTGTACTGTGACTTTTTTTTCCTTCACTGTGATCAATTCTGAATCCGTCGCCGCTTTGCCGCCTGCAGTTTCCCCGCCAATAATCTATAGGCTTTTATTCCAGGGCAATGATAc             >  1:209728/1‑150 (MQ=255)
                                                                                                                                          gCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATAc             >  1:125630/1‑150 (MQ=255)
                                                                                                                                            aGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTc           <  2:557705/150‑1 (MQ=255)
                                                                                                                                            aGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTc           <  1:903563/150‑1 (MQ=255)
                                                                                                                                             gggCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTcc          >  1:1403491/1‑150 (MQ=255)
                                                                                                                                              ggCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCCCCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCg         >  1:1016350/1‑150 (MQ=255)
                                                                                                                                              ggCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCg         <  1:497825/150‑1 (MQ=255)
                                                                                                                                              ggCCCCTTTACCGGGGCCCTGAAGAGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCAAGGGCAATGATACTCCg         >  1:912845/1‑150 (MQ=255)
                                                                                                                                                  ccTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAgg     >  2:1454240/1‑150 (MQ=255)
                                                                                                                                                  ccTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAgg     >  2:1711144/1‑150 (MQ=255)
                                                                                                                                                    tttACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGcc   <  1:545219/150‑1 (MQ=255)
                                                                                                                                                     ttACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGccc  <  1:463578/150‑1 (MQ=255)
                                                                                                                                                     ttACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGccc  <  1:186581/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCC  >  NC_000913/4537491‑4537789

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: