Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,159,336 A→C 5.6% intergenic (+34/+26) 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,3360AC 5.6% 266.1 / 5.9 107intergenic (+34/+26)ptsG/fhuEfused glucose‑specific PTS enzymes: IIB component/IIC component/ferric‑rhodotorulic acid outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (54/46);  new base C (3/3);  total (57/50)
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.12e-01

Only 100 of 108 total aligned reads displayed.

GAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAG  >  NC_000913/1159188‑1159482
                                                                                                                                                    |                                                                                                                                                  
gaagaaACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGa                                                                                                                                                    >  1:184758/1‑149 (MQ=255)
 aagaaACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGAt                                                                                                                                                   <  1:255833/149‑1 (MQ=255)
 aagaaACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGAt                                                                                                                                                   >  1:1215041/1‑149 (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)
                          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)
                                                                     cGATAACCTGAAAACCGAGATGGAGGAGTACATCCGTAACCACTAATCCGCAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGc                                                                               <  2:567934/149‑1 (MQ=255)
                                                                      gATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGcc                                                                              >  1:1458952/1‑149 (MQ=255)
                                                                      gATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGcc                                                                              <  2:1220868/149‑1 (MQ=255)
                                                                      gATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTATTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGcc                                                                              >  1:575811/1‑149 (MQ=255)
                                                                           ccTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCTAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTATGCGGTGTGCCGTAgc                                                                         >  2:195885/1‑148 (MQ=255)
                                                                           ccTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGa                                                                         <  1:502055/149‑1 (MQ=255)
                                                                           ccTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGa                                                                         >  2:1518013/1‑149 (MQ=255)
                                                                            cTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGAc                                                                        >  2:423286/1‑149 (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)
                                                                                         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)
                                                                                               aGTACATCCGTAACAACTAATACGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCTGAATTGATACGTGCCGGTAATGATGAAATTACGCGGTGTGCCGTAGACGATATAACCTTCAAAGTTg                                                     >  2:603345/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)
                                                                                                         tAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCAGGTAATGCTGAAATTACGCGCTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAgg                                           >  1:1279840/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)
                                                                                                             cACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGtttt                                       >  1:442822/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)
                                                                                                                  aTCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGa                                  >  2:1137516/1‑149 (MQ=255)
                                                                                                                     cGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACa                               <  1:1369033/149‑1 (MQ=255)
                                                                                                                     cGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTAATTCAGAATTGATACGTTGCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACTTTAGAACCTTCCACGTTTGTATCGTAGGTTTTGTCGATCa                               <  2:1123879/149‑1 (MQ=255)
                                                                                                                         aGACTTTGGGGAGACTAAGGCAGCCAGAGGGATGCCTTTTTTACAGCTGTTAATCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGtt                           <  1:1382987/149‑1 (MQ=255)
                                                                                                                         aGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCGGTAGACGATAGAACCTTCGACGTTGGTATCGTAGGTTTTCTCGAACAgggt                           >  2:861692/1‑147 (MQ=255)
                                                                                                                           aCGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCAc                                                         <  1:730613/117‑1 (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                     <  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)
                                                                                                                                       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)
                                                                                                                                                    |                                                                                                                                                  
GAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAG  >  NC_000913/1159188‑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: