Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,537,639 T→G 5.3% intergenic (‑45/+20) 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,6390TG 5.3% 350.8 / 4.8 117intergenic (‑45/+20)nanS/nanM9‑O‑acetyl N‑acetylneuraminic acid esterase/N‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base T (54/56);  new base G (3/4);  total (57/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

TCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCC  >  NC_000913/4537491‑4537788
                                                                                                                                                    |                                                                                                                                                     
tCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGCCCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTGCATATAAATTTCGCGTAATTGAAAGGACAAATTGATAGCAGGGCCCCtt                                                                                                                                                      >  2:202565/1‑150 (MQ=255)
      cTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCgg                                                                                                                                                >  1:1330658/1‑150 (MQ=255)
         tCCGGTAATGGCAGTCATTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTATGAACTTAGTAATAATCTGTCGATATTTTTGCGTTCATATAAATTTCGCGTATTTGAAATGACAAATTGATATCAGGGCCCCTTTACCGGGGc                                                                                                                                             <  1:752410/150‑1 (MQ=255)
         gccGGTAATGGCATTCCTTCGCCATACGCAATGGCATTAGACTGACCAGCAACGGTAAGATCATATTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGc                                                                                                                                             <  2:938619/149‑1 (MQ=255)
           cGGTAATGGCAGTCCTTCGCCATACGCCATTGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGccc                                                                                                                                           <  1:557150/150‑1 (MQ=255)
           cGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGATATGACAAATTGATAGCAGGGCCCCTTTACCGGGGccc                                                                                                                                           >  2:305552/1‑150 (MQ=255)
              tAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCGGTAAAGGGGCCCTgc                                                                                                                                        <  2:1519077/150‑2 (MQ=255)
                  ggCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGATGCg                                                                                                                                    >  1:285003/1‑150 (MQ=255)
                  ggCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCg                                                                                                                                    >  2:1691646/1‑150 (MQ=255)
                  ggCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCg                                                                                                                                    >  2:1706064/1‑150 (MQ=255)
                  ggCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCg                                                                                                                                    <  1:1402771/150‑1 (MQ=255)
                  ggCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCGGTAAAGGGGCCCTgctatc                                                                                                                                    <  2:1415703/150‑6 (MQ=255)
                  ggCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAACAAAATAGTAATAATCGGGCTATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCg                                                                                                                                    <  2:1157603/150‑1 (MQ=255)
                   gCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTACTTGAAATGACAAATTGATCGCAGGGCCCCGGTAACGGGGCCCTgctatca                                                                                                                                   <  1:1019229/150‑7 (MQ=255)
                   gCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGt                                                                                                                                   <  2:808850/150‑1 (MQ=255)
                      gTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTCTTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAg                                                                                                                                <  1:1005398/150‑1 (MQ=255)
                       tCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAATAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGt                                                                                                                               <  2:1573611/150‑1 (MQ=255)
                        ccTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGtt                                                                                                                              >  2:455067/1‑150 (MQ=255)
                        acTTCGACAAACGCCATGGCATTAGACTGACCAGCAACTGTAAGAACATACTAATAATCGGGCGATACTATTGCGATCATATACATTTCGCGTAACTGATATGACTAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGtt                                                                                                                              <  1:81120/149‑1 (MQ=255)
                         cTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGTGCCCCTTTACCGGGGCCCTGAAGCGTTAGttt                                                                                                                             >  1:1151512/1‑150 (MQ=255)
                         cTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCACCTTTACCGGGGACCTGAAGCGTTAGttt                                                                                                                             >  1:1840333/1‑150 (MQ=255)
                         cTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGAGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGttt                                                                                                                             >  1:1258172/1‑150 (MQ=255)
                           tCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATACTCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTg                                                                                                                           <  1:336666/150‑1 (MQ=255)
                            cgacATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATAATTGCGTTCATATAAATTACGCGTAATTGAAATGACAAATTGATAGCAGGGCAACTTTAACGGTGCCCTGAAGCGATAGTATTGt                                                                                                                          >  1:462767/4‑150 (MQ=255)
                            cGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGt                                                                                                                          >  1:462759/1‑150 (MQ=255)
                               cATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGATGCGTTAGTTTTGTagt                                                                                                                       >  1:296848/1‑148 (MQ=255)
                                aTACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACtt                                                                                                                      >  1:1578254/1‑149 (MQ=255)
                                aTACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACtg                                                                                                                      <  1:1818581/150‑1 (MQ=255)
                                  aCGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACtgtg                                                                                                                    <  2:1433747/150‑1 (MQ=255)
                                    gCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGAc                                                                                                                  <  1:1401422/150‑1 (MQ=255)
                                    gCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGAc                                                                                                                  <  1:1576232/150‑1 (MQ=255)
                                              tAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATcc                                                                                                        <  1:1802051/150‑1 (MQ=255)
                                              tAGACTGACCAGCAACGGGAAGAACCTAGTAATAATCGTGCGATTTTATTGCGGTCATATAAATTTCGCGTAATTGAATTGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATcc                                                                                                        <  1:265820/150‑1 (MQ=255)
                                                  cTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTca                                                                                                    >  1:52694/1‑150 (MQ=255)
                                                   tGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTcac                                                                                                   >  1:566997/1‑150 (MQ=255)
                                                   tGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGAGAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTcac                                                                                                   <  2:131375/150‑1 (MQ=255)
                                                    gACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTcaca                                                                                                  >  1:1610401/1‑150 (MQ=255)
                                                        aGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGTCAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGa                                                                                              <  1:966122/150‑1 (MQ=255)
                                                        aGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGa                                                                                              <  1:1519871/150‑1 (MQ=255)
                                                         gCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGAt                                                                                             <  1:1099295/150‑1 (MQ=255)
                                                          cAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATc                                                                                            <  1:1375055/150‑1 (MQ=255)
                                                           aaCGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCATGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTTATCa                                                                                           >  1:1626623/1‑150 (MQ=255)
                                                              ggTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCaaaa                                                                                        >  2:503948/1‑150 (MQ=255)
                                                               gTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAAc                                                                                       >  2:688795/1‑150 (MQ=255)
                                                                tAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACt                                                                                      >  2:537714/1‑150 (MQ=255)
                                                                tAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACt                                                                                      <  1:1727357/150‑1 (MQ=255)
                                                                 aaGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTg                                                                                     <  2:1817476/150‑1 (MQ=255)
                                                                 aaGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTg                                                                                     <  2:1363317/150‑1 (MQ=255)
                                                                     aCATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATc                                                                                 >  1:892918/1‑150 (MQ=255)
                                                                      cATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATcc                                                                                <  2:1656636/150‑1 (MQ=255)
                                                                        tAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGt                                                                              >  1:937800/1‑150 (MQ=255)
                                                                         aGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTc                                                                             <  2:1261833/150‑1 (MQ=255)
                                                                           taataaTCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCAc                                                                           >  1:1725972/1‑150 (MQ=255)
                                                                            aataatCGGGCGATATTATTGCGTTCATATAAATTTCGAGTAATTGAAATGACAAATTGATATCAGGGCCCCTTTACCGTGTCCCTGAAGCGTTAGTTTTGTACTTTTACTTTATTATCCTTCACAGTGATAAAAACTGAATCCGTAAcc                                                                          >  1:1754460/1‑150 (MQ=255)
                                                                                       gATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTgccgcc                                                               >  1:784471/1‑150 (MQ=255)
                                                                                          attattGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGc                                                            <  2:1280098/150‑1 (MQ=255)
                                                                                             attGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGt                                                         >  1:1394468/1‑150 (MQ=255)
                                                                                              ttGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGtt                                                        <  1:886059/150‑1 (MQ=255)
                                                                                              ttGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGtt                                                        <  1:458380/150‑1 (MQ=255)
                                                                                                 cgagcATATAAATTTCTCGTAATTGAAATGACAAATTGATAGCCGGGCCCCTTTACCGGGGCCCTGAAGAGGTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCa                                                     <  1:305552/146‑1 (MQ=255)
                                                                                                      atatAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGcc                                                <  1:505797/150‑1 (MQ=255)
                                                                                                       tataAATTTCGCGTAATTGATATGACAAATAGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGAGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACGGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCa                                               >  1:227980/1‑150 (MQ=255)
                                                                                                          aaaTTTCGCGTAATTGAAATGACAAATTGATAGCAGGTCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCCCAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCaata                                            <  2:898683/150‑1 (MQ=255)
                                                                                                             tttCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATc                                         <  1:1550851/150‑1 (MQ=255)
                                                                                                             tttCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATc                                         <  1:1184254/150‑1 (MQ=255)
                                                                                                              ttCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCa                                        <  1:1348756/150‑1 (MQ=255)
                                                                                                                cgcgTAATTGAAATGACTAATAGATAGCAGGGCCCCTTAACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAAt                                      <  1:283624/150‑1 (MQ=255)
                                                                                                                cgcgTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAAt                                      >  2:1675/1‑150 (MQ=255)
                                                                                                                 gcgTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATa                                     >  1:1589132/1‑150 (MQ=255)
                                                                                                                         gAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACtatta                             >  2:1637480/1‑150 (MQ=255)
                                                                                                                         gAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGATTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGATTTGCCGCCTGCAGTTTCACCGCCAAAAATAAATAGACtataa                             >  2:232068/1‑148 (MQ=255)
                                                                                                                         gAAATGACAAATTGATAGCACGGCCCCTATACCGGGGCCCGGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACtatta                             <  1:131284/150‑1 (MQ=255)
                                                                                                                              gACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCa                        <  1:503948/150‑1 (MQ=255)
                                                                                                                               aCAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAg                       <  1:711297/150‑1 (MQ=255)
                                                                                                                               aCAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAg                       <  2:720785/150‑1 (MQ=255)
                                                                                                                                cAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAgg                      >  2:1071067/1‑150 (MQ=255)
                                                                                                                                cAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAgg                      >  2:1037902/1‑150 (MQ=255)
                                                                                                                                cAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGACTGCAGTTTCACCGCCAATAATCAATAGACTATTATTACAgg                      >  1:420062/1‑150 (MQ=255)
                                                                                                                                cAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACAGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAgg                      >  1:93094/1‑150 (MQ=255)
                                                                                                                                cAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTATTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAgg                      >  1:116599/1‑150 (MQ=255)
                                                                                                                                 aaaTTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAggg                     <  2:1313698/150‑1 (MQ=255)
                                                                                                                                 aaaTTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAggg                     >  2:419112/1‑150 (MQ=255)
                                                                                                                                  aaTTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGc                    <  2:1617728/150‑1 (MQ=255)
                                                                                                                                    ttGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCaa                  >  1:964523/1‑150 (MQ=255)
                                                                                                                                      gATAGCAGGGCCCCTTTACCGGGGCCCTGATGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCTTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCaagg                >  2:1020003/1‑148 (MQ=255)
                                                                                                                                       atatCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTAATATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGa               <  2:1191528/150‑1 (MQ=255)
                                                                                                                                       aTAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGa               <  1:66845/150‑1 (MQ=255)
                                                                                                                                        cttcAGGGCCCCGGTAAAGGGGCCCTGAAGAGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACAGAATCCGTCACCGCTTTGCCGCCTGCAGTTTAACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGAt              >  1:432639/4‑150 (MQ=255)
                                                                                                                                         aGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATa             >  1:977548/1‑150 (MQ=255)
                                                                                                                                         aGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATa             >  1:26029/1‑150 (MQ=255)
                                                                                                                                          tcAGGGCCCCGGTAAAGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTTTTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGACGCCTGCAGTTTCCCCGCCAATAATCAATAGACGATTATTCCCGGGCAATGATAc            >  2:1308721/2‑150 (MQ=255)
                                                                                                                                           cAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGCATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACt           >  2:458580/1‑150 (MQ=255)
                                                                                                                                            aGGGCCCCGGTAAAGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACta          >  1:141523/1‑149 (MQ=255)
                                                                                                                                              ggCCCCGTTAGCGGGGCCCTGAAGCGTTAGTGTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCg        <  2:749751/150‑1 (MQ=255)
                                                                                                                                               gCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCGCCGCCAATAATCAATAGACTATTATTCCAGGGCTATGATACTCCGt       >  2:532623/1‑150 (MQ=255)
                                                                                                                                                ccccTTTACCGGGGCCCTGAAGCGTTGGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTa      <  2:1347822/150‑1 (MQ=255)
                                                                                                                                                ccccTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTa      <  2:565469/150‑1 (MQ=255)
                                                                                                                                                  ccTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAgg    >  1:488271/1‑150 (MQ=255)
                                                                                                                                                   cTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGc   >  1:765721/1‑150 (MQ=255)
                                                                                                                                                    tttACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGAcc  <  1:1496476/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCC  >  NC_000913/4537491‑4537788

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