Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,639,343 T→A 5.4% intergenic (+202/+42) pitA → / ← uspB phosphate transporter, low‑affinity; tellurite importer/universal stress (ethanol tolerance) protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,639,3430TA 5.4% 372.4 / 8.8 129intergenic (+202/+42)pitA/uspBphosphate transporter, low‑affinity; tellurite importer/universal stress (ethanol tolerance) protein B
Reads supporting (aligned to +/- strand):  ref base T (57/65);  new base A (2/5);  total (59/70)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.36e-01

Only 100 of 134 total aligned reads displayed.

GTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATC  >  NC_000913/3639194‑3639492
                                                                                                                                                     |                                                                                                                                                     
gTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGtt                                                                                                                                                       >  2:1142680/1‑150 (MQ=255)
  gcgcTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATACGGCCTACAAAAACATGCAAATTAAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATAAGGCAATTTTGCGTTTGTAATCAGACTCGAGCGGGTCAGTTTa                                                                                                                                                     >  2:1212117/1‑150 (MQ=255)
    gcTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACt                                                                                                                                                   <  2:415614/150‑1 (MQ=255)
    gcTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACt                                                                                                                                                   <  1:619518/150‑1 (MQ=255)
     cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTg                                                                                                                                                  <  2:438520/150‑1 (MQ=255)
          ttCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCTAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCg                                                                                                                                             <  2:49985/150‑1 (MQ=255)
           tCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGCCAGTAAACTGACCCGc                                                                                                                                            <  1:749488/150‑1 (MQ=255)
            cATGCCGGCTGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCATTAGATTGCAGAGATTATGTAGGACTGATAAGCGTAGAGCATCAGGCAATTTTGCGGTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCt                                                                                                                                           <  2:151414/150‑1 (MQ=21)
            cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGCCAGTAAACTGACCCGCt                                                                                                                                           <  1:184021/150‑1 (MQ=255)
            cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGCCAGTAAACTGACCCGCt                                                                                                                                           <  1:824019/150‑1 (MQ=255)
               gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCCGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCtttt                                                                                                                                        >  1:311280/1‑150 (MQ=21)
                ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGCCAGTAAACTGACCCGCtcgag                                                                                                                                       <  1:1013277/150‑5 (MQ=18)
                ccGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCttttt                                                                                                                                       >  2:657630/1‑150 (MQ=25)
                  ggataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGAATATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCTGGTCAGTTTACTGGCCCGCttttttt                                                                                                                                     >  2:1379755/7‑150 (MQ=21)
                  ggATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGCGCGGGTCAGTTTACTGGCCCGCttttttt                                                                                                                                     >  1:354865/1‑150 (MQ=255)
                    aTGCGGCTGGAACGCCTTATCCGGCCTACAAAAAACTGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGAAAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTAt                                                                                                                                   <  1:77276/150‑1 (MQ=25)
                     tGCGGCGTGAACGTCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTTCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATg                                                                                                                                  <  2:1622139/150‑1 (MQ=255)
                     tGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATg                                                                                                                                  <  1:1324215/150‑1 (MQ=255)
                     tGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATg                                                                                                                                  <  1:398919/150‑1 (MQ=255)
                     tGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATg                                                                                                                                  <  2:1609067/150‑1 (MQ=255)
                      tcggcgTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGc                                                                                                                                 <  1:1096510/149‑1 (MQ=255)
                      gcggcgTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAAGCAATTTTGCGTTTGTAATCAGTCTCGAGCCGGTCAGTTTACTGGCCCGCTTTTTTTATGc                                                                                                                                 >  2:684746/1‑150 (MQ=255)
                        ggcgTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGt                                                                                                                               >  1:1378332/1‑150 (MQ=255)
                        ggcgTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGAATGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGt                                                                                                                               >  1:1384989/1‑150 (MQ=255)
                           gTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGttat                                                                                                                            <  2:1399705/150‑1 (MQ=255)
                              aaCGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTtattat                                                                                                                         >  2:504014/1‑150 (MQ=255)
                               aCGCCTTATCCGGCATTAAAAAACATGCAAATTCAATAGTTTGCATAGATTATGTAGGCCTGAGAAGCGTAGGGCATGAGGCAATTTTTCGTTTTTAATCAGTATCGAGCGGGTCTGTTTACTGGCTAGCTTTTTTTATGCGTTattatt                                                                                                                        <  2:443414/150‑1 (MQ=255)
                                cGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTg                                                                                                                       <  1:988968/150‑1 (MQ=255)
                                 gCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGc                                                                                                                      >  2:414876/1‑150 (MQ=255)
                                  ccTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCt                                                                                                                     >  1:1482088/1‑150 (MQ=255)
                                  ccTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCt                                                                                                                     >  2:415050/1‑150 (MQ=255)
                                  ccTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCt                                                                                                                     >  1:750537/1‑150 (MQ=255)
                                   cTTATCCGGCCTACAAAAACATGCAAATTCTATAGATTGGAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTTCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTg                                                                                                                    <  1:718732/150‑1 (MQ=255)
                                   cTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTg                                                                                                                    >  2:1309119/1‑150 (MQ=255)
                                   cTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTg                                                                                                                    <  1:1731012/150‑1 (MQ=255)
                                       tCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGt                                                                                                                <  1:319091/150‑1 (MQ=255)
                                             cTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTaa                                                                                                          >  2:17190/1‑150 (MQ=255)
                                             cTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTaa                                                                                                          >  2:1407032/1‑150 (MQ=255)
                                              tACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAAt                                                                                                         <  1:1150431/150‑1 (MQ=255)
                                              tACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATg                                                                                                        >  1:211507/1‑150 (MQ=255)
                                               aCAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATg                                                                                                        <  1:1270302/150‑1 (MQ=255)
                                                  aaaaCATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCa                                                                                                     >  2:1067283/1‑150 (MQ=255)
                                                  aaaaCATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCCCATCAGGCAATTTAGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATCAGTTATTATTGCTGAAGTGCTTAATGCCa                                                                                                     >  1:807413/1‑150 (MQ=255)
                                                      cATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAatc                                                                                                 >  2:1447052/1‑150 (MQ=255)
                                                         gCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAatcata                                                                                              >  1:1505096/1‑149 (MQ=255)
                                                         gCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGGCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAatcatc                                                                                              <  2:857965/150‑1 (MQ=255)
                                                          cAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAAtcatca                                                                                             <  1:1171147/150‑1 (MQ=255)
                                                            aaTTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAAt                                                                                           >  2:228202/1‑150 (MQ=255)
                                                             aTTCAATAGATTGCAGAGATTATGTAGGGCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGGCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATg                                                                                          <  1:1066697/150‑1 (MQ=255)
                                                             aTTCAATAGATTGCAGAGATTATGTAGGCCTGGTAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATg                                                                                          <  2:1398901/150‑1 (MQ=255)
                                                             aTTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATg                                                                                          <  1:1160142/150‑1 (MQ=255)
                                                                cAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATAATTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAAAGCCAAATCATCAATGaaa                                                                                       >  2:778834/1‑150 (MQ=255)
                                                                   tAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCa                                                                                    <  2:105903/150‑1 (MQ=255)
                                                                    aGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAg                                                                                   <  1:838550/150‑1 (MQ=255)
                                                                    aGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAg                                                                                   <  2:660876/150‑1 (MQ=255)
                                                                     gATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAgg                                                                                  >  2:1106361/1‑150 (MQ=255)
                                                                     gATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGGTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAgg                                                                                  >  1:682022/1‑150 (MQ=255)
                                                                       ttGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCGCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCt                                                                                <  2:675395/150‑1 (MQ=255)
                                                                         gCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGa                                                                              >  2:1251316/1‑150 (MQ=255)
                                                                           agagATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGacc                                                                            >  2:1692749/1‑150 (MQ=255)
                                                                             agaTTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGaccac                                                                          <  2:802639/150‑1 (MQ=255)
                                                                             agaTTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGaccac                                                                          >  2:299898/1‑150 (MQ=255)
                                                                                ttATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGaccaccac                                                                       <  2:1378332/150‑1 (MQ=255)
                                                                                  aTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTTCTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGAccaccacca                                                                     >  1:1373217/1‑150 (MQ=255)
                                                                                     tAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGAc                                                                  <  1:1142680/150‑1 (MQ=255)
                                                                                      aGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGAcc                                                                 <  1:642214/150‑1 (MQ=255)
                                                                                       ggCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACca                                                                >  2:885852/1‑150 (MQ=255)
                                                                                        gCCTGATAAGCGTAGCGCAACAGGCAATTTTGCGTTTGTAAACAGTCTCGAGCGGGTCACTTTACTGGCCCGCTTTTATTATGCGTTATTATTGCTGAACTCCTAAATGCCAAATCATCAATGCAATCAGTCTGACCCCCACCAGAcaaa                                                               >  2:1295711/1‑147 (MQ=38)
                                                                                              tAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATg                                                         >  1:1570330/1‑150 (MQ=255)
                                                                                               aaGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATgc                                                        <  2:1612651/150‑1 (MQ=255)
                                                                                                aGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATgcg                                                       >  2:30312/1‑150 (MQ=255)
                                                                                                    tAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTgg                                                   <  2:850326/150‑1 (MQ=255)
                                                                                                    tAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTgg                                                   >  2:105896/1‑150 (MQ=255)
                                                                                                     aGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGTGCTGGt                                                  <  1:1155522/150‑1 (MQ=255)
                                                                                                      gcgcATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTc                                                 <  1:400898/150‑1 (MQ=255)
                                                                                                         cATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGa                                              <  1:1251316/150‑1 (MQ=255)
                                                                                                         cATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGa                                              <  1:1379672/150‑1 (MQ=255)
                                                                                                         cATCAGGCAATTTTGCGTTTGTAATCAGTATCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGa                                              >  2:820186/1‑150 (MQ=255)
                                                                                                          aTCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGaa                                             <  2:829055/150‑1 (MQ=255)
                                                                                                           tCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATTCGCTGGTCAGAAt                                            <  1:871058/150‑1 (MQ=255)
                                                                                                               gCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAAc                                        <  1:975398/150‑1 (MQ=255)
                                                                                                                cAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACt                                       <  1:1211769/150‑1 (MQ=255)
                                                                                                                  aTTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGc                                     <  1:530620/150‑1 (MQ=255)
                                                                                                                      tGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGAc                                 <  1:692516/150‑1 (MQ=255)
                                                                                                                           ttGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCAcc                            <  2:53258/150‑1 (MQ=255)
                                                                                                                               aaTCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCt                        <  1:89844/150‑1 (MQ=255)
                                                                                                                                     tctcGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAgc                  <  1:1581154/150‑1 (MQ=255)
                                                                                                                                        cGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAgcggc               >  1:155152/1‑150 (MQ=255)
                                                                                                                                         gcccGGGCCCGTTTACTGCCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATTATCAATGCAATCAGGCTGACCACCACCATACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAgcggcg              <  2:295937/147‑1 (MQ=255)
                                                                                                                                          aGCGGGCCAGTAAACTGACCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGTAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGa             >  1:994571/1‑150 (MQ=255)
                                                                                                                                          aGCGGGCCAGTAAACTGACCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGa             >  1:602069/1‑150 (MQ=255)
                                                                                                                                           gCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAAATGCCGACGCACCCGCTCACAGCGGCGaa            >  1:1217546/1‑150 (MQ=255)
                                                                                                                                           gCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAAAGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGaa            >  1:1278003/1‑150 (MQ=255)
                                                                                                                                              ggTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATaa         <  1:437400/150‑1 (MQ=255)
                                                                                                                                                tCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAAt       <  2:1482088/150‑1 (MQ=255)
                                                                                                                                                 cAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAAtt      <  1:467412/150‑1 (MQ=255)
                                                                                                                                                   ttttACTGGCCCGATTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATTTTTTATTAATTAATGCTGTACAACAGCTTTCTACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATtca    <  2:1314050/149‑1 (MQ=255)
                                                                                                                                                   gTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATtca    >  2:439534/1‑150 (MQ=255)
                                                                                                                                                   gTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATtca    <  1:630619/150‑1 (MQ=255)
                                                                                                                                                     ttACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATtcatc  <  1:1485613/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATC  >  NC_000913/3639194‑3639492

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: