Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,663,858 A→G 14.8% S84G (AGT→GGT)  ynfH → oxidoreductase, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,663,8580AG14.8% 158.7 / 2.1 75S84G (AGT→GGT) ynfHoxidoreductase, membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (34/25);  new base G (5/7);  total (43/32)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GCCACTGGTTATCTTTACCGTATTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTG  >  NC_000913/1663632‑1664085
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                    
gCCACTGGTTATCTTTACCGTATTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGGGCTTCTGGTCTGa                                                                                                                                                                                                                                      >  1:483149/1‑227 (MQ=255)
        ttATCTTTACCGTATTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAAATCCCTGAATCGAATC‑GGAGCTTCTGGTCTGGGTAACggga                                                                                                                                                                                                                             >  1:391207/1‑228 (MQ=255)
                  cGTATTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGACAAAAATGACGCAAGCGGTCAGCGTATTGTCCGTGGCATGTTTTTCATCTGGCTGTTGATGGGCGTCGGGTTTATTGCATCTGTCATGCACCTGGGATCGCCACTGGGGGCGTTTAATTCCCTGACTCGCATCGGGCGTTTCT‑GTCTGGGGAACGAGATTGCTccccg                                                                                                                                                                                                                   >  2:308045/1‑226 (MQ=255)
                         ggTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCCCTCCGGGCGTTTAACTCCCTGAATCGCATC‑GTGGCTTCTGGCTTGAGTAACGCAATTGCTGCCGGTTCCa‑ta                                                                                                                                                                                                           >  1:664777/1‑228 (MQ=255)
                         ggTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCa‑ta                                                                                                                                                                                                           >  1:14658/1‑228 (MQ=255)
                         ggTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGGGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATTGCCCCTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGGGCTTCTGGTCTGAGAAACGAGATTGCTTCCGGTTCCaaa                                                                                                                                                                                                            >  1:314019/1‑226 (MQ=255)
                              aTGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCCCCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTTGGTGACGAGATTGCTGCGGGTTCCA‑TAttttt                                                                                                                                                                                                      >  1:296209/1‑228 (MQ=255)
                              aTGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGGGATTGCTGCCGGTTCCA‑TAttttt                                                                                                                                                                                                      >  1:625094/1‑228 (MQ=255)
                               tGCGTAGTTGGCGCGCTGATTGTTAGCGTAATTGGCTGGTTCGCCGTCAAAAATGACGCAGACAGTCAGCGTACTGTCCGCGGCCTGTTTTTCCTCTGGCTGTTGATGGGCCTCGGGATTAGTGCCCCTTTCATACACCCGGGATACCCACTGGTGGGGTTTAAATCCCCAACTCGAATC‑GGAGCTTCTGGTCTGAGTAACGCAAGTGCCCGAGGTTCCA‑tttttttt                                                                                                                                                                                                     >  2:133288/1‑228 (MQ=255)
                               tGCGTAGTAGGCGCGCTGATTGTTAGCGGACTTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCCCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCTTGCACCTGGGATAGCCACAGCGGGCGTTTAATTCCCTCATTCGAATC‑GGAGTCTCCGGCCTGCTCACCGAGACTCCTTCCGGTTCCT‑TAtttttt                                                                                                                                                                                                     >  2:466857/1‑228 (MQ=255)
                                   tagtagGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTGCCTGAATCGGATC‑CGAGCTTCTGTCCTTGGTAACGGGGTGTCTGCGGGGTCCACTAtattttcct                                                                                                                                                                                                  >  2:381983/1‑225 (MQ=255)
                                    agtagGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGGGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGt                                                                                                                                                                                                >  1:488655/1‑228 (MQ=255)
                                       tggtgtgcTGTTTTTTAGGGGAATGGGCTGGTTCGCCGCCAAAAATGGCGCAGACCGGCAGCGTATTGGCCGTGGCCTGTTTTGCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCGGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTg                                                                                                                                                                                              <  2:175084/221‑1 (MQ=255)
                                       agGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTg                                                                                                                                                                                              <  1:194845/227‑1 (MQ=255)
                                       agGCGCGCGGATTGTTAGCGGAATTGGCCGGTTCGCCGTCAAAAATGATGCAGTCCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTg                                                                                                                                                                                              <  1:466857/227‑1 (MQ=255)
                                           gcgcTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGGGATTGCTGCCGGTTCCA‑TATTTTTTGCCGtt                                                                                                                                                                                               >  1:175084/1‑222 (MQ=255)
                                                gATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTCTTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGTATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTAt                                                                                                                                                                                     >  1:549087/1‑227 (MQ=255)
                                                           gAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAAGGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGCAGGGTTAtggtggctgttgg                                                                                                                                                                         >  1:483416/1‑228 (MQ=255)
                                                           gAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGCCTGCGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTg                                                                                                                                                                                              >  1:74923/1‑207 (MQ=255)
                                                            atttGGCTTGTTCGCCGCCAAAAATGACGCGGACCGTCAGCGTGTTGGCCGTGGCATGTTCTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGTCGTTTATCTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTgg                                                                                                                                                                                             <  2:74923/205‑1 (MQ=255)
                                                             aTTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTTAGTAACGTGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGCTTTTGGGGGCTGGTGGCa                                                                                                                                                                       >  2:413236/1‑228 (MQ=255)
                                                                          ccggcAAAAATGGCGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTaaaa                                                                                                                                                          <  1:593881/228‑1 (MQ=255)
                                                                               aaaaaTGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATTGCCACGGGGGGCGTTTAACTCCCCTAATCGAATC‑GGCGCTTCTGGTCGGAGTAACGACATTGCTGCCGGTTCAC‑TATTTTTTGCCGTTGTCGGATTATGTGTGATGGTGGGAGTCATCGTTgaaagacgg                                                                                                                                                     >  2:64617/1‑223 (MQ=255)
                                                                                               gtcgtctttttGTCCTTGGCATGTTTTTCCCGTGGCTGTTGATGGGCGTCGGGTTTATTGCCTTTGTCATGCCCCTGGGATTGGCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAAc                                                                                                                                     <  2:255725/219‑1 (MQ=255)
                                                                                                     ttttttttcttGGCTTTTTTTTCCTCGGGCTGTGGATGGGCGTGGGGTTTATTCCCTTCGGCATGCCCCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGt                                                                                                                               <  2:14658/218‑1 (MQ=255)
                                                                                                           tCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGGGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGACTGCGTAACGGAATTGCCGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATAATGGGGGCTGGGGGCAGTCCTCGGTTAAATGCCGGAACAGTTGGGGAACATCGGGTTACtt                                                                                                                          >  2:69615/1‑227 (MQ=255)
                                                                                                           tCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGTGTGCTGGTGCCAGTCAGCGGTAAAATGCCCCAACCGTTGGGGAAACTCTGGTTACtt                                                                                                                          >  1:369149/1‑227 (MQ=255)
                                                                                                           tCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACtt                                                                                                                          >  1:452573/1‑227 (MQ=255)
                                                                                                              gTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTCATTTTCa                                                                                                                      >  1:359930/1‑228 (MQ=255)
                                                                                                                ggCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCCAGCGTTGGGGAAACTCTGGTTACTTTTTAGc                                                                                                                    >  1:144254/1‑228 (MQ=255)
                                                                                                                            cttggggtGTTGAGGGGGGTCGGTTTTATTGCCCCTGGCATGCCCCTGGGTTCGCCCCTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGc                                                                                                        <  2:20315/221‑1 (MQ=255)
                                                                                                                                ggCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TCTTTTTTGCCCTTGGCGGCTTATGGTGGCTGGCGGGGGTCTTTGGTAAACTGCCCAACGCGGTGGGGGAATTCTGGTTAATTTTCCGCATGGGGCCCcgcggcc                                                                                                    >  2:92729/1‑225 (MQ=255)
                                                                                                                                  cgtttgttgGGCGGCGGGTTTATTGCCTCGGCTATGCACCGGGGTGGGCCTGTGCGGGCGTTTAACTCCCTGAAACGAATC‑GGGGCTTCTGGTCTGGGTAACGAGTTTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCAtt                                                                                                  <  2:695413/225‑1 (MQ=255)
                                                                                                                                       gaggtgcttCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGt                                                                                             <  1:164784/220‑1 (MQ=255)
                                                                                                                                       gATGGGCGTCGGGTTTATTGCCTATGTCATGCCCCTGGGATCGCCACTGCGGGCGGTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGATTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGt                                                                                             <  1:409561/228‑1 (MQ=255)
                                                                                                                                         tGGGCGTCGGGTTTATTGCCTCTGTTATGGCTCTGGGCTCGGCACTGGGGGCGTTTACCTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGATTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCt                                                                                           <  1:308045/228‑1 (MQ=255)
                                                                                                                                         gtggggtggggTTTATTGCCTCTGTCCGTCCCCTGGGATCGCCACTGGGGGCGTTTAACTCCCTGAATCGAATC‑GGGGCTTCTGGGTTGGGTAAAGGGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGGTGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCt                                                                                           <  2:483149/220‑1 (MQ=255)
                                                                                                                                           ggCGTCGGGTTTATTGCCTCTGTCATGCACTTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGGCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCGAGCGTTGGGGAAATTCTGGTTGATTTTCAGCATGGCGTTCGGCGGCATTTTGGTCTgg                                                                                         >  1:645976/1‑228 (MQ=255)
                                                                                                                                           ggCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTgg                                                                                         <  1:92729/228‑1 (MQ=255)
                                                                                                                                               tCGGGTTTATTGCCTCTGTCATGCACCTTGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGCGTAACGAGATTGCTGCCGGTTCCA‑TCTTTTTTGCCGTTGGCGGATTATGGAGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGGAACTCTGGGTCATTTTAACCAAGGCGCCCGACGCCATTTTGGCCTGgttga                                                                                     >  2:167013/1‑228 (MQ=255)
                                                                                                                                               tCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCAACGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCAAGGCGCTCGGCGTCATTTTCGTCTGgatgc                                                                                     >  1:375420/1‑227 (MQ=255)
                                                                                                                                               tCGGGGTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑CGAGCTTCTGGTCTGAGGACCGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGAGTATGGTGGCTGGTGGCAGACATCGGCAAAATGCCGCAGGCGTTGGGGGAACTCTGGTTACTTTTCACAATGCGGCCCGGGGTAATTTTCGGCTggtgga                                                                                     >  2:164784/1‑224 (MQ=255)
                                                                                                                                                        ttccctcttgtaTGCCCCGGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCGGGTCTGAGTAGCGGGAGTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGc                                                                             <  2:369149/216‑1 (MQ=255)
                                                                                                                                                             tctGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGCCGTCAttt                                                                                                 >  1:186968/1‑202 (MQ=255)
                                                                                                                                                              ttgttcggccccTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTAt                                                                       <  1:413236/217‑1 (MQ=255)
                                                                                                                                                              ttGTCATGCCCCGGGGATCGCCCCTTCGGGCGTTTAACTCCCTGAATCGAATC‑GGGTCTTCTGGTCTGAGTAACGGGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTAt                                                                       <  2:452573/226‑1 (MQ=255)
                                                                                                                                                              ctGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCCTCGGTAAAATGCCGCAAGCGTTTGGGAAACTCTGGTTACTTTTCAGCATGGCGCCCCGCGCCATTTTTGTCTGGATGATGACCTGCGTGTATc                                                                      >  2:358386/1‑228 (MQ=255)
                                                                                                                                                              ctGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCTTTTTCGGCTGGATATTGACCttccggga                                                                        >  2:368170/1‑219 (MQ=255)
                                                                                                                                                              ctGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTCTCGTCTGGATGATGACCTGCGTGTAt                                                                       <  1:561984/227‑1 (MQ=255)
                                                                                                                                                              aggTCATGCCCCTGGGTTCGCCACGGCGGGCTTTTAACCCCCTGAATTGAATC‑GGAGCTTCTGGTCGGGGTAACGGGGTTGGGGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTAt                                                                       <  2:476679/225‑1 (MQ=255)
                                                                                                                                                                       cccggggATCGCCCCTGCGGGCGTTTAACTCCTTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGaca                                                             <  1:381983/224‑1 (MQ=255)
                                                                                                                                                                        ccTGGGTTCGCCACTGCGGGCGTTTAACGCCCTGAATCGAATC‑GGATCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGacac                                                            <  1:69615/228‑1 (MQ=255)
                                                                                                                                                                         cTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGcc                                                                                                                                                                                                                     >  2:670447/1‑75 (MQ=255)
                                                                                                                                                                          tGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCg                                                                                                                                                                                                                    <  1:670447/75‑1 (MQ=255)
                                                                                                                                                                           gggATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTa                                                                        >  1:476679/1‑213 (MQ=255)
                                                                                                                                                                           gggATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGt                                                         >  1:565263/1‑228 (MQ=255)
                                                                                                                                                                           gggATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGAGGATGACCTTCGTGTTTCAAATCgcatccgt                                                         >  1:249830/1‑221 (MQ=255)
                                                                                                                                                                           gggATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCCCGGCGTCCTTTTGGTCTGGTGGGTGACCTGCGTGTTTAAAATctggctct                                                          >  2:166318/1‑220 (MQ=255)
                                                                                                                                                                            ggATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAAct                                                                                                                                    <  1:144151/152‑1 (MQ=255)
                                                                                                                                                                            ggATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGc                                                                             <  1:166318/207‑1 (MQ=255)
                                                                                                                                                                                 gCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGCAATTTTCGTCTGGATGGTTACTTTGGTGTATCAAAACGAACCCGTGCCacc                                                   >  2:687356/1‑226 (MQ=255)
                                                                                                                                                                                 gCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGGGGCTGGTGGCAGTCATCGGTACAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTTCGTCTGTCAAATCGACACCGGTCCAAc                                                   >  2:115288/1‑228 (MQ=255)
                                                                                                                                                                                  ccACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGTTGACCTGCGTGTATCAAATCGACACCTGGCCAAcc                                                  >  1:183485/1‑228 (MQ=255)
                                                                                                                                                                                  ccACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACCTTGGTTACTTTTCAGCATGGCGCCAGGCGTCATTTTCGGTTGGATGGTTACCTGCGGGGTTAAAAGCGACACCGTGCCAAcc                                                  >  2:32888/1‑228 (MQ=255)
                                                                                                                                                                                        cGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGGCCCAACCTGCCAt                                            >  1:370165/1‑228 (MQ=255)
                                                                                                                                                                                         ggggGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCGGAGTAACGGGATTGCGGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTGTGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCAt                                            <  2:664777/227‑1 (MQ=255)
                                                                                                                                                                                         gggCGTTTAACTCCCTGAATCGATGC‑GGAGCTTCGTGTCTGAGTTACTAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTGCTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCAt                                            <  1:687356/227‑1 (MQ=255)
                                                                                                                                                                                         cggtgTTTAACTCCCGTGATTGAATT‑GGAGCTTCTTGTCTGGGTAACGAGATTGCTGGCGGTTCCA‑TATTTTTTGCCGGTGGCGGAGTATGGTGGCTTGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCAt                                            <  1:133288/223‑1 (MQ=255)
                                                                                                                                                                                          ggCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAACCTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGAACCCGTGCCAACCTGTCATa                                           >  1:143849/1‑227 (MQ=255)
                                                                                                                                                                                           gctttttACCCCCTTAATTGAATC‑GTGTCTTCTGGTGTGGGTAACGGGATTGCTGTCGGTTCCA‑TATTTTTTGCCGTTGGGGGATTATGGTGGCTGGTGTCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATaa                                          <  2:382600/224‑1 (MQ=255)
                                                                                                                                                                                             ttttATCTCCTTGAATCGATTC‑GGAGCTTCTGGTCTGGGTAACGGGATTGCTGCCGGTTCCA‑TTTTTTTTGCCGGTGGCGGATTATGGCGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGAAGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACgg                                       <  2:625094/227‑1 (MQ=255)
                                                                                                                                                                                             gTTTATCTCCCTGTTTTGATTC‑GGGGTTTTTGGTCTGGGTAAGGGGATTGCTGCGGGGTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACgg                                       <  2:359930/228‑1 (MQ=255)
                                                                                                                                                                                                       cTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAACCGACACCGTGCCAACCTGGCATAACGGTTCCACAAc                              >  1:693048/1‑227 (MQ=255)
                                                                                                                                                                                                         gaatcgaatc‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTGGTCTGGCTGATGACTTGCGTGGTTCAAATCGACACCCTGCCAACCTGGCATAACGGttcaaccaccca                           >  2:88669/1‑218 (MQ=255)
                                                                                                                                                                                                                      atcTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCt               <  1:336916/226‑1 (MQ=255)
                                                                                                                                                                                                                       gCTTTTGGTCTGAGTAACGTGTTGTCTGCTGGTTCCT‑TATTGTTTGCCGTTGGCGGATTATGGTGACTGGTGGCATGTATCGGGAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGTATGGCGCTCGGCGTCATTTTCGTCTGGATGATGGCCTGCGTGTATCAAATCGACACCGTGCCATCCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCt               <  1:9688/227‑1 (MQ=255)
                                                                                                                                                                                                                                   aGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCCAATCGACACCGTGCCAACCCGGCATACCGGTTCAACCACGCTTGCTTTCTTCCTTACgggtttcttt  >  2:77800/1‑220 (MQ=255)
                                                                                                                                                                                                                                   aGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCTTAACGGTTACACTACGCTGGCTTTATTCCTGACGGTAttgtt   >  1:20360/1‑227 (MQ=255)
                                                                                                                                                                                                                                   aGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAAATCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTAtggttg  >  1:96673/1‑228 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                    
GCCACTGGTTATCTTTACCGTATTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATC‑GGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCA‑TATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTG  >  NC_000913/1663632‑1664085

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑