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

GCCACTGGTTATCTTTACCGTATTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTG  >  NC_000913/1663632‑1664085
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                    
gCCACTGGTTATCTTTACCGTATTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGGGCTTCTGGTCTGa                                                                                                                                                                                                                                      >  1:483149/1‑227 (MQ=255)
        ttATCTTTACCGTATTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAAATCCCTGAATCGAATCGGAGCTTCTGGTCTGGGTAACggga                                                                                                                                                                                                                             >  1:391207/1‑228 (MQ=255)
                  cGTATTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGACAAAAATGACGCAAGCGGTCAGCGTATTGTCCGTGGCATGTTTTTCATCTGGCTGTTGATGGGCGTCGGGTTTATTGCATCTGTCATGCACCTGGGATCGCCACTGGGGGCGTTTAATTCCCTGACTCGCATCGGGCGTTTCTGTCTGGGGAACGAGATTGCTccccg                                                                                                                                                                                                                   >  2:308045/1‑226 (MQ=255)
                         ggTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCCCTCCGGGCGTTTAACTCCCTGAATCGCATCGTGGCTTCTGGCTTGAGTAACGCAATTGCTGCCGGTTCCata                                                                                                                                                                                                           >  1:664777/1‑228 (MQ=255)
                         ggTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCata                                                                                                                                                                                                           >  1:14658/1‑228 (MQ=255)
                         ggTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGGGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATTGCCCCTGCGGGCGTTTAACTCCCTGAATCGAATCGGGGCTTCTGGTCTGAGAAACGAGATTGCTTCCGGTTCCaaa                                                                                                                                                                                                            >  1:314019/1‑226 (MQ=255)
                              aTGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCCCCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTTGGTGACGAGATTGCTGCGGGTTCCATAttttt                                                                                                                                                                                                      >  1:296209/1‑228 (MQ=255)
                              aTGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGGGATTGCTGCCGGTTCCATAttttt                                                                                                                                                                                                      >  1:625094/1‑228 (MQ=255)
                               tGCGTAGTTGGCGCGCTGATTGTTAGCGTAATTGGCTGGTTCGCCGTCAAAAATGACGCAGACAGTCAGCGTACTGTCCGCGGCCTGTTTTTCCTCTGGCTGTTGATGGGCCTCGGGATTAGTGCCCCTTTCATACACCCGGGATACCCACTGGTGGGGTTTAAATCCCCAACTCGAATCGGAGCTTCTGGTCTGAGTAACGCAAGTGCCCGAGGTTCCAtttttttt                                                                                                                                                                                                     >  2:133288/1‑228 (MQ=255)
                               tGCGTAGTAGGCGCGCTGATTGTTAGCGGACTTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCCCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCTTGCACCTGGGATAGCCACAGCGGGCGTTTAATTCCCTCATTCGAATCGGAGTCTCCGGCCTGCTCACCGAGACTCCTTCCGGTTCCTTAtttttt                                                                                                                                                                                                     >  2:466857/1‑228 (MQ=255)
                                   tagtagGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTGCCTGAATCGGATCCGAGCTTCTGTCCTTGGTAACGGGGTGTCTGCGGGGTCCACTAtattttcct                                                                                                                                                                                                  >  2:381983/1‑225 (MQ=255)
                                    agtagGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGGGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGt                                                                                                                                                                                                >  1:488655/1‑228 (MQ=255)
                                       tggtgtgcTGTTTTTTAGGGGAATGGGCTGGTTCGCCGCCAAAAATGGCGCAGACCGGCAGCGTATTGGCCGTGGCCTGTTTTGCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCGGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTg                                                                                                                                                                                              <  2:175084/221‑1 (MQ=255)
                                       agGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTg                                                                                                                                                                                              <  1:194845/227‑1 (MQ=255)
                                       agGCGCGCGGATTGTTAGCGGAATTGGCCGGTTCGCCGTCAAAAATGATGCAGTCCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTg                                                                                                                                                                                              <  1:466857/227‑1 (MQ=255)
                                           gcgcTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGGGATTGCTGCCGGTTCCATATTTTTTGCCGtt                                                                                                                                                                                               >  1:175084/1‑222 (MQ=255)
                                                gATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTCTTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGTATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTAt                                                                                                                                                                                     >  1:549087/1‑227 (MQ=255)
                                                           gAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAAGGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGCAGGGTTAtggtggctgttgg                                                                                                                                                                         >  1:483416/1‑228 (MQ=255)
                                                           gAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGCCTGCGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTg                                                                                                                                                                                              >  1:74923/1‑207 (MQ=255)
                                                            atttGGCTTGTTCGCCGCCAAAAATGACGCGGACCGTCAGCGTGTTGGCCGTGGCATGTTCTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGTCGTTTATCTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTgg                                                                                                                                                                                             <  2:74923/205‑1 (MQ=255)
                                                             aTTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTTAGTAACGTGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGCTTTTGGGGGCTGGTGGCa                                                                                                                                                                       >  2:413236/1‑228 (MQ=255)
                                                                          ccggcAAAAATGGCGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTaaaa                                                                                                                                                          <  1:593881/228‑1 (MQ=255)
                                                                               aaaaaTGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATTGCCACGGGGGGCGTTTAACTCCCCTAATCGAATCGGCGCTTCTGGTCGGAGTAACGACATTGCTGCCGGTTCACTATTTTTTGCCGTTGTCGGATTATGTGTGATGGTGGGAGTCATCGTTgaaagacgg                                                                                                                                                     >  2:64617/1‑223 (MQ=255)
                                                                                               gtcgtctttttGTCCTTGGCATGTTTTTCCCGTGGCTGTTGATGGGCGTCGGGTTTATTGCCTTTGTCATGCCCCTGGGATTGGCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAAc                                                                                                                                     <  2:255725/219‑1 (MQ=255)
                                                                                                     ttttttttcttGGCTTTTTTTTCCTCGGGCTGTGGATGGGCGTGGGGTTTATTCCCTTCGGCATGCCCCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGt                                                                                                                               <  2:14658/218‑1 (MQ=255)
                                                                                                           tCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGGGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGACTGCGTAACGGAATTGCCGCCGGTTCCATATTTTTTGCCGTTGGCGGATAATGGGGGCTGGGGGCAGTCCTCGGTTAAATGCCGGAACAGTTGGGGAACATCGGGTTACtt                                                                                                                          >  2:69615/1‑227 (MQ=255)
                                                                                                           tCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGTGTGCTGGTGCCAGTCAGCGGTAAAATGCCCCAACCGTTGGGGAAACTCTGGTTACtt                                                                                                                          >  1:369149/1‑227 (MQ=255)
                                                                                                           tCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACtt                                                                                                                          >  1:452573/1‑227 (MQ=255)
                                                                                                              gTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTCATTTTCa                                                                                                                      >  1:359930/1‑228 (MQ=255)
                                                                                                                ggCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCCAGCGTTGGGGAAACTCTGGTTACTTTTTAGc                                                                                                                    >  1:144254/1‑228 (MQ=255)
                                                                                                                            cttggggtGTTGAGGGGGGTCGGTTTTATTGCCCCTGGCATGCCCCTGGGTTCGCCCCTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGc                                                                                                        <  2:20315/221‑1 (MQ=255)
                                                                                                                                ggCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATCTTTTTTGCCCTTGGCGGCTTATGGTGGCTGGCGGGGGTCTTTGGTAAACTGCCCAACGCGGTGGGGGAATTCTGGTTAATTTTCCGCATGGGGCCCcgcggcc                                                                                                    >  2:92729/1‑225 (MQ=255)
                                                                                                                                  cgtttgttgGGCGGCGGGTTTATTGCCTCGGCTATGCACCGGGGTGGGCCTGTGCGGGCGTTTAACTCCCTGAAACGAATCGGGGCTTCTGGTCTGGGTAACGAGTTTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCAtt                                                                                                  <  2:695413/225‑1 (MQ=255)
                                                                                                                                       gaggtgcttCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGt                                                                                             <  1:164784/220‑1 (MQ=255)
                                                                                                                                       gATGGGCGTCGGGTTTATTGCCTATGTCATGCCCCTGGGATCGCCACTGCGGGCGGTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGATTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGt                                                                                             <  1:409561/228‑1 (MQ=255)
                                                                                                                                         tGGGCGTCGGGTTTATTGCCTCTGTTATGGCTCTGGGCTCGGCACTGGGGGCGTTTACCTCCCTGAATCGAATCGGAGCTTCTGGTCTGATTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCt                                                                                           <  1:308045/228‑1 (MQ=255)
                                                                                                                                         gtggggtggggTTTATTGCCTCTGTCCGTCCCCTGGGATCGCCACTGGGGGCGTTTAACTCCCTGAATCGAATCGGGGCTTCTGGGTTGGGTAAAGGGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGGTGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCt                                                                                           <  2:483149/220‑1 (MQ=255)
                                                                                                                                           ggCGTCGGGTTTATTGCCTCTGTCATGCACTTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGGCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCGAGCGTTGGGGAAATTCTGGTTGATTTTCAGCATGGCGTTCGGCGGCATTTTGGTCTgg                                                                                         >  1:645976/1‑228 (MQ=255)
                                                                                                                                           ggCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTgg                                                                                         <  1:92729/228‑1 (MQ=255)
                                                                                                                                               tCGGGTTTATTGCCTCTGTCATGCACCTTGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGCGTAACGAGATTGCTGCCGGTTCCATCTTTTTTGCCGTTGGCGGATTATGGAGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGGAACTCTGGGTCATTTTAACCAAGGCGCCCGACGCCATTTTGGCCTGgttga                                                                                     >  2:167013/1‑228 (MQ=255)
                                                                                                                                               tCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCAACGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCAAGGCGCTCGGCGTCATTTTCGTCTGgatgc                                                                                     >  1:375420/1‑227 (MQ=255)
                                                                                                                                               tCGGGGTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCCGAGCTTCTGGTCTGAGGACCGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGAGTATGGTGGCTGGTGGCAGACATCGGCAAAATGCCGCAGGCGTTGGGGGAACTCTGGTTACTTTTCACAATGCGGCCCGGGGTAATTTTCGGCTggtgga                                                                                     >  2:164784/1‑224 (MQ=255)
                                                                                                                                                        ttccctcttgtaTGCCCCGGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCGGGTCTGAGTAGCGGGAGTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGc                                                                             <  2:369149/216‑1 (MQ=255)
                                                                                                                                                             tctGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGCCGTCAttt                                                                                                 >  1:186968/1‑202 (MQ=255)
                                                                                                                                                              ttgttcggccccTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTAt                                                                       <  1:413236/217‑1 (MQ=255)
                                                                                                                                                              ttGTCATGCCCCGGGGATCGCCCCTTCGGGCGTTTAACTCCCTGAATCGAATCGGGTCTTCTGGTCTGAGTAACGGGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTAt                                                                       <  2:452573/226‑1 (MQ=255)
                                                                                                                                                              ctGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCCTCGGTAAAATGCCGCAAGCGTTTGGGAAACTCTGGTTACTTTTCAGCATGGCGCCCCGCGCCATTTTTGTCTGGATGATGACCTGCGTGTATc                                                                      >  2:358386/1‑228 (MQ=255)
                                                                                                                                                              ctGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCTTTTTCGGCTGGATATTGACCttccggga                                                                        >  2:368170/1‑219 (MQ=255)
                                                                                                                                                              ctGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTCTCGTCTGGATGATGACCTGCGTGTAt                                                                       <  1:561984/227‑1 (MQ=255)
                                                                                                                                                              aggTCATGCCCCTGGGTTCGCCACGGCGGGCTTTTAACCCCCTGAATTGAATCGGAGCTTCTGGTCGGGGTAACGGGGTTGGGGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTAt                                                                       <  2:476679/225‑1 (MQ=255)
                                                                                                                                                                       cccggggATCGCCCCTGCGGGCGTTTAACTCCTTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGaca                                                             <  1:381983/224‑1 (MQ=255)
                                                                                                                                                                        ccTGGGTTCGCCACTGCGGGCGTTTAACGCCCTGAATCGAATCGGATCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGacac                                                            <  1:69615/228‑1 (MQ=255)
                                                                                                                                                                         cTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGcc                                                                                                                                                                                                                     >  2:670447/1‑75 (MQ=255)
                                                                                                                                                                          tGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCg                                                                                                                                                                                                                    <  1:670447/75‑1 (MQ=255)
                                                                                                                                                                           gggATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTa                                                                        >  1:476679/1‑213 (MQ=255)
                                                                                                                                                                           gggATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGt                                                         >  1:565263/1‑228 (MQ=255)
                                                                                                                                                                           gggATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGAGGATGACCTTCGTGTTTCAAATCgcatccgt                                                         >  1:249830/1‑221 (MQ=255)
                                                                                                                                                                           gggATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCCCGGCGTCCTTTTGGTCTGGTGGGTGACCTGCGTGTTTAAAATctggctct                                                          >  2:166318/1‑220 (MQ=255)
                                                                                                                                                                            ggATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAAct                                                                                                                                    <  1:144151/152‑1 (MQ=255)
                                                                                                                                                                            ggATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGc                                                                             <  1:166318/207‑1 (MQ=255)
                                                                                                                                                                                 gCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGCAATTTTCGTCTGGATGGTTACTTTGGTGTATCAAAACGAACCCGTGCCacc                                                   >  2:687356/1‑226 (MQ=255)
                                                                                                                                                                                 gCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGGGGCTGGTGGCAGTCATCGGTACAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTTCGTCTGTCAAATCGACACCGGTCCAAc                                                   >  2:115288/1‑228 (MQ=255)
                                                                                                                                                                                  ccACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGTTGACCTGCGTGTATCAAATCGACACCTGGCCAAcc                                                  >  1:183485/1‑228 (MQ=255)
                                                                                                                                                                                  ccACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACCTTGGTTACTTTTCAGCATGGCGCCAGGCGTCATTTTCGGTTGGATGGTTACCTGCGGGGTTAAAAGCGACACCGTGCCAAcc                                                  >  2:32888/1‑228 (MQ=255)
                                                                                                                                                                                        cGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGGCCCAACCTGCCAt                                            >  1:370165/1‑228 (MQ=255)
                                                                                                                                                                                         ggggGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCGGAGTAACGGGATTGCGGCCGGTTCCATATTTTTTGCCGTTGGCGGATTGTGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCAt                                            <  2:664777/227‑1 (MQ=255)
                                                                                                                                                                                         gggCGTTTAACTCCCTGAATCGATGCGGAGCTTCGTGTCTGAGTTACTAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTGCTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCAt                                            <  1:687356/227‑1 (MQ=255)
                                                                                                                                                                                         cggtgTTTAACTCCCGTGATTGAATTGGAGCTTCTTGTCTGGGTAACGAGATTGCTGGCGGTTCCATATTTTTTGCCGGTGGCGGAGTATGGTGGCTTGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCAt                                            <  1:133288/223‑1 (MQ=255)
                                                                                                                                                                                          ggCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAACCTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGAACCCGTGCCAACCTGTCATa                                           >  1:143849/1‑227 (MQ=255)
                                                                                                                                                                                           gctttttACCCCCTTAATTGAATCGTGTCTTCTGGTGTGGGTAACGGGATTGCTGTCGGTTCCATATTTTTTGCCGTTGGGGGATTATGGTGGCTGGTGTCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATaa                                          <  2:382600/224‑1 (MQ=255)
                                                                                                                                                                                             ttttATCTCCTTGAATCGATTCGGAGCTTCTGGTCTGGGTAACGGGATTGCTGCCGGTTCCATTTTTTTTGCCGGTGGCGGATTATGGCGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGAAGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACgg                                       <  2:625094/227‑1 (MQ=255)
                                                                                                                                                                                             gTTTATCTCCCTGTTTTGATTCGGGGTTTTTGGTCTGGGTAAGGGGATTGCTGCGGGGTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACgg                                       <  2:359930/228‑1 (MQ=255)
                                                                                                                                                                                                       cTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAACCGACACCGTGCCAACCTGGCATAACGGTTCCACAAc                              >  1:693048/1‑227 (MQ=255)
                                                                                                                                                                                                         gaatcgaatcGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTGGTCTGGCTGATGACTTGCGTGGTTCAAATCGACACCCTGCCAACCTGGCATAACGGttcaaccaccca                           >  2:88669/1‑218 (MQ=255)
                                                                                                                                                                                                                      atcTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCt               <  1:336916/226‑1 (MQ=255)
                                                                                                                                                                                                                       gCTTTTGGTCTGAGTAACGTGTTGTCTGCTGGTTCCTTATTGTTTGCCGTTGGCGGATTATGGTGACTGGTGGCATGTATCGGGAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGTATGGCGCTCGGCGTCATTTTCGTCTGGATGATGGCCTGCGTGTATCAAATCGACACCGTGCCATCCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCt               <  1:9688/227‑1 (MQ=255)
                                                                                                                                                                                                                                   aGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCCAATCGACACCGTGCCAACCCGGCATACCGGTTCAACCACGCTTGCTTTCTTCCTTACgggtttcttt  >  2:77800/1‑220 (MQ=255)
                                                                                                                                                                                                                                   aGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCTTAACGGTTACACTACGCTGGCTTTATTCCTGACGGTAttgtt   >  1:20360/1‑227 (MQ=255)
                                                                                                                                                                                                                                   aGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAAATCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTAtggttg  >  1:96673/1‑228 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                    
GCCACTGGTTATCTTTACCGTATTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGTTACTTTTCAGCATGGCGCTCGGCGTCATTTTCGTCTGGATGATGACCTGCGTGTATCAAATCGACACCGTGCCAACCTGGCATAACGGTTACACTACGCTGGCTTTCTTCCTGACGGTATTGTTG  >  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: