Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,164,154 T→C 7.9% Y172Y (TAT→TAC btuB → vitamin B12/cobalamin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,164,1540TC 7.9% 234.9 / 10.1 76Y172Y (TAT→TACbtuBvitamin B12/cobalamin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (30/40);  new base C (3/3);  total (33/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

GTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAA  >  NC_000913/4164017‑4164293
                                                                                                                                         |                                                                                                                                           
gtTGAATATCTCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGTAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTatga                                                                                                                                           <  2:1559270/140‑1 (MQ=255)
gtTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTatga                                                                                                                                           <  2:302224/140‑1 (MQ=255)
gtTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTatga                                                                                                                                           <  2:302531/140‑1 (MQ=255)
gtTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTatga                                                                                                                                           <  1:388493/140‑1 (MQ=255)
gtTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTatga                                                                                                                                           <  1:164777/140‑1 (MQ=255)
      tatatCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtctc                                                                                                                                     <  1:277996/140‑1 (MQ=255)
       atatCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTctct                                                                                                                                    >  1:982812/1‑140 (MQ=255)
         atCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGATATCAGAACTATGATGTCTCTAc                                                                                                                                  >  1:619055/1‑140 (MQ=255)
          tCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACg                                                                                                                                 <  1:1077150/140‑1 (MQ=255)
          tCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACg                                                                                                                                 <  1:1448875/140‑1 (MQ=255)
           ccGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgc                                                                                                                                <  2:1577064/140‑1 (MQ=255)
           ccGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgc                                                                                                                                >  1:383063/1‑140 (MQ=255)
           ccGGGGGCCGCGCTCCGCTTTTTATGGTACGGATGAAATAGGCGGGGTGGTGAATATCATCTCGACGCGCGATGAACCCGGAACGGAAATTTCTGCGGGGTGGGGACGCATTAGTTAGCAGAACTATGATGTCTCTAGgc                                                                                                                                >  2:52491/1‑140 (MQ=255)
              tGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCTGAACTATGATGTCTCTACgcagc                                                                                                                             <  1:969625/140‑1 (MQ=255)
              tGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcagc                                                                                                                             <  2:1360789/140‑1 (MQ=255)
              tGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcagc                                                                                                                             >  2:1595870/1‑140 (MQ=255)
              tGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcagc                                                                                                                             >  1:1272159/1‑140 (MQ=255)
              tGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcagc                                                                                                                             >  1:38328/1‑140 (MQ=255)
              tGGGCCGCGCGCCGCTGTTTATGGTTCCGATGCAAGAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcagc                                                                                                                             <  2:1360017/140‑1 (MQ=255)
                    gcgcTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCaacaac                                                                                                                       <  1:1450692/140‑1 (MQ=255)
                    gcgcTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCaacaac                                                                                                                       <  1:637469/140‑1 (MQ=255)
                    gcgcTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCaacaac                                                                                                                       <  1:1450699/140‑1 (MQ=255)
                     cgcTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACt                                                                                                                      >  1:521543/1‑140 (MQ=255)
                      gcTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGATTATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTg                                                                                                                     <  2:566608/140‑1 (MQ=255)
                           gCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGa                                                                                                                <  2:351140/140‑1 (MQ=255)
                            cTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGAt                                                                                                               <  1:873313/140‑1 (MQ=255)
                            cTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGAt                                                                                                               <  1:614230/140‑1 (MQ=255)
                                 tATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGac                                                                                                          <  1:5895/140‑1 (MQ=255)
                                 tATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGac                                                                                                          >  2:1048634/1‑140 (MQ=255)
                                 tATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTATACGCAGCAACAACTGGGGGATAAGac                                                                                                          <  1:1361261/140‑1 (MQ=255)
                                         cGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTaa                                                                                                  >  1:883968/1‑140 (MQ=255)
                                           aTGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACg                                                                                                <  1:515131/140‑1 (MQ=255)
                                           aTGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACg                                                                                                >  1:166998/1‑140 (MQ=255)
                                           aTGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACg                                                                                                >  1:26146/1‑140 (MQ=255)
                                           aTGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACg                                                                                                <  2:1517955/140‑1 (MQ=255)
                                                   ggCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTggg                                                                                        <  1:863816/140‑1 (MQ=255)
                                                   ggCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTggg                                                                                        <  1:1029580/140‑1 (MQ=255)
                                                   ggCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTggg                                                                                        <  1:1393143/140‑1 (MQ=255)
                                                   ggCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTggg                                                                                        <  1:1176514/140‑1 (MQ=255)
                                                   ggCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTggg                                                                                        <  2:521543/140‑1 (MQ=255)
                                                       gggTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGAt                                                                                    <  1:432945/140‑1 (MQ=255)
                                                       gggTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGAt                                                                                    <  1:1591784/140‑1 (MQ=255)
                                                        ggtggtGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGAtt                                                                                   >  2:741087/1‑140 (MQ=255)
                                                             tgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGcc                                                                              <  1:1501333/140‑1 (MQ=255)
                                                             tgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGcc                                                                              <  1:45912/140‑1 (MQ=255)
                                                                     atcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCa                                                                      <  2:588145/140‑1 (MQ=255)
                                                                      tcaCGTCGCGCGATGAACCCGGAACGGTAATTTCAGCAGGGTGGGGAAGCAATAGTTTTCAGAACTATGATGTCTCTTCGCAGCAACATCTGGTGGATAAGACACGGGTAACGCTGTTGGGCGTTTTTGCCCTTACTCag                                                                     >  2:913039/1‑139 (MQ=255)
                                                                      tcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCAt                                                                     <  1:1507284/140‑1 (MQ=255)
                                                                      tcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCAt                                                                     >  2:886725/1‑140 (MQ=255)
                                                                      tcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCAt                                                                     <  1:76037/140‑1 (MQ=255)
                                                                      tcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCAt                                                                     >  2:886624/1‑140 (MQ=255)
                                                                      tcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCAt                                                                     >  2:256093/1‑140 (MQ=255)
                                                                      tcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCAt                                                                     <  2:38328/140‑1 (MQ=255)
                                                                       caCGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATg                                                                    >  1:574430/1‑140 (MQ=255)
                                                                         cgacgCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGt                                                                  >  1:820265/1‑140 (MQ=255)
                                                                            cgcgcgATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTat                                                               >  1:566013/1‑140 (MQ=255)
                                                                             gcgcgATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTatg                                                              <  2:1420681/140‑1 (MQ=255)
                                                                             gcgcgATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTatg                                                              <  2:1420819/140‑1 (MQ=255)
                                                                                   tGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGAtgttg                                                        <  2:820265/140‑1 (MQ=255)
                                                                                        ccGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGcc                                                   <  1:1536482/140‑1 (MQ=255)
                                                                                        ccGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGcc                                                   <  2:852413/140‑1 (MQ=255)
                                                                                        ccGGAACGGAAATTACAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGcc                                                   <  2:1270603/140‑1 (MQ=255)
                                                                                           gaacggaaATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTACGATGTCTCTACGCAGCAACAACTGGGGGATaa                                                                                                             <  1:451878/79‑1 (MQ=255)
                                                                                           gaacggaaATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTACGATGTCTCTACGCAGCAACAACTGGGGGATaa                                                                                                             >  2:450883/1‑79 (MQ=255)
                                                                                           gaacggaaATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTACGATGTCTCTACGCAGCAACAACTGGGGGATaa                                                                                                             >  2:450894/1‑79 (MQ=255)
                                                                                           gaacggaaATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTACGATGTCTCTACGCAGCAACAACTGGGGGATaa                                                                                                             >  2:451878/1‑79 (MQ=255)
                                                                                           gaacggaaATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTACGATGTCTCTACGCAGCAACAACTGGGGGATaa                                                                                                             <  1:450894/79‑1 (MQ=255)
                                                                                           gaacggaaATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTACGATGTCTCTACGCAGCAACAACTGGGGGATaa                                                                                                             <  1:450883/79‑1 (MQ=255)
                                                                                                 aaaTTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTTCGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATTGTAAt                                          >  2:323743/1‑140 (MQ=255)
                                                                                                         gcagGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAAc                                  >  2:551369/1‑140 (MQ=255)
                                                                                                         gcagGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAAc                                  >  2:551400/1‑140 (MQ=255)
                                                                                                         gcagGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAAc                                  >  2:562563/1‑140 (MQ=255)
                                                                                                              gTGGGGAAGCAATAGCTCTCAGAACCATGAGTTCTCTACGCAGCATCACCTGGGGGATATCACACGGGTAACGCTGATGGGAGATTATGCCCATACTCATGGTTATGATGATGTTGGCTATTTGAATACCGGAACGCAAg                             <  2:875363/140‑1 (MQ=255)
                                                                                                                      gCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGcagaca                     <  1:1111511/140‑1 (MQ=255)
                                                                                                                       cAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGcagacag                    >  1:110803/1‑140 (MQ=255)
                                                                                                                              tATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGa             >  1:922944/1‑140 (MQ=255)
                                                                                                                               aTCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGAt                                                                                    <  2:559182/68‑1 (MQ=255)
                                                                                                                               aTCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGAt                                                                                    >  1:559182/1‑68 (MQ=255)
                                                                                                                                      cTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACTGGTAACGCTGTTGGGCGTTTATGCCCATACTCATGTTTATGATGTTGTAGCCTTTGGTAATAACGTAAATCAAGCTCACACAGATAAAGAAGGttttt     >  2:1673858/1‑140 (MQ=255)
                                                                                                                                      cTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGAGGGttttt     >  2:478540/1‑140 (MQ=255)
                                                                                                                                       tATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGtttttt    >  2:167878/1‑140 (MQ=255)
                                                                                                                                         tgatgTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTaa  >  2:114920/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                           
GTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAA  >  NC_000913/4164017‑4164293

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: