Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 371,588 (G)6→7 coding (350/867 nt) mhpC → 2‑hydroxy‑6‑ketonona‑2,4‑dienedioate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913371,5821.G98.5% 230.2 / ‑3.4 65coding (344/867 nt)mhpC2‑hydroxy‑6‑ketonona‑2,4‑dienedioate hydrolase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (36/28);  total (37/28)
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.41e-01

ATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑‑GGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAGGTATTAAGCGACTGAAT  >  NC_000913/371443‑371715
                                                                                                                                            |                                                                                                                                       
aTCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAt                                                                                                                                          >  2:559553/1‑140 (MQ=255)
aTCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAt                                                                                                                                          >  1:275147/1‑140 (MQ=255)
   ctgctgGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑‑ggg                                                                                                                                     >  1:999248/1‑140 (MQ=255)
   ctgctgGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCATTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑‑ggg                                                                                                                                     <  2:246996/140‑1 (MQ=255)
       tgGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑‑ggggggg                                                                                                                                <  1:517864/140‑2 (MQ=255)
           ttGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑‑gggggggccat                                                                                                                            >  2:883745/1‑135 (MQ=255)
           ttGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑‑gggggggccat                                                                                                                            >  1:474484/1‑135 (MQ=255)
              tCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGt                                                                                                                          >  1:1092973/1‑140 (MQ=255)
                   gTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑‑gggg                                                                                                                                    >  2:933025/1‑125 (MQ=255)
                   gTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑‑gggg                                                                                                                                    <  1:933025/125‑1 (MQ=255)
                        ggCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGt                                                                                                                <  2:312472/140‑1 (MQ=255)
                        ggCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGt                                                                                                                >  2:116403/1‑140 (MQ=255)
                         gCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGtt                                                                                                               <  1:151052/140‑1 (MQ=255)
                         gCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGtt                                                                                                               <  2:1092973/140‑1 (MQ=255)
                                  aTTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTaa                                                                                                      >  1:1138749/1‑140 (MQ=255)
                                  aTTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTaa                                                                                                      >  2:412647/1‑140 (MQ=255)
                                            tAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAg                                                                                            <  1:116403/140‑1 (MQ=255)
                                              aTAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAgcg                                                                                          <  2:208870/140‑1 (MQ=255)
                                                  tGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTc                                                                                      >  1:532184/1‑140 (MQ=255)
                                                    gTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCgg                                                                                    >  1:997135/1‑140 (MQ=255)
                                                    gTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCgg                                                                                    >  1:1173096/1‑140 (MQ=255)
                                                       cgcgATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCaa                                                                                 <  1:559553/140‑1 (MQ=255)
                                                        gcgATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGTAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCaaa                                                                                <  2:1138749/140‑1 (MQ=255)
                                                        gcgATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCaaa                                                                                <  2:532184/140‑1 (MQ=255)
                                                          gATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑‑GGGGAGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACt                                                                              >  1:838764/1‑140 (MQ=255)
                                                            tCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTgg                                                                            >  1:716362/1‑140 (MQ=255)
                                                                      aTGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAAc                                                                               >  2:824185/1‑127 (MQ=255)
                                                                      aTGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAAc                                                                               <  1:824185/127‑1 (MQ=255)
                                                                      aTGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATggg                                                                  >  1:1073156/1‑140 (MQ=255)
                                                                      aTGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATggg                                                                  <  1:883745/140‑1 (MQ=255)
                                                                          aCGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGggcggc                                                              >  2:549601/1‑140 (MQ=255)
                                                                          aCGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGggcggc                                                              >  2:345877/1‑140 (MQ=255)
                                                                          aCGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGggcggc                                                              <  2:999248/140‑1 (MQ=255)
                                                                          aCGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGggcggc                                                              >  1:726015/1‑140 (MQ=255)
                                                                           cGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGgcggcg                                                             >  2:641480/1‑140 (MQ=255)
                                                                           cGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGgcggcg                                                             >  2:930727/1‑140 (MQ=255)
                                                                               tCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTAc                                                         >  1:1157894/1‑140 (MQ=255)
                                                                                 cTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACgg                                                       <  2:1073156/140‑1 (MQ=255)
                                                                                    aaaaGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGggcg                                                    >  1:594951/1‑140 (MQ=255)
                                                                                     aaaGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGggcgg                                                   <  1:549601/140‑1 (MQ=255)
                                                                                      aaGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGggcggc                                                  <  2:997135/140‑1 (MQ=255)
                                                                                         cGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGggcgg                                                               >  1:632497/1‑124 (MQ=255)
                                                                                         cGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGggcgg                                                               <  2:632497/124‑1 (MQ=255)
                                                                                         cGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGggcg                                                    <  2:1048257/135‑1 (MQ=255)
                                                                                         cGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGggcg                                                    >  1:1048257/1‑135 (MQ=255)
                                                                                         cGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATg                                               >  2:317137/1‑140 (MQ=255)
                                                                                           tggtggATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGa                                              >  1:997715/1‑140 (MQ=255)
                                                                                            ggtggATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAgt                                            <  2:570124/140‑1 (MQ=255)
                                                                                            ggtggATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAgt                                            <  1:641480/140‑1 (MQ=255)
                                                                                                  tCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAgtttgttt                                      >  2:1104168/1‑140 (MQ=255)
                                                                                                   cAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTa                                     >  2:363200/1‑140 (MQ=255)
                                                                                                          atatCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGAt                              >  1:1170072/1‑140 (MQ=255)
                                                                                                           tatCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATg                             >  2:152289/1‑140 (MQ=255)
                                                                                                              cGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAAGTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCa                          >  2:1067958/1‑140 (MQ=255)
                                                                                                              cGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAgt                                            <  1:696913/122‑1 (MQ=255)
                                                                                                              cGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAgt                                            >  2:696913/1‑122 (MQ=255)
                                                                                                              cGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCa                          <  1:1104168/140‑1 (MQ=255)
                                                                                                                ccAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAAc                        >  1:1066001/1‑140 (MQ=255)
                                                                                                                  aaaaTCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCg                      <  1:345877/140‑1 (MQ=255)
                                                                                                                  aaaaTCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCg                      <  2:1170072/140‑1 (MQ=255)
                                                                                                                    aaTCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGaa                    >  2:1037982/1‑140 (MQ=255)
                                                                                                                      tCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAgg                  <  1:1037982/140‑1 (MQ=255)
                                                                                                                            tgctgGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAGGTATTaa            <  1:502278/140‑1 (MQ=255)
                                                                                                                             gctgGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCg                      <  1:258082/129‑1 (MQ=255)
                                                                                                                             gctgGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCg                      >  2:258082/1‑129 (MQ=255)
                                                                                                                              ctgGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAGGTATTAAGc          >  1:816950/1‑140 (MQ=255)
                                                                                                                                gggCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCg                      <  2:405172/126‑1 (MQ=255)
                                                                                                                                gggCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCg                      >  1:405172/1‑126 (MQ=255)
                                                                                                                                gggCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAGGTATTAAGCGa        <  1:1153404/140‑1 (MQ=255)
                                                                                                                                gggCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAGGTATTAAGCGa        >  2:1202491/1‑140 (MQ=255)
                                                                                                                                   cAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAGGTATTAAGCGACTg     >  1:1174767/1‑140 (MQ=255)
                                                                                                                                      cTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAGGTATTAAGCGAc       >  2:929367/1‑135 (MQ=255)
                                                                                                                                      cTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAGGTATTAAGCGAc       <  1:929367/135‑1 (MQ=255)
                                                                                                                                      cTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAGGTATTAAGCGACTGAAt  >  1:539163/1‑140 (MQ=255)
                                                                                                                                            |                                                                                                                                       
ATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑‑GGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAGGTATTAAGCGACTGAAT  >  NC_000913/371443‑371715

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