Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 12,484 C→A 20.3% Q108K (CAG→AAG)  dnaK → chaperone Hsp70, with co‑chaperone DnaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091312,4840CA20.3% 106.0 / 3.2 56Q108K (CAG→AAG) dnaKchaperone Hsp70, with co‑chaperone DnaJ
Reads supporting (aligned to +/- strand):  ref base C (41/3);  new base A (2/10);  total (43/13)
Fisher's exact test for biased strand distribution p-value = 4.69e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAG  >  NC_000913/12258‑12712
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                      
cGATCGCACCACGCCTTCTATCATTGCCTATGCCCAGGATGGTAAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCa                                                                                                                                                                                                                                       <  1:688979/228‑1 (MQ=255)
    cGCACGACGTCTTCTATCATTGCCTATATCCAGGATCGGGAAACTCTAGTTGGTCAGCCGGCTAAACGTCTGGCAGTGACGAACCACCAAAACACTCTGGTTGCGATTAAACGCCTGATTGGGCGCCGCGTCCAGGACCCAGAAGGACAGACGGAAGTTTCCATACTGCCGCTCAAAATTATTGCGTCGTATAACGGCGCCACAAGGGTCGAAGTTAAAAGCCACaaa                                                                                                                                                                                                                                   >  2:498516/1‑228 (MQ=255)
         cacGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAAAGGCa                                                                                                                                                                                                                            >  2:517683/1‑229 (MQ=255)
                        tGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAGACAATCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGGTTCCATCATGCCGTTCCAAATTATTGCTGCTGATAACGGCGACCCCTGGGTCGAAGTTAAAGTCCAGAAAATGGGACCGCCGGCAATTTct                                                                                                                                                                                                            >  2:331740/1‑229 (MQ=255)
                        tGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGAttt                                                                                                                                                                                                              >  1:258592/1‑227 (MQ=255)
                          ccTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGAGGCATGGGTCGAGGTTAAAGGCCGAAAAATGGCACCGCGCCAAATTTcggc                                                                                                                                                                                                          >  2:86078/1‑229 (MQ=255)
                            tataCCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTCCAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCCTGGGCCGAAGTTAAAGGCCAGAAAATGGCACCACCCCGGATTTctgcgg                                                                                                                                                                                                        >  2:189853/1‑227 (MQ=255)
                            tataCCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCCGATTTcttctg                                                                                                                                                                                                        >  1:392156/1‑229 (MQ=255)
                              taCCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGTTGATAACGGCGACGAATGTGTCGACGTCAAAGGACAGAAAATGGAAACCCCGCAGATATATGCTGaa                                                                                                                                                                                                      >  2:312712/1‑229 (MQ=255)
                                  cAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGTTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGCAGTGc                                                                                                                                                                                                  >  1:186169/1‑229 (MQ=255)
                                    ggATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGCTAACGGCAACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGAGGATTTCTGCTGAAGTGCTg                                                                                                                                                                                                >  1:63047/1‑229 (MQ=255)
                                     gATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTGGAGGTTAAAGGCCAGAAAATGGAACCGCCGCAGATTTCTGCTGAAGTGCTGa                                                                                                                                                                                               >  1:679708/1‑229 (MQ=255)
                                        ggTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCACCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCCTGCCGTTCCAAATTATTGCTGCTGATAACGGCGACGCCTGGGTCGAAGTTAAAAGCCAGAAAAGGGCACCTCCCGAGATTTTTGCTGaaaagggtaaac                                                                                                                                                                                            >  2:75483/1‑219 (MQ=255)
                                          tGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGaaaaa                                                                                                                                                                                           >  1:80683/1‑228 (MQ=255)
                                                  tAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCCGATTTTTGCTGAAGTGCTGAAAAAAAAaaagaa                                                                                                                                                                                  >  2:606135/1‑229 (MQ=255)
                                                   aGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGCCGCATGGCTCGAAGTTAAAGGCAAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAAAGaaaaaa                                                                                                                                                                                 >  2:723555/1‑229 (MQ=255)
                                                   aGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGACAATGGCACCGCCGCAGATTTCGGCTGAAGTGCTGAAAAAAAAAAGaaaa                                                                                                                                                                                >  2:597945/1‑229 (MQ=255)
                                                         caccgcGGCAAAACGCCAGGCAGTGCCGACCCCGCAAAACATTTGGTTGAGGATAAAACGCTTGGTTGGCCGCACGTTCCGGGACAAAAAAGTACAGGGTGAGTTTTCCAGCAGCCCGTACAAATTAATGGCGCTGAAAACCGGCGACGCTGGGGCAGAAGTAAAAGACAAAAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCt                                                                                                                                                                           <  2:645265/224‑1 (MQ=255)
                                                          cAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCTTCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCAAGAAAATGGCACCGCCGCAGATTTATGCGGAATTGCTTAAAAAAATGAAGAAAACCg                                                                                                                                                                             >  1:143244/1‑226 (MQ=255)
                                                              cGGCAAAACGTAAGGCAGTGACCACCCCGAAAACCACTCTGTTTGCGATTAACCGCCTGATTGGTCGCCGCTTCCAGGACGAAAAAGTACAGCGTGATGTTTCCATCATGCCGTTAAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCAAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGa                                                                                                                                                                      <  1:597946/229‑1 (MQ=255)
                                                                          aGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAAGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTCTTGCTGCTGATCACGGCGACGCCTGGGTCCAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGGGCTGCAAAAAATTAAGAAAACCGCTGAAGCTAACCTGGGTGa                                                                                                                                                          >  1:189713/1‑229 (MQ=255)
                                                                                tGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCCCCGCCGCAGATTTCTGCTGAGGGGCTGGAAAAAATGAAGAAAACCGCTGACGATTACCTGGGTGAACCggg                                                                                                                                                    >  1:346458/1‑228 (MQ=255)
                                                                                      aCCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACACCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCGGCAGATTTCTCCTGAAGGGCTGAAAAAAATGAAGAAAACAGCTGACAGTTACCTGGTGGAACCGGTAACTg                                                                                                                                               >  1:72725/1‑228 (MQ=255)
                                                                                          gCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGACGATTTCTGCTGAAGTGCTGGAAACACTTAAGAAAACCGCGGAAGATTAACTGGGTGAACCGGTAACTGaaac                                                                                                                                           >  1:331625/1‑226 (MQ=255)
                                                                                          gCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCAACGCCGCCGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAAGTTACCCGGGTGAACCCGGCACTTAAGca                                                                                                                                          >  1:347086/1‑228 (MQ=255)
                                                                                                  ctctGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATAAAGAAAACCGCTGAAGATTACCTGGGGGAACCGGAAACTAAAGCGGTTaaaac                                                                                                                                  >  2:119008/1‑225 (MQ=255)
                                                                                                   tctGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGGGACCCGGTAACTGAAGCTGTTAtaaac                                                                                                                                 >  1:510022/1‑225 (MQ=255)
                                                                                                          gCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCg                                                                                                                           >  1:623303/1‑228 (MQ=255)
                                                                                                                    gCCTGATGGGCCCCCGTTTCAGGGACGAAGAAGTAGAGGGGGATTTTTCCATCATGCCGTTAAAAATTATTGCTGCTGATAACGGCGACGTAGGGGTAGAAGTAAAGGGCAAGAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTaa                                                                                                                <  1:119008/229‑1 (MQ=255)
                                                                                                                      cTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAACAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAACCTGTTATCACCGTACCAGCAAACTTTaaa                                                                                                               >  2:259967/1‑227 (MQ=255)
                                                                                                                           tggcccccgtTTCAAGGAGCAAAAAATACAGGGTGATTTTTCCATCATGCCTTACAAATTTATTGCTGTTGAAAAGGGCGCGGAAGGGGCCAAAGTAAAAGGCAAAAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCt                                                                                                         <  1:150810/223‑1 (MQ=255)
                                                                                                                           gggTCGCCGTTTCCAGGAAGAAAGAATACAGGGTGATGTTTCAAATGTGCCTTTCAAAATTATTGAGGAGGATACGGGCGACGCGTGGGTCAAAGTTAAAGGCAAAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCt                                                                                                         <  1:501073/228‑1 (MQ=255)
                                                                                                                                gccgcTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCCGATTTCTGCTGAAGGGCTGAAAAAAATGAAGAAAACCGCTGAAGCTTACCTGGGTGAAACGGTAACTGAAGCGGTTAAAACCGAACCGGCATCCTTTAACGAttttcaagc                                                                                                    >  2:407237/1‑221 (MQ=255)
                                                                                                                                     ttCAAGGAGAAAGAAGCAGAGGGGTATGTTTCCATCGTGCCTTTAAAAATTATTGCCGGTGAAGACGGGGAGGGAGGGGCAGATGTAAAGGCCAAAAAAAGGGACCCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAgg                                                                                               <  2:49290/229‑1 (MQ=255)
                                                                                                                                     gtcgaGGACGAAAAAGGACAGCGTTATGTTCCGTACAGGCCGGAAAAATTTATGGCGGTGATAAACGGCGACGAAGGGTGCGAAGTAAAAGGCAAGAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAgg                                                                                               <  1:440333/225‑1 (MQ=255)
                                                                                                                                            aCGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGGGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCa                                                                                                                        >  1:57455/1‑197 (MQ=255)
                                                                                                                                                     tACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATGTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTTGGTGAACCGGTAACGGAAGCTGTTATCACCGTACCGGCATACTTTAACGATTCTCAGCGTCAGGAAACCAAAGACACAgg                                                                               >  2:743912/1‑229 (MQ=255)
                                                                                                                                                      aCAGCGTGATGTTTCCATCATGCCGTAAAAAATTATTGCTGCTGAAAACGGCGCCGAATGGGTAAAATTAAAAGCCAAAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTGCCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCa                                                                                 <  1:86078/226‑1 (MQ=255)
                                                                                                                                                       cAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGCAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGACAATTACCTGGGTGAACCGGTACATGAAGCTGTAATCACCGTAACGGCATACTTTAACGCATCTCAACGCATGGCAACAAAAAACGCAGGc                                                                              >  2:74146/1‑228 (MQ=255)
                                                                                                                                                         gCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATAACCGTACCGGCATACTTTAACGATGCCCAGCGTCAGTCAACCAAAGACGCAGGcc                                                                             >  1:420252/1‑227 (MQ=255)
                                                                                                                                                          cGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGCAAAACATGAAGAAAACCGCTGAAGATTACCTCGGTGAACCGGGAACTGAAGCCGTTACCAACGTACCGGCATACTTTAACCATGCTCAGCGTCACGCACCCAAAACCGCAGccc                                                                             >  2:299299/1‑226 (MQ=255)
                                                                                                                                                          cGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCt                                                                                                                                                                           >  1:748810/1‑132 (MQ=255)
                                                                                                                                                          cGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGcc                                                                             >  1:387407/1‑226 (MQ=255)
                                                                                                                                                          cGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACATGGTTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCAGACTTTAACGAGGCTCAGCGTCCGGCAACCCAAGACTCAGGccatc                                                                          >  1:611579/1‑226 (MQ=255)
                                                                                                                                                            tGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCa                                                                                                                        >  1:513023/1‑181 (MQ=255)
                                                                                                                                                                   tCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGAGGCTCAGAGTCAGGCAACAAAAAACACAGGCCGTAACCCTggac                                                                 >  2:387884/1‑227 (MQ=255)
                                                                                                                                                                   cccATCATGCCTTACAAAATTATGGCTGGTGATAACGGCTGGGATGGGATGGAATTTAAAGGGAAGAAATTGGTACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTc                                                                 <  1:712406/228‑1 (MQ=255)
                                                                                                                                                                       tcatGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGa                                                             >  1:616846/1‑229 (MQ=255)
                                                                                                                                                                           gCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGGAACTGAAGCTGTTCTAACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGGTTGGAAGTa                                                         >  1:739961/1‑229 (MQ=255)
                                                                                                                                                                               ttCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGGGAACCGGTAACTGAAGCTGTTATCACCGAACCGCCATCATTTAACGATCCTAAGCGTCAGGCAACCAAAGACGCAGGCCGCATGGCTGGTCTGGAAGTaaaa                                                      >  1:108594/1‑228 (MQ=255)
                                                                                                                                                                                 cAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGt                                                   >  1:370212/1‑229 (MQ=255)
                                                                                                                                                                                                              cGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTggcc                                                                 >  1:384470/1‑186 (MQ=255)
                                                                                                                                                                                                                     ggttcagagtAAAGGCCAAAAAACGCCCCCCGCCGCAGATTTCAGCTGAAGTGCGTAAAAAAATGAAGAAAACCTCTGAAGATTACCGGGTGGAACCGGTAACGGAAGCCGTTATCCCCGTACCGGAATACTTTAACTATGCTCAGCCTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCt                      <  2:580881/217‑1 (MQ=255)
                                                                                                                                                                                                                       ccaaaGTAAAAGGCAAAAAAATGCCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCtt                 <  1:570656/225‑1 (MQ=255)
                                                                                                                                                                                                                        caaatttAAAGGCAAAAAAGTGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTa                <  1:517683/223‑1 (MQ=255)
                                                                                                                                                                                                                                aaGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGGATCATCAACGAACCAACCGCAACTGCGCTGGCTTAAGGACTAGa       >  2:156873/1‑229 (MQ=255)
                                                                                                                                                                                                                                     cAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGCACGCAGCTGCGCTGGCTTACGGACGGGAAAAAg  >  2:737516/1‑229 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                      
CGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAG  >  NC_000913/12258‑12712

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