Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,537,643 C→A 5.3% intergenic (‑49/+16) 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,6430CA 5.3% 413.0 / 9.6 132intergenic (‑49/+16)nanS/nanM9‑O‑acetyl N‑acetylneuraminic acid esterase/N‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base C (54/71);  new base A (2/5);  total (56/76)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-01

Only 100 of 134 total aligned reads displayed.

CCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGA  >  NC_000913/4537495‑4537791
                                                                                                                                                    |                                                                                                                                                    
cccTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTAcc                                                                                                                                                     >  2:1515669/1‑150 (MQ=255)
  cTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCgg                                                                                                                                                   <  1:225665/150‑1 (MQ=255)
    gTCCGGTAATGGCAGTCCTTCGCCTTACGCCATGGCATTAGACTGACCAGCAAAGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCgggg                                                                                                                                                 <  2:1438375/150‑1 (MQ=255)
     tactgTAATTTCAGGCCTTCGCAATACGCCATTGCATTAGACTGACCAGCAACTGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGc                                                                                                                                                <  2:189417/146‑1 (MQ=255)
     tCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACTGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGc                                                                                                                                                <  2:804835/150‑1 (MQ=255)
       cGGTAATGGCAGTCCTTCGCCATACTCCATGGCATGAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCATGGCCCCTTTACCGGGGccc                                                                                                                                              <  1:1465005/150‑1 (MQ=255)
       cGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCGGTAAAGGGGccc                                                                                                                                              <  2:742287/150‑1 (MQ=255)
       cGGTAATGGCACTCCTTCGCCCTACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGccc                                                                                                                                              <  2:1011716/150‑1 (MQ=255)
        ggTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCt                                                                                                                                             >  2:1011288/1‑150 (MQ=255)
        ggTAATGGCAGTCCTTCGCCATACACCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCt                                                                                                                                             <  1:1143622/150‑1 (MQ=255)
          tAATGGCAGTCCTTCGCCATACGCGATGGCATTAGACTGAGCAGCTACGGTAAGAAGATAGTTATAATCGGGGGATATTTTTGCGTTCATATAAATTTCGCGTAATTTAAATGACAAATTGATAGCAGGGCCCCGGTAAAGGGGCCCTgc                                                                                                                                           <  2:893265/150‑2 (MQ=255)
          tAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGa                                                                                                                                           <  1:502947/150‑1 (MQ=255)
          tAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGa                                                                                                                                           <  1:1366268/150‑1 (MQ=255)
          tAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCGGTAAAGGGGCCCTgc                                                                                                                                           <  2:1081302/150‑2 (MQ=255)
          tAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCGGTAAAGGGGCCCTgc                                                                                                                                           <  2:901143/150‑2 (MQ=255)
           aaTGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATACTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCGGTAAAGGGGCCCTgct                                                                                                                                          <  2:179452/150‑3 (MQ=255)
            aTGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAg                                                                                                                                         <  1:930873/150‑1 (MQ=255)
              ggCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCg                                                                                                                                       <  2:888819/150‑1 (MQ=255)
              ggCAGTCCTTCGCCATACGCCATGACATTATACTGACCAGCAAGGGTAAGAACATAGTAATAATCGTGCGATATTATTGCGTTCATATAAATATCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCg                                                                                                                                       <  1:1333577/150‑1 (MQ=255)
                cAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGtt                                                                                                                                     >  1:260369/1‑150 (MQ=255)
                     cTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGttt                                                                                                                                <  2:1476920/150‑1 (MQ=255)
                     cTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGttt                                                                                                                                >  1:1165171/1‑150 (MQ=255)
                     cTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGttt                                                                                                                                <  1:542411/150‑1 (MQ=255)
                       tCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTg                                                                                                                              <  2:1778273/150‑1 (MQ=255)
                        cGCCATACGCCATGTCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAATTGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGt                                                                                                                             <  1:2009668/150‑1 (MQ=255)
                         gcaaTACGCCATGTCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTa                                                                                                                            <  2:592050/147‑1 (MQ=255)
                          ccATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTCTTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTAc                                                                                                                           <  1:277915/150‑1 (MQ=255)
                           cATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACt                                                                                                                          >  1:680767/1‑150 (MQ=255)
                              aCGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACtgtg                                                                                                                       <  2:520366/150‑1 (MQ=255)
                               cGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGa                                                                                                                      >  1:826195/1‑150 (MQ=255)
                                gCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGAc                                                                                                                     >  1:562641/1‑150 (MQ=255)
                                 ccATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACt                                                                                                                    <  2:45859/150‑1 (MQ=255)
                                 ccATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACt                                                                                                                    <  2:1331896/150‑1 (MQ=255)
                                  cATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACtt                                                                                                                   >  2:1330786/1‑150 (MQ=255)
                                       cATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCACTGAAGCGTTAGTTTTGTACTGTGACTttatta                                                                                                              >  1:232829/1‑150 (MQ=255)
                                        aTTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTtattat                                                                                                             >  1:1012167/1‑150 (MQ=255)
                                          tAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATcc                                                                                                           <  1:1185913/150‑1 (MQ=255)
                                            gACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCtt                                                                                                         <  2:815294/150‑1 (MQ=255)
                                             aCTGACAAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACGGTGACTTTATTATCCTTc                                                                                                        >  2:490763/1‑150 (MQ=255)
                                              cTGCCCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAAATTCGCGTAATTGAAATGACAAATTGATATCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTTGTTTTGTACTGTGACTTTATTATCCTTca                                                                                                       <  2:129590/150‑1 (MQ=255)
                                              cTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTca                                                                                                       >  1:413639/1‑150 (MQ=255)
                                              cTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTca                                                                                                       >  2:6022/1‑150 (MQ=255)
                                                gACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTcaca                                                                                                     >  2:1201785/1‑150 (MQ=255)
                                                    aGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGa                                                                                                 >  2:1675571/1‑150 (MQ=255)
                                                    aGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGa                                                                                                 >  1:1520831/1‑150 (MQ=255)
                                                     gCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGAt                                                                                                >  1:1980291/1‑150 (MQ=255)
                                                        aCGTTAAGAACATAGTAATAATCTGGCGATATTATTTCTTTCATATAAATTTCGAGAAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCat                                                                                             <  1:1225751/150‑2 (MQ=255)
                                                          ggTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCaaaa                                                                                           >  1:1736696/1‑150 (MQ=255)
                                                           gTACGAACATAGTAAAAATCTGGCGCTATTATAGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAAc                                                                                          <  2:289515/150‑1 (MQ=255)
                                                           gTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAAc                                                                                          <  1:1068429/150‑1 (MQ=255)
                                                              aGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGa                                                                                       >  1:551502/1‑150 (MQ=255)
                                                              aGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGa                                                                                       >  1:462092/1‑150 (MQ=255)
                                                                  cATAGTAATAATCGGGCGATATTATTGCGTTCATAAAAATTTCGCGTAAATGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCATGAAGCGTTAGTTTTGTACTGTGACTTTATTATCATTCACAGTGATCAAAACTGAATcc                                                                                   >  2:323042/1‑150 (MQ=255)
                                                                  aaTAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATcc                                                                                   <  1:1737940/149‑1 (MQ=255)
                                                                     aGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTc                                                                                <  2:232771/150‑1 (MQ=255)
                                                                      gTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCa                                                                               >  2:774342/1‑150 (MQ=255)
                                                                         ataatCGGGCGATATTATTGCTTTCATATAAATTTCTCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACAGGGGCCCTGAAGCGTTAGTTTTGTAATGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCAGTCAcag                                                                            >  2:207675/1‑148 (MQ=255)
                                                                            atCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCtt                                                                         >  1:151318/1‑150 (MQ=255)
                                                                              cGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTg                                                                       <  1:1790218/150‑1 (MQ=255)
                                                                                   gATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTgccgcc                                                                  <  2:955009/150‑1 (MQ=255)
                                                                                      attattGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGc                                                               <  2:629266/150‑1 (MQ=255)
                                                                                       ttattGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCa                                                              <  1:1056862/150‑1 (MQ=255)
                                                                                         attGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGt                                                            <  1:1916229/150‑1 (MQ=255)
                                                                                          ttGCGTTCATATAAATTTCTCTAAATTGAAATGACAAATTGATAGCAGCGCCCCTTTACCGGGGCCCTGAAGCGATAGTGTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGGCACCGCTTTGCCGCCTGCAGtt                                                           <  2:139418/150‑1 (MQ=255)
                                                                                          ttGCGTTCATATAAATTTCGCGTAATTGAAAAGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGtt                                                           <  1:867185/150‑1 (MQ=255)
                                                                                           tGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGttt                                                          <  2:1130740/150‑1 (MQ=255)
                                                                                             cGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCa                                                        <  2:1392424/150‑1 (MQ=255)
                                                                                               ttCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAAATGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCAcc                                                      >  1:1322361/1‑150 (MQ=255)
                                                                                                tCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCg                                                     <  2:945530/150‑1 (MQ=255)
                                                                                                tCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCg                                                     <  1:976641/150‑1 (MQ=255)
                                                                                                 cATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGc                                                    >  2:1150856/1‑150 (MQ=255)
                                                                                                  atatAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGCCTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGcc                                                   >  1:367609/1‑150 (MQ=255)
                                                                                                  atatAAAATTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCTGGGCCCTGAAGCGTTAGTTATGTACTGTGACTTTATTATCCTACACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCTcc                                                   >  1:362149/1‑150 (MQ=255)
                                                                                                   tataAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCa                                                  <  1:1734379/150‑1 (MQ=255)
                                                                                                   tataAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGAACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCa                                                  >  1:1853634/1‑150 (MQ=255)
                                                                                                    ataAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCaa                                                 >  2:1648497/1‑150 (MQ=255)
                                                                                                     taATTTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGATTTGTACTGTGACTTAATTATCCTTCACAGAGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCaaa                                                >  2:1451423/1‑149 (MQ=255)
                                                                                                       aaTTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCaataa                                              >  1:58401/1‑150 (MQ=255)
                                                                                                        aTTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCaataat                                             >  1:1551320/1‑150 (MQ=255)
                                                                                                         tttCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATc                                            <  1:1206750/150‑1 (MQ=255)
                                                                                                          ttCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCa                                           <  2:1992258/150‑1 (MQ=255)
                                                                                                             gcgTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATa                                        <  1:777655/150‑1 (MQ=255)
                                                                                                              cgTAATTGTAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAATCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAg                                       <  2:177739/150‑1 (MQ=255)
                                                                                                                 aaTTGAAATGTCAAATTTTTAGCAGGGCCCCTTTCCCGGCGCCCTGAAGCGTTAGTTTTGTCCTGTGACTTTATTCTCCTTCACAGTGATCAAAACTGAATCCGTCACGGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACt                                    <  2:112951/150‑1 (MQ=255)
                                                                                                                 aaTTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACt                                    <  2:1257620/150‑1 (MQ=255)
                                                                                                                     gAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGGTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACtatta                                >  2:1542774/1‑150 (MQ=255)
                                                                                                                     gAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACtatta                                <  2:1520831/150‑1 (MQ=255)
                                                                                                                         tGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTcc                            >  1:1204122/1‑150 (MQ=255)
                                                                                                                           aCAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAg                          <  1:1609738/150‑1 (MQ=255)
                                                                                                                            cAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAgg                         >  1:1807576/1‑150 (MQ=255)
                                                                                                                               aTTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTACTATTCCAGGGCa                      >  1:1873621/1‑150 (MQ=255)
                                                                                                                                ttGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCaa                     >  2:655277/1‑150 (MQ=255)
                                                                                                                                   aTAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGa                  <  1:1995688/150‑1 (MQ=255)
                                                                                                                                    cttcAGGGCCCCGGTAAAGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGAt                 >  1:694711/4‑150 (MQ=255)
                                                                                                                                        aGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGGATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTc             <  2:656088/150‑1 (MQ=255)
                                                                                                                                         gggCCCCGGTAAAGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTcc            >  2:1634804/1‑150 (MQ=255)
                                                                                                                                             cccTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAg        <  1:1330786/150‑1 (MQ=255)
                                                                                                                                             cccTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACAGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAg        >  2:1452801/1‑150 (MQ=255)
                                                                                                                                              ccTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTAACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAgg       >  2:1360497/1‑150 (MQ=255)
                                                                                                                                                   aCCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGACCGa  >  1:1786817/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGATATTATTGCGTTCATATAAATTTCGCGTAATTGAAATGACAAATTGATAGCAGGGCCCCTTTACCGGGGCCCTGAAGCGTTAGTTTTGTACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTGCCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGA  >  NC_000913/4537495‑4537791

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: