Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,117,430 (C)5→6 coding (975/1281 nt) wcaK ← putative colanic acid biosynthesis pyruvyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,117,4251.C98.6% 255.0 / ‑3.5 71coding (980/1281 nt)wcaKputative colanic acid biosynthesis pyruvyl transferase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (31/39);  total (31/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

ATTGCCATCTCCGGTAGTCCCAGCTGCTGCATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAGGATCGCTGATGTGCTGGCGCAGGTTGAGCGCCACCATGCGGTCGTCTTTGTTATAGCTGTCAATGCCCGTA  >  NC_000913/2117290‑2117559
                                                                                                                                        |                                                                                                                                       
aTTGCCATCTCCGGTAGTCCCAGCTGCTGCATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGcccc                                                                                                                                    >  2:993612/1‑140 (MQ=255)
 ttGCCATCTCCGGTAGTCCCAGCTGCTGCATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGccccc                                                                                                                                   <  1:602511/140‑1 (MQ=255)
    ccATCTCCGGTAGTCCCAGCTGCTGCATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGccccccag                                                                                                                                >  1:689403/1‑137 (MQ=255)
    ccATCTCCGGTAGTCCCAGCTGCTGCATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGccccccag                                                                                                                                >  2:1066270/1‑137 (MQ=255)
        ctcCGGTAGTCCCAGCTGCTGCATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAAtt                                                                                                                             >  1:149927/1‑140 (MQ=255)
             gTAGTCCCAGCTGCTGCATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGcc                                                                                                                        <  1:726770/140‑1 (MQ=255)
             gTAGTCCCAGCTGCTGCATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGcc                                                                                                                        <  2:689403/140‑1 (MQ=255)
             gTAGTCCCAGCTGCTGCATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGcc                                                                                                                        <  1:368707/140‑1 (MQ=255)
             gTAGTCCCAGCTGCTGCATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGcc                                                                                                                        <  2:444696/140‑1 (MQ=255)
             gTAGTCCCAGCTGCTGCATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGcc                                                                                                                        >  1:666824/1‑140 (MQ=255)
             gTAGTCCCAGCTGCTGCATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGcc                                                                                                                        <  2:318103/140‑1 (MQ=255)
              tAGTCCCAGCTGCTGCATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGccc                                                                                                                       >  2:434589/1‑140 (MQ=255)
              tAGTCCCAGCTGCTGCATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGccc                                                                                                                       <  2:17444/140‑1 (MQ=255)
              tAGTCCCAGCTGCTGCATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGccc                                                                                                                       <  2:261638/140‑1 (MQ=255)
              tAGTCCCAGCTGCTGCATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGccc                                                                                                                       <  1:1066270/140‑1 (MQ=255)
                 tCCCAGCTGCTGCATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCAtt                                                                                                                    >  2:438246/1‑140 (MQ=255)
                   ccAGCTGCTGCATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTc                                                                                                                  >  2:532582/1‑140 (MQ=255)
                   ccAGCTGCTGCATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTc                                                                                                                  >  2:540922/1‑140 (MQ=255)
                   ccAGCTGCTGCATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTc                                                                                                                  >  2:590882/1‑140 (MQ=255)
                          ctgcATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGccc                                                                                                                       <  2:946433/128‑1 (MQ=255)
                          ctgcATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGccc                                                                                                                       >  1:946433/1‑128 (MQ=255)
                          ctgcATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCg                                                                                                           >  1:770246/1‑140 (MQ=255)
                          ctgcATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCg                                                                                                           <  1:300713/140‑1 (MQ=255)
                           tgcATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGt                                                                                                          >  1:100651/1‑140 (MQ=255)
                            gcATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGtt                                                                                                         >  1:200047/1‑140 (MQ=255)
                                aaTCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGaat                                                                                                     <  2:100651/140‑3 (MQ=255)
                                aaTCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGt                                                                                                     <  1:540922/140‑1 (MQ=255)
                                  tCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGt                                                                                                          >  1:804004/1‑133 (MQ=255)
                                  tCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGt                                                                                                          <  2:804004/133‑1 (MQ=255)
                                  tCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTc                                                                                                   >  2:13018/1‑140 (MQ=255)
                                   cccGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTca                                                                                                  <  2:666824/140‑1 (MQ=255)
                                   cccGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTca                                                                                                  <  1:438246/140‑1 (MQ=255)
                                    ccGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTcat                                                                                                 <  2:921270/140‑1 (MQ=255)
                                    ccGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTcat                                                                                                 <  2:200047/140‑1 (MQ=255)
                                     cggcggATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTcatc                                                                                                <  2:11853/140‑1 (MQ=255)
                                       gcggATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTcatcca                                                                                              >  2:4252/1‑140 (MQ=255)
                                       gcggATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTcatcca                                                                                              <  1:590882/140‑1 (MQ=255)
                                        cggATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTcatccat                                                                                             <  1:210425/140‑1 (MQ=255)
                                        cggATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTcatccat                                                                                             <  2:149927/140‑1 (MQ=255)
                                          gATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCa                                                                                           <  2:238182/140‑1 (MQ=255)
                                          gATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCa                                                                                           <  1:586286/140‑1 (MQ=255)
                                             ttATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCa                                                                                           <  1:1004823/137‑1 (MQ=255)
                                             ttATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCa                                                                                           >  2:1004823/1‑137 (MQ=255)
                                               aTGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACg                                                                                      <  1:13018/140‑1 (MQ=255)
                                               aTGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACg                                                                                      <  1:909615/140‑1 (MQ=255)
                                                    cGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTa                                                                                 <  1:549994/140‑1 (MQ=255)
                                                    cGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGAGGTa                                                                                 >  1:625155/1‑140 (MQ=255)
                                                       aGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACg                                                                              >  1:1120505/1‑140 (MQ=255)
                                                            aTGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAg                                                                         <  2:761215/140‑1 (MQ=255)
                                                              ggCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAgg                                                                        <  1:497890/139‑1 (MQ=255)
                                                              ggCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAgg                                                                        >  2:497890/1‑139 (MQ=255)
                                                              ggCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGAGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAGGa                                                                       >  1:1084342/1‑140 (MQ=255)
                                                                         gAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAGGATCGCTGAtgtg                                                            <  1:610197/140‑1 (MQ=255)
                                                                         gAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAGGATCGCTGAtgtg                                                            >  2:610197/1‑140 (MQ=255)
                                                                              gCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAg                                                                         >  1:339238/1‑122 (MQ=255)
                                                                              gCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAg                                                                         <  2:339238/122‑1 (MQ=255)
                                                                                        tCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAGGATCGCTGATGTGCTGGCGCAGGTTGAg                                             <  1:532582/140‑1 (MQ=255)
                                                                                                 tGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAGGATCGCTGATGTGCTGGCGCAGGTTGAGCGCCACCAt                                    >  2:129321/1‑140 (MQ=255)
                                                                                                    cAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAGGATCGCTGATGTGCTGGCGCAGGTTGAGCGCCACCATGCg                                 <  1:129321/140‑1 (MQ=255)
                                                                                                     agagTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAGGATCGCTGATGTGCTGGCGCAGGTTGAGCGCCACCATGCgg                                <  1:434589/140‑1 (MQ=255)
                                                                                                       agTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAGGATCGCTGATGTGCTGGCGCAGGTTGAGCGCCACCATGCGgtc                              >  1:204656/1‑140 (MQ=255)
                                                                                                       agTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAGGATCGCTGATGTGCTGGCGCAGGTTGAGCGCCACCATGCGgtc                              <  1:4252/140‑1 (MQ=255)
                                                                                                       agTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAGGATCGCTGATGTGCTGGCGCAGGTTGAGCGCCACCATGCGgtc                              <  1:174649/140‑1 (MQ=255)
                                                                                                       agTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAGGATCGCTGATGTGCTGGCGCAGGTTGAGCGCCACCATGCGgtc                              >  2:439313/1‑140 (MQ=255)
                                                                                                       agTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAGGATCGCTGATGTGCTGGCGCAGGTTGAGCGCCACCATGCGgtc                              >  1:79223/1‑140 (MQ=255)
                                                                                                         tGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGCTCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAGGATCGCTGATGTGCTGGCGCAGGTTGAGCGCCACCATGCGgtcgt                            <  2:999531/140‑1 (MQ=255)
                                                                                                         tGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAGGATCGCTGATGTGCTGGCGCAGGTTGAGCGCCACCATGCGgtcgt                            >  1:154736/1‑140 (MQ=255)
                                                                                                         tGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAGGATCGCTGATGTGCTGGCGCAGGTTGAGCGCCACCATGCGgtcgt                            <  2:660771/140‑1 (MQ=255)
                                                                                                         tGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAGGATCGCTGATGTGCTGGCGCAGGTTGAGCGCCACCATGCGgtcgt                            >  1:328526/1‑140 (MQ=255)
                                                                                                                 cgTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAGGATCGCTGATGTGCTGGCGCAGGTTGAGCGCCACCATGCGgtcgtc                           <  2:999432/133‑1 (MQ=255)
                                                                                                                 cgTACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAGGATCGCTGATGTGCTGGCGCAGGTTGAGCGCCACCATGCGgtcgtc                           >  1:999432/1‑133 (MQ=255)
                                                                                                                   tACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAGGATCGCTGATGTGCTGGCGCAGGTTGAGCGCCACCATGCGGTCGTCTTTGTtata                  >  2:165142/1‑140 (MQ=255)
                                                                                                                   tACCGACGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAGGATCGCTGATGTGCTGGCGCAGGTTGAGCGCCACCATGCGGTCGTCTTTGTtata                  >  2:694003/1‑140 (MQ=255)
                                                                                                                        aCGGTGAGTTCACAGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAGGATCGCTGATGTGCTGGCGCAGGTTGAGCGCCACCATGCGGTCGTCTTTGTTATAGCTGt             >  2:687551/1‑140 (MQ=255)
                                                                                                                               gTTCACAGGCCTCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAGGATCGCTGATGTGCTGGCGCAGGTTGAGCGCCACCATGTGGTCGTCTTTGTTATAGCTGTCAATGcc      <  1:27308/140‑1 (MQ=255)
                                                                                                                                   acaGGCCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAGGATCGCTGATGTGCTGGCGCAGGTTGAGCGCCACCATGCGGTCGTCTTTGTTATAGCTGTCAATGCCCGTa  <  1:439313/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
ATTGCCATCTCCGGTAGTCCCAGCTGCTGCATAATCCCGGCGGATTTATGTTCGTAGTTGATGGCAATTGCCGGAGTGGCAAAATTCATCGAGATAATGGCAGAGTGCAGGCGCGTACCGACGGTGAGTTCACAGGCCCCCAGAATTTTGCCCATTTCCAGATCGTTGAGTTCATCCATCACTACGTGGTAACGGGCAGGATCGCTGATGTGCTGGCGCAGGTTGAGCGCCACCATGCGGTCGTCTTTGTTATAGCTGTCAATGCCCGTA  >  NC_000913/2117290‑2117559

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