Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,639,165 T→C 5.5% intergenic (+24/+220) 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,1650TC 5.5% 301.5 / 7.2 113intergenic (+24/+220)pitA/uspBphosphate transporter, low‑affinity; tellurite importer/universal stress (ethanol tolerance) protein B
Reads supporting (aligned to +/- strand):  ref base T (52/51);  new base C (4/2);  total (57/56)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.81e-01

Only 100 of 113 total aligned reads displayed.

GTGGCGGCTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTG  >  NC_000913/3639026‑3639313
                                                                                                                                           |                                                                                                                                                    
gtggCGGCTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGccc                                                                                                                                            >  2:1510957/1‑150 (MQ=255)
gtggCGGCTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGccc                                                                                                                                            >  2:1467107/1‑150 (MQ=255)
 tggCGGCTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCg                                                                                                                                           >  2:1027806/1‑150 (MQ=255)
  ggcggcTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGc                                                                                                                                          <  1:1323024/150‑1 (MQ=255)
    cggcTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGCCAGCTGACCCGCtc                                                                                                                                        <  1:656321/150‑2 (MQ=255)
    cggcTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGACAGCTGGCCCGCtt                                                                                                                                        >  1:184660/1‑150 (MQ=255)
      gcTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCa                                                                                                                                      <  1:897285/150‑1 (MQ=255)
      gcTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCa                                                                                                                                      <  1:936581/150‑1 (MQ=255)
      gcTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCa                                                                                                                                      >  1:992245/1‑150 (MQ=255)
      gcTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGCCAGCTGACCCGCtcgt                                                                                                                                      <  1:1608284/150‑4 (MQ=255)
         tACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGAt                                                                                                                                   <  2:1377689/150‑1 (MQ=255)
           cAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATtg                                                                                                                                 >  1:621382/1‑150 (MQ=255)
             gCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATtgtg                                                                                                                               >  1:1129816/1‑150 (MQ=255)
               gTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGAc                                                                                                                             >  2:734525/1‑150 (MQ=255)
               gTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGAc                                                                                                                             >  2:412377/1‑150 (MQ=255)
               gTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGAc                                                                                                                             <  2:1679671/150‑1 (MQ=255)
                tAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACa                                                                                                                            <  2:163373/150‑1 (MQ=255)
                 aaaaCCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACAt                                                                                                                           <  1:1027806/150‑1 (MQ=255)
                    aCCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGt                                                                                                                        >  1:1623838/1‑150 (MQ=255)
                          aCCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCtt                                                                                                                  <  2:1021362/150‑1 (MQ=255)
                             aGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGt                                                                                                               >  2:928364/1‑150 (MQ=255)
                              gCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCGGCGGTACTGCTTGCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGAAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGtt                                                                                                              <  1:1588197/150‑1 (MQ=255)
                                aTTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCa                                                                                                            <  2:495438/150‑1 (MQ=255)
                                 ttcgtaTGGGCTGGGTGTTTACCCTTCCGGCTGCGGGACTGCTTTCCGGCGGGCTGTAATGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCAt                                                                                                           <  2:1471111/145‑1 (MQ=255)
                                  tCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATg                                                                                                          <  2:621382/150‑1 (MQ=255)
                                   cTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGc                                                                                                         <  1:1636448/150‑1 (MQ=255)
                                   cTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGc                                                                                                         >  1:1140772/1‑150 (MQ=255)
                                            tgagtgtTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATgcgg                                                                                                <  2:962280/147‑1 (MQ=255)
                                             gggTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATgcggc                                                                                               >  1:324646/1‑150 (MQ=255)
                                             gggTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATgcggc                                                                                               <  1:1561020/150‑1 (MQ=255)
                                                 gtTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGa                                                                                           >  2:1315876/1‑150 (MQ=255)
                                                   ttACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAAc                                                                                         <  2:894458/150‑1 (MQ=255)
                                                   ttACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAAc                                                                                         <  2:1699360/150‑1 (MQ=255)
                                                    tACCCTTCCGGCTGCGGTACTGCTTTCCAGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACg                                                                                        >  1:616769/1‑150 (MQ=255)
                                                      cccTTGCGGCTGCGGTACTGCGTTCCGGCGGGCTGTACTGGCTCGCCTTGCAGTTCCGGTAATCGTACGCACCTAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACTTAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                                      <  1:140653/150‑1 (MQ=39)
                                                      cccTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGCGCGGGTCAGCGGGCCCGCTTCAGATTGTGAAATAGTGCGGTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                                      >  1:1189953/1‑150 (MQ=255)
                                                       ccTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                                     >  1:1458212/1‑150 (MQ=255)
                                                       ccTTCCGGCGGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                                     <  2:1036099/150‑1 (MQ=255)
                                                          tCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt                                                                                  <  2:405817/150‑1 (MQ=255)
                                                             ggCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                               >  2:114307/1‑150 (MQ=255)
                                                             ggCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCg                                                                               <  2:973632/150‑1 (MQ=255)
                                                                ggCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATTGTACGCACCAAAACGCGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGGGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                            <  1:1036836/149‑1 (MQ=255)
                                                                 gCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTATTCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCt                                                                           <  2:1600706/150‑1 (MQ=255)
                                                                 gCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCACAACCAGCGGGGCAGCGCGCCCGCTTCAGATTGTGACAAAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCt                                                                           <  1:858382/150‑1 (MQ=38)
                                                                 gCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCt                                                                           <  2:787296/150‑1 (MQ=255)
                                                                  cGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                          >  1:1130163/1‑150 (MQ=255)
                                                                      aCTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaa                                                                      >  2:130521/1‑150 (MQ=255)
                                                                      aCTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaa                                                                      >  1:955196/1‑150 (MQ=255)
                                                                        tGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaaaa                                                                    >  1:165180/1‑150 (MQ=255)
                                                                         gCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAAc                                                                   >  1:1705832/1‑150 (MQ=255)
                                                                         gCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAAc                                                                   >  2:1225730/1‑150 (MQ=255)
                                                                          cTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACa                                                                  >  1:1126454/1‑150 (MQ=255)
                                                                          cTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACa                                                                  <  1:1467107/150‑1 (MQ=255)
                                                                          cTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACa                                                                  <  2:1129816/150‑1 (MQ=255)
                                                                          cTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACa                                                                  >  1:54260/1‑150 (MQ=255)
                                                                          cTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACa                                                                  <  1:734525/150‑1 (MQ=255)
                                                                            ttCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACAtt                                                                >  2:523527/1‑149 (MQ=255)
                                                                             tCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGc                                                               >  2:401800/1‑150 (MQ=255)
                                                                              ccGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGAGAACGCCTTATCCGGCCTACAAAAACATGCa                                                              >  2:1635984/1‑150 (MQ=255)
                                                                               gggcggGCTGTACTGGCTCTCCTTGCAGATCCTGTACTCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCaa                                                             <  2:184660/149‑1 (MQ=255)
                                                                                gtcggGCTGTACTGGCTCACCGTGCAGTTCCTGCAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCaaa                                                            <  1:1262107/148‑1 (MQ=255)
                                                                                    ggCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCa                                                        >  2:295829/1‑150 (MQ=255)
                                                                                       tGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGACGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATa                                                     >  2:983056/1‑150 (MQ=255)
                                                                                        gTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAg                                                    <  1:1428740/150‑1 (MQ=255)
                                                                                         tACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGa                                                   >  1:489993/1‑150 (MQ=255)
                                                                                         tACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGa                                                   >  1:1207937/1‑150 (MQ=255)
                                                                                          aCTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGAt                                                  >  2:468193/1‑150 (MQ=255)
                                                                                             ggCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTga                                               >  2:147697/1‑149 (MQ=255)
                                                                                             ggCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGc                                               <  2:1179376/150‑1 (MQ=255)
                                                                                              gCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCa                                              >  1:978314/1‑150 (MQ=255)
                                                                                                tctcCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCaga                                            <  2:733437/150‑1 (MQ=255)
                                                                                                tctcCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCaga                                            >  2:806994/1‑150 (MQ=255)
                                                                                                        cAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGt                                    >  2:782936/1‑150 (MQ=255)
                                                                                                           gtCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGTTGGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAgg                                 <  2:152431/149‑1 (MQ=255)
                                                                                                             ccTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGcc                               >  1:129676/1‑150 (MQ=255)
                                                                                                             ccTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGGGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGcc                               <  2:322606/150‑1 (MQ=255)
                                                                                                              cTGTTATCGTCCGCACCAACACTAGCGGGTCAGCGGGCCTGCTTCAGATTGTGACATAGTGAGCTTTGTTCCTGCATGATTCGGCGTGATATAATTATCCGGACTACAAAAACATTCAAATTCAATAGATTGCAGAGATTATGTAGGCCt                              <  1:71687/150‑1 (MQ=255)
                                                                                                              cTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCt                              <  1:412377/150‑1 (MQ=255)
                                                                                                                 tAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGAt                           >  2:210880/1‑150 (MQ=255)
                                                                                                                  aaTCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATa                          <  1:1510957/150‑1 (MQ=255)
                                                                                                                  aaTCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATa                          >  2:1456301/1‑150 (MQ=255)
                                                                                                                      gTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCg                      >  1:422190/1‑150 (MQ=255)
                                                                                                                        aCGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTa                    >  1:880840/1‑150 (MQ=255)
                                                                                                                          gCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAgc                  <  2:1701322/150‑1 (MQ=255)
                                                                                                                            aCCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAgcgc                >  2:795536/1‑150 (MQ=255)
                                                                                                                             ccAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCa               >  2:915678/1‑150 (MQ=255)
                                                                                                                              aatctgaaGCGGGCCAGCTGACCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCAt              >  1:1331270/8‑150 (MQ=255)
                                                                                                                               atctgaaGCGGGCCAGCTGACCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATc             >  2:1523579/7‑150 (MQ=255)
                                                                                                                               aaaaCGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATc             <  2:129676/150‑1 (MQ=255)
                                                                                                                               aaaaCGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATc             >  1:652751/1‑150 (MQ=255)
                                                                                                                                  aCGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAgg          <  2:989493/150‑1 (MQ=255)
                                                                                                                                  aCGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAgg          <  2:1126454/150‑1 (MQ=255)
                                                                                                                                   gaaGCGGGCCAGCTGACCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGc         >  1:557526/3‑150 (MQ=255)
                                                                                                                                   gaaGCGGGCCAGCTGACCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGc         >  2:1406453/3‑150 (MQ=255)
                                                                                                                                   cgcgCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGc         <  1:334578/147‑1 (MQ=255)
                                                                                                                                     aGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCaa       >  2:801028/1‑150 (MQ=255)
                                                                                                                                      gCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAAt      <  1:637540/150‑1 (MQ=255)
                                                                                                                                      gCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAAt      <  1:114307/150‑1 (MQ=255)
                                                                                                                                      gCGGGTCAGCTGGCCCGCTGCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAAt      <  1:32696/150‑1 (MQ=255)
                                                                                                                                          gTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTg  <  2:1107796/150‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
GTGGCGGCTTACAGCGTAAAACCGTGACCAGCATTCTGATGGCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTG  >  NC_000913/3639026‑3639313

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: