Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,639,338 T→C 5.1% intergenic (+197/+47) 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,3380TC 5.1% 248.7 / 3.3 78intergenic (+197/+47)pitA/uspBphosphate transporter, low‑affinity; tellurite importer/universal stress (ethanol tolerance) protein B
Reads supporting (aligned to +/- strand):  ref base T (32/42);  new base C (2/2);  total (34/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

GCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAA  >  NC_000913/3639196‑3639485
                                                                                                                                              |                                                                                                                                                   
gcgcTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGttt                                                                                                                                               <  2:306689/149‑1 (MQ=255)
gcgcTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGttt                                                                                                                                               >  1:450128/1‑149 (MQ=255)
  gcTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTAc                                                                                                                                             <  2:360402/149‑1 (MQ=255)
    tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTg                                                                                                                                           <  3:16073/149‑1 (MQ=255)
      tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGc                                                                                                                                         <  2:210242/149‑1 (MQ=255)
        ttCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGCCAGTAAACTGAccc                                                                                                                                       <  4:159884/149‑1 (MQ=255)
         tCATGCCGAATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCg                                                                                                                                      >  4:50911/1‑149 (MQ=255)
          cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGCCAGTAAACTGACCCGc                                                                                                                                     <  2:75729/149‑1 (MQ=255)
             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCttt                                                                                                                                  <  2:28406/149‑1 (MQ=21)
               cGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCttttt                                                                                                                                <  1:338386/149‑1 (MQ=255)
               cGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCttttt                                                                                                                                <  2:550738/149‑1 (MQ=255)
               cGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGAGTCAGTTTACTGGCCCGCttttt                                                                                                                                <  3:424165/149‑1 (MQ=255)
                   tGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTAt                                                                                                                            <  2:176837/149‑1 (MQ=255)
                    gcggcgTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATg                                                                                                                           >  3:102309/1‑149 (MQ=255)
                     cggcgTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGc                                                                                                                          <  4:102309/149‑1 (MQ=255)
                        cgTGAACGCCGTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCATGGATTATGTAGGCCTGATAAGCGTAGCGCATCAGCCAATTTTGCGTTTGTAATCCGTCTCCAGCGGGTCAGTTTACTGGCCCACTTTTTTTATGCAtt                                                                                                                       <  2:243068/149‑1 (MQ=21)
                         gTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGtta                                                                                                                      <  2:64245/149‑1 (MQ=255)
                          tGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGttat                                                                                                                     >  1:58811/1‑149 (MQ=255)
                            aaCGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGttatta                                                                                                                   >  4:329355/1‑149 (MQ=255)
                              cGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTattatt                                                                                                                 <  3:210918/149‑1 (MQ=255)
                               gCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTg                                                                                                                <  4:158688/149‑1 (MQ=255)
                               gCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTg                                                                                                                >  2:186165/1‑149 (MQ=255)
                                 cTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCt                                                                                                              >  1:274371/1‑149 (MQ=255)
                                  tgagccGGCCTACAAAAACATTCAAATTCAATAGATTGCAGAGTTTATGGAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGATCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTg                                                                                                             >  2:30323/5‑149 (MQ=25)
                                     tCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAg                                                                                                          >  2:581771/1‑149 (MQ=255)
                                         gCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCt                                                                                                      <  3:453458/149‑1 (MQ=255)
                                            tACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAGTCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTCaa                                                                                                   <  3:116357/149‑1 (MQ=255)
                                            tACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTaa                                                                                                   >  2:16237/1‑149 (MQ=255)
                                            tACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTaa                                                                                                   <  1:518705/149‑1 (MQ=255)
                                             aCAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCGGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGATTTTTTTATGCGTTATTATTGCTGAAGTGCTTAAt                                                                                                  <  3:30271/149‑1 (MQ=255)
                                              cAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATg                                                                                                 >  2:407479/1‑149 (MQ=255)
                                               aaaaaCATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAAt                                                                                                  <  1:20407/147‑1 (MQ=255)
                                               aaaaaCATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAAt                                                                                                  >  2:20407/1‑147 (MQ=255)
                                                 aaaCATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCa                                                                                              >  3:483870/1‑149 (MQ=255)
                                                        cAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAatcatc                                                                                       <  3:329355/149‑1 (MQ=255)
                                                          aaTTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCaa                                                                                     >  3:389711/1‑149 (MQ=255)
                                                           aTTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAAt                                                                                    <  4:261615/149‑1 (MQ=255)
                                                                 tAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATc                                                                              <  1:27064/149‑1 (MQ=255)
                                                                   gATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAg                                                                            >  1:452701/1‑149 (MQ=255)
                                                                        cAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGa                                                                       <  2:230560/149‑1 (MQ=255)
                                                                         agagATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGac                                                                      >  3:59084/1‑149 (MQ=255)
                                                                         agagATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGac                                                                      >  1:291660/1‑149 (MQ=255)
                                                                           agaTTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGacca                                                                    >  1:398317/1‑149 (MQ=255)
                                                                            gaTTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGaccac                                                                   <  2:274371/149‑1 (MQ=255)
                                                                              ttATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGaccacca                                                                 <  2:398317/149‑1 (MQ=255)
                                                                                 tGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGAccaccacca                                                              <  4:119866/149‑1 (MQ=255)
                                                                                   tAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGa                                                            <  3:158507/149‑1 (MQ=255)
                                                                                    aGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGAccaccacca                                                              >  1:64514/1‑146 (MQ=255)
                                                                                    aGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGAccaccacca                                                              <  2:64514/146‑1 (MQ=255)
                                                                                    aGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGAc                                                           <  3:270108/149‑1 (MQ=255)
                                                                                          gATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTg                                                                                                                >  1:461164/1‑90 (MQ=255)
                                                                                          gATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTg                                                                                                                <  2:461164/90‑1 (MQ=255)
                                                                                              aGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATgc                                                 <  1:30323/149‑1 (MQ=255)
                                                                                                       cATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAg                                        >  1:347801/1‑149 (MQ=255)
                                                                                                        aTCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGa                                       <  1:581771/149‑1 (MQ=255)
                                                                                                           aGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATa                                    >  3:98823/1‑149 (MQ=255)
                                                                                                             gCAATTTTGCGTTTGTAATGAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATaaa                                  >  3:484208/1‑149 (MQ=255)
                                                                                                             gCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATaaa                                  <  3:322400/149‑1 (MQ=255)
                                                                                                             gCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATaaa                                  <  2:291660/149‑1 (MQ=255)
                                                                                                                aTTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTg                               <  3:232975/149‑1 (MQ=255)
                                                                                                                        tttGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCa                       >  4:329556/1‑149 (MQ=255)
                                                                                                                         ttGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCAc                      <  1:16237/149‑1 (MQ=255)
                                                                                                                            tAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATgcgc                                               >  1:415962/1‑121 (MQ=255)
                                                                                                                            tAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATgcgc                                               <  2:415962/121‑1 (MQ=255)
                                                                                                                              aTCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCt                 >  2:274874/1‑149 (MQ=255)
                                                                                                                               tCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTc                >  1:211168/1‑149 (MQ=255)
                                                                                                                                cAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGac                                                                      >  1:44126/1‑94 (MQ=255)
                                                                                                                                cAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGac                                                                      <  2:44126/94‑1 (MQ=255)
                                                                                                                                   tctcGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGt                                                                                                         <  1:566318/56‑1 (MQ=255)
                                                                                                                                   tctcGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGt                                                                                                         >  2:566318/1‑56 (MQ=255)
                                                                                                                                   tctcGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAg            >  2:7128/1‑149 (MQ=255)
                                                                                                                                     tcGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAgcg          >  4:227613/1‑149 (MQ=255)
                                                                                                                                       aaGCGGGCCAGTAAACTGACCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAgcggc        >  3:578922/2‑149 (MQ=255)
                                                                                                                                         gCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGa      <  2:51126/149‑1 (MQ=255)
                                                                                                                                          cGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGaa     <  4:200927/149‑1 (MQ=255)
                                                                                                                                           gggCCAGTAAACTGACCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAAt    >  3:227926/1‑149 (MQ=255)
                                                                                                                                             gTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATaa  <  2:58811/149‑1 (MQ=255)
                                                                                                                                             gTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATaa  >  2:185788/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAA  >  NC_000913/3639196‑3639485

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: