Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,159,335 G→T 6.0% intergenic (+33/+27) ptsG → / ← fhuE fused glucose‑specific PTS enzymes: IIB component/IIC component/ferric‑rhodotorulic acid outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,159,3350GT 6.0% 268.2 / 5.0 108intergenic (+33/+27)ptsG/fhuEfused glucose‑specific PTS enzymes: IIB component/IIC component/ferric‑rhodotorulic acid outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (55/46);  new base T (3/4);  total (58/50)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

Only 100 of 109 total aligned reads displayed.

TGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAG  >  NC_000913/1159187‑1159482
                                                                                                                                                    |                                                                                                                                                   
tGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAg                                                                                                                                                     >  2:1369033/1‑149 (MQ=255)
  aagaaACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGAt                                                                                                                                                   >  1:1215041/1‑149 (MQ=255)
  aagaaACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGAt                                                                                                                                                   <  1:255833/149‑1 (MQ=255)
   agaaACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTATTCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATg                                                                                                                                                  >  2:82229/1‑149 (MQ=255)
       aCTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATAAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTg                                                                                                                                              >  1:270032/1‑149 (MQ=255)
       aCTGGGCGCAGCGGGCGTAGTGGTGGCTGGTGCTGGTGTTCAGGCGATTTTCGGTACTAAATCCTAAAACCTGAAAACCGAGATCGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTg                                                                                                                                              <  2:503168/149‑1 (MQ=255)
         tGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTGCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGGTTGGCAGCCTAAGGCAGCCAGATCGCTGcg                                                                                                                                            >  2:1382987/1‑148 (MQ=255)
         ggggCGCAGCGGGCGTAGTGGTTGCTGGTGCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGcc                                                                                                                                            <  1:1222055/148‑1 (MQ=255)
            gcgcAGCGGGCGTGGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCttt                                                                                                                                         <  1:603768/149‑1 (MQ=255)
            gcgcAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCTATTTTCGGTACTAAATCCGATAACCTGAAAAGCGAGATGGATGAGTTCATCCGTAACCACTAATCCGTAAGTCGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCtta                                                                                                                                         <  1:1446783/149‑2 (MQ=255)
              gcAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGAAGGCTGCCttttt                                                                                                                                       >  1:1063650/1‑149 (MQ=255)
               cAGCGGGCGTAGTGGTTGCGTGGTCTGGTGTTCAGGCGAGATTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCttagtc                                                                                                                                      <  2:1083673/149‑5 (MQ=255)
                aGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTa                                                                                                                                     >  1:436869/1‑149 (MQ=255)
                 gCGGGCGTAGTGGTTGCTGGTTCTGGTGGTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCTCTAATCCGTAAGAGGTGGGGGAGACTAAGGCAGCCAGAAGCCTGCCTTTTTTAc                                                                                                                                    >  2:1031846/1‑149 (MQ=255)
                 ccGGGCGTAGTGGGTGCTGGGTCTGGTGGGCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGCGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTAc                                                                                                                                    <  1:1088971/148‑1 (MQ=255)
                  cGGTCGTATTGGTTGCTGGTTCTGGTTTTCAGGCGATTTTCGTTACTAAATCCGATAACCTGAAAACCGATATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGTGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACa                                                                                                                                   <  1:991570/149‑1 (MQ=255)
                  cGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACa                                                                                                                                   <  2:38143/149‑1 (MQ=255)
                  cGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCGGCCTTTTTTACa                                                                                                                                   >  2:874807/1‑149 (MQ=255)
                   gggCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTATCGGTACTAAATCCGATAACCTGAATACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTACGGCAGCCAGATGGCTGCCTTTTTTACAg                                                                                                                                  <  2:350269/149‑1 (MQ=255)
                   gggCGTAGTGGGTGCTGGGTCTGGTGTGCAGGCGATTTTCGGGACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAg                                                                                                                                  <  2:487750/149‑1 (MQ=255)
                       gTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGgtgt                                                                                                                              >  1:1220868/1‑149 (MQ=255)
                          gTGGTTGCTGGTTCTGTTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGCCGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTAt                                                                                                                           >  2:563159/1‑149 (MQ=255)
                           tGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTAtt                                                                                                                          <  2:1063383/149‑1 (MQ=255)
                           gggTTGCTGGTTCTGGTGTTCAGGCTATTTTCGGTACTAAATCCGATAAACTGAAAACCGAGATGTATGAGTACATCCGTAACCACTAATCCGTAAGAAGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTAtt                                                                                                                          <  2:1547072/148‑1 (MQ=255)
                             gtggcaggTTCTGGTGTTCCGGCGCTTTTCGGTACTAACACCGATAACCTGAAAACCCAGATGGATGAGAACATCCGTAACCACTACTCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCa                                                                                                                        <  2:287209/143‑1 (MQ=255)
                             gTTGCTGGTTCTGGTCTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGGACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCa                                                                                                                        <  1:1167154/149‑1 (MQ=255)
                                 ctggttctggttTTCAGGCTATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAAt                                                                                                                    <  2:1429645/149‑1 (MQ=255)
                                    gTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGa                                                                                                                 >  1:144708/1‑149 (MQ=255)
                                    gTTATGGTGATCATGCGATTTTCGTTACTAAATCCGATAAGCTGAAATCCGAGATGGATGAGTACATCCTTAACCACAAATCAGTAAGACGTATGGGAGACTAATGCAGCAAGATTCATAACTTATTTACATTAGTTATTCAGAAttca                                                                                                                 >  2:1212954/1‑147 (MQ=255)
                                           tgtTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGc                                                                                                          >  1:25655/1‑149 (MQ=255)
                                             gtCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCg                                                                                                        <  2:839989/148‑1 (MQ=255)
                                              tCAGGAGATATTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGTGGAGACTAAGGCAGACAGATGGATGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCgg                                                                                                       <  1:563159/149‑1 (MQ=255)
                                                  gCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAAt                                                                                                   >  1:409012/1‑149 (MQ=255)
                                                  gCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTAAGAATTGATACGTGCCGGTAAt                                                                                                   >  2:936743/1‑149 (MQ=255)
                                                     aTTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCt                                                                                                >  1:1463779/1‑149 (MQ=255)
                                                      ttttCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTg                                                                                               >  2:462649/1‑149 (MQ=255)
                                                         tCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTCCCGGTAATGCTGaaa                                                                                            >  1:700820/1‑149 (MQ=255)
                                                            gTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATCCTGAAATTa                                                                                         >  1:845361/1‑149 (MQ=255)
                                                                aaattCCGATAAAGTGAAAACCGAGATGGTTTAGTACATCAGTATCCACTAATCCGTAAGACGTAGTGGAGACTATGGCAGCCAGATGGCTGCCTTTTTTACCGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGATATTAcgcg                                                                                     <  1:1445154/145‑1 (MQ=255)
                                                                    tCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGgtgt                                                                                 >  1:894109/1‑149 (MQ=255)
                                                                      cGATAACCTGAAAACCGAGATGTATGAGTACATGCGTAACGACAAATCAGTAAGTCTTTGGGGAGACTAAGGCAGGCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGc                                                                               <  2:432114/149‑1 (MQ=255)
                                                                       gATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGcc                                                                              <  2:1220868/149‑1 (MQ=255)
                                                                       gATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGcc                                                                              >  1:1458952/1‑149 (MQ=255)
                                                                       gATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTATTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGcc                                                                              >  1:575811/1‑149 (MQ=255)
                                                                            ccTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCTAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTATGCGGTGTGCCGTAgc                                                                         >  2:195885/1‑148 (MQ=255)
                                                                            ccTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGa                                                                         <  1:502055/149‑1 (MQ=255)
                                                                             cTGAAAACCGAGATGGATGAGTACATCCGTAACAACTAATACGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCAGGTAATGCTGAAATTACGCGTTGTGCCGTAGaa                                                                        >  1:70876/1‑148 (MQ=255)
                                                                             cTGAAAACCGACATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGTGGAGACTAAGGCAGCCAGATGCCTGCCTTTATTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATCCTGAAATTACGAGGTGTGCTGTAGAc                                                                        >  1:567609/1‑149 (MQ=255)
                                                                               taaaaCCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGa                                                                      <  1:115230/148‑1 (MQ=255)
                                                                                 aaaCCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATa                                                                    <  1:767570/149‑1 (MQ=255)
                                                                                       tgatgtatgCGTACATCCGTAAGCACTCATCCGTAAGACGTTGGGGAGACTAAGGCAGCCCGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCt                                                              <  2:1132471/148‑1 (MQ=255)
                                                                                       agaTGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCt                                                              >  1:1306066/1‑149 (MQ=255)
                                                                                          tggatgATTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTcc                                                           <  2:951444/149‑1 (MQ=255)
                                                                                            gatgAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCGGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCAc                                                         >  2:539999/1‑149 (MQ=255)
                                                                                            gatgAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCAc                                                         >  2:1317820/1‑149 (MQ=255)
                                                                                             atgAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGATCCTTCCACg                                                        >  1:1225466/1‑149 (MQ=255)
                                                                                                aGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTg                                                     <  1:962577/149‑1 (MQ=255)
                                                                                                aGTACATCCGTAACAACTAATACGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCTGAATTGATACGTGCCGGTAATGATGAAATTACGCGGTGTGCCGTAGACGATATAACCTTCAAAGTTg                                                     >  2:603345/1‑149 (MQ=255)
                                                                                                 gTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTgg                                                    >  2:1426733/1‑149 (MQ=255)
                                                                                                  tACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACgggggt                                                   >  2:1561448/1‑144 (MQ=255)
                                                                                                  tACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGt                                                   <  1:362935/149‑1 (MQ=255)
                                                                                                       ccGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGt                                              >  1:495396/1‑149 (MQ=255)
                                                                                                         gTAACCACTAATCCGTTAGACGATGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAg                                            <  2:270032/149‑1 (MQ=255)
                                                                                                         gTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAg                                            >  1:1359998/1‑149 (MQ=255)
                                                                                                           aaCCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGt                                          >  2:1480898/1‑149 (MQ=255)
                                                                                                            aCCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGTTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGtt                                         <  1:82229/149‑1 (MQ=255)
                                                                                                            aCCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGtt                                         >  1:406247/1‑149 (MQ=255)
                                                                                                            aCCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGtt                                         <  1:390610/149‑1 (MQ=255)
                                                                                                              cACTACTCCGTAAGACGTTGGGGAGACTATGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCATAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGTTAGAACCTTCCACGTTGGTATCGTAGGtttt                                       <  2:47709/149‑1 (MQ=255)
                                                                                                              cACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGtttt                                       >  2:615288/1‑149 (MQ=255)
                                                                                                              cACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGGTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTGGGTATCGTAGGtttt                                       <  1:1126150/149‑1 (MQ=255)
                                                                                                              cACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTAAAGGTGTTAATCAGAATTGAAACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGattt                                       >  1:1155967/1‑149 (MQ=255)
                                                                                                               aCTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTg                                      <  2:25655/149‑1 (MQ=255)
                                                                                                                 tAATCCGTAAGACGTTGGGGAGTCTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGATAGTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTc                                    <  1:1031846/149‑1 (MQ=255)
                                                                                                                   aTCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGa                                  >  2:1437858/1‑149 (MQ=255)
                                                                                                                      cGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACa                               <  1:1369033/149‑1 (MQ=255)
                                                                                                                      cGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTAATTCAGAATTGATACGTTGCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACTTTAGAACCTTCCACGTTTGTATCGTAGGTTTTGTCGATCa                               <  2:1123879/149‑1 (MQ=255)
                                                                                                                        tAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAgg                             <  1:502526/149‑1 (MQ=255)
                                                                                                                          aGACTTTGGGGAGACTAAGGCAGCCAGAGGGATGCCTTTTTTACAGCTGTTAATCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGtt                           <  1:1382987/149‑1 (MQ=255)
                                                                                                                          aGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCGGTAGACGATAGAACCTTCGACGTTGGTATCGTAGGTTTTCTCGAACAgggt                           >  2:861692/1‑147 (MQ=255)
                                                                                                                            aCGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCAc                                                         >  2:730613/1‑117 (MQ=255)
                                                                                                                              gTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTAttt                       >  2:450309/1‑148 (MQ=255)
                                                                                                                               ttGGGGAGACTAAGGCAGGCAGAGGGCTGGCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGa                      <  1:1343872/149‑1 (MQ=255)
                                                                                                                                tgtggAGTCTAAGGCAGCCTGTTGGCTGCCTTTTTTACAGCTGTTATTAAGAATTGAAACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCATCCACGGTGGTATCGTAGGTTTTGTCTAACATGTTATTGAc                     <  2:429405/146‑1 (MQ=255)
                                                                                                                                   ggAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGtt                  <  2:436869/149‑1 (MQ=255)
                                                                                                                                   ggAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGtt                  >  2:927021/1‑149 (MQ=255)
                                                                                                                                    gagaCTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGAc                     >  2:309794/1‑145 (MQ=255)
                                                                                                                                    gagaCTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGAc                     <  1:309794/145‑1 (MQ=255)
                                                                                                                                      gaCTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTccc               <  2:184758/149‑1 (MQ=255)
                                                                                                                                       aaaaaGGCAGCCATCTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTTATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGATCAGGTTATTGACGTTcccc              >  2:1548799/4‑149 (MQ=255)
                                                                                                                                       aCTAAGGCAGCCAGATGGCTGCCGTATTTACAGGTGTTATTCAGAATTGATACGTGACGGTAATGCTGAAATTATGCGGTGTGCCGTAGACGATAGAACGTTCCACGTTGGTATCGTAGCTTGTGTCTAACAGGGTATTGACGTTcccc              >  2:850909/1‑149 (MQ=255)
                                                                                                                                        cTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCt             <  1:221751/149‑1 (MQ=255)
                                                                                                                                        aaaaGGCAGCCATCTGGCTGCCTTTTTTACAGGTGATATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGCTTATTGCCGTTCCCCt             >  1:460687/3‑149 (MQ=255)
                                                                                                                                         tAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTg            >  1:588427/1‑149 (MQ=255)
                                                                                                                                         tAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGAAGACGATAGAACCGTACACGTTGGTATCGTAGGGTTTGTCGAACAGGTTATTGTCGTTccccgg            >  1:1015456/1‑147 (MQ=255)
                                                                                                                                          aaGGCAGCCATATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATCCGTGACGGTAATGCTGAAATTACGCGGGGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGt           <  2:1341107/149‑1 (MQ=255)
                                                                                                                                           aGGCAGCCATCTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTATGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTa          >  2:47325/1‑149 (MQ=255)
                                                                                                                                           aGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTAAGAATTGATACGTGCCGGTAATGATGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTa          >  1:21775/1‑149 (MQ=255)
                                                                                                                                               aGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCACTGTAACga      >  1:1417020/1‑149 (MQ=255)
                                                                                                                                                   aGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCGGAAATTACGCGGTGTGCCTTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACTAGAAg  <  2:1215041/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAG  >  NC_000913/1159187‑1159482

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: