Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 512,643 A→G 57.8% intergenic (+711/‑150) recR → / → rrs1 recombination protein RecR/16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730512,6430AG57.8% 26.2 / 84.3 78intergenic (+711/‑150)recR/rrs1recombination protein RecR/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (18/15);  new base G (24/21);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-02

GTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACAT  >  CP000730/512507‑512778
                                                                                                                                        |                                                                                                                                       
gtAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaata                                                                                                                                      >  2:131862/1‑140 (MQ=255)
    tttttACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAt                                                                                                                                  <  1:264254/140‑1 (MQ=255)
      aagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGt                                                                                                                                <  2:103806‑M1/33‑1 (MQ=255)
        tACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGtat                                                                                                                              <  2:165768/140‑1 (MQ=255)
        tACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATCTTAATAAAATGtat                                                                                                                              <  1:186538/140‑1 (MQ=255)
         ctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGtata                                                                                                                             <  1:229849‑M1/36‑1 (MQ=255)
          cTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaa                                                                                                                            >  1:57837/1‑140 (MQ=255)
          cTGTTGTTAAGCAGGTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaa                                                                                                                            >  2:88588/1‑140 (MQ=255)
              tgttAAGCAGTTTGAAAACCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACATGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaattaa                                                                                                                        <  1:191256/140‑1 (MQ=255)
                 cgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTc                                                                                                                     <  1:212937‑M1/44‑1 (MQ=255)
                  cagcactgttcaaGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGCATTACTCACGTAATTTAATTTTAAAAATTTATTGACCTAAATGTTAATAAAATGTATACTTAATTCt                                                                                                                    <  2:21850/129‑1 (MQ=255)
                    gCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTg                                                                                                                  >  2:233616/1‑140 (MQ=255)
                    cgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTg                                                                                                                  <  1:126175‑M1/47‑1 (MQ=255)
                     gatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaatcaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAAttctggt                                                                                                                 <  2:80038‑M1/48‑4 (MQ=255)
                     gatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGAATTAAATGTTGATAAAATGTATAATTAAttctcgt                                                                                                                 <  2:71713‑M1/48‑4 (MQ=255)
                      aGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTc                                                                                                                >  1:146977/1‑140 (MQ=255)
                           agtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTaa                                                                                                           >  2:201462‑M1/87‑140 (MQ=255)
                            gtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacatgttaaAAAAAGTTATTGACTTAAATGTTGATAAAATGAATAATTAATTCTTGTCGGTAAg                                                                                                          <  1:101848‑M1/55‑1 (MQ=255)
                            gtcgaaattaagttccgttagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAg                                                                                                          >  1:162518‑M1/86‑140 (MQ=255)
                                aacttacgttccgctagagtagaacgttgccaggcagttttttaatcagattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAagaagg                                                                                                      >  2:148591‑M1/82‑138 (MQ=255)
                                aacttacgttccgcaagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGaaaa                                                                                                      >  1:291008‑M1/82‑140 (MQ=255)
                                 cTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTTTCGGTAAGaaaaa                                                                                                     >  2:304863/1‑140 (MQ=255)
                                 cTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGaaaaa                                                                                                     >  1:9344/1‑140 (MQ=255)
                                      cgttacgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAAGGAAc                                                                                                <  2:169678‑M1/65‑1 (MQ=255)
                                          ccgtttggatagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg                                                                                            >  2:46879‑M1/72‑140 (MQ=255)
                                          ccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg                                                                                            <  2:320207‑M1/69‑1 (MQ=255)
                                              atttGTTGAGGCAAACAAAACAACTCAACTTAAGACATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGACGGTAAGAAAAATGAACATTAaaaa                                                                                        >  1:35988/1‑140 (MQ=255)
                                              atttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa                                                                                        <  1:131862/140‑1 (MQ=255)
                                                     gaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagagaaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATg                                                                                 <  1:67889‑M1/80‑1 (MQ=255)
                                                        cgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTACATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGGAAAATGCACATTGAAAACTGAATGACa                                                                              >  2:235832‑M1/58‑140 (MQ=255)
                                                        agttgacaggcagttttttaattaaatttttgttaaaagatgaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTGTTGTAGGTAAGAAAAATGAAAATTGAAAACTGAATGACa                                                                              <  2:90871‑M1/83‑1 (MQ=255)
                                                         gttgccaggcagttttttaatcaaattttggttaaaaactgtacgggagaaattaaAAAAAGTTATTGACTTAAATGTTGATACAATGTATAATTAATTCGTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACaa                                                                             >  1:203230‑M1/57‑140 (MQ=255)
                                                         gttgccaggcagttttttaatcaaattttggttaaaaaataaaatggactagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACaa                                                                             >  1:273805‑M1/57‑140 (MQ=255)
                                                         gttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTAACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACaa                                                                             <  1:9545‑M1/84‑1 (MQ=255)
                                                            acacaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTAGTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAattt                                                                          >  1:256116/1‑138 (MQ=255)
                                                              caggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGt                                                                        <  2:326088‑M1/89‑1 (MQ=255)
                                                                 acaacTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGGAAAATGAACATTGAAAACTGAATGACAATATGTCaa                                                                     >  2:120064/1‑140 (MQ=255)
                                                                 acaacTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCaa                                                                     >  2:306488/1‑140 (MQ=255)
                                                                  cagtgttgatatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAc                                                                    >  2:88489‑M1/48‑140 (MQ=255)
                                                                    actcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGt                                                                  >  1:148354/1‑140 (MQ=255)
                                                                        tttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAt                                                              <  1:336892‑M1/99‑1 (MQ=255)
                                                                         ttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTAACTTAAATGATGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAAGATAAtt                                                             >  2:201750‑M1/41‑140 (MQ=255)
                                                                         actTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAACGTTATTGACTTACATGTTAATAAAATGTATAATTAATTAAAGTAGATAAGAAAAATGAACATTAAAAACTGAATCACAATATGTCAACCTTAAtt                                                             >  1:54607/1‑140 (MQ=255)
                                                                         actTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAtt                                                             >  1:281460/1‑140 (MQ=255)
                                                                           aatccaattttggttaaaaaataaaatggacacgataaAAAAATTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTcc                                                           >  1:206772‑M1/39‑140 (MQ=255)
                                                                            atcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCa                                                          >  2:126175‑M1/38‑140 (MQ=255)
                                                                             tcaaagtttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaa                                                         >  2:66246‑M1/37‑140 (MQ=255)
                                                                             aaGAAATAACTTGAATTACTAACTAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaa                                                         <  2:83101/140‑1 (MQ=255)
                                                                              caaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTGGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAAAATGTCAACATTAATTCCaaa                                                        <  1:127310‑M1/105‑1 (MQ=255)
                                                                                aaaTAACTTTAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaa                                                      <  1:715/140‑1 (MQ=255)
                                                                                 aaTAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaaa                                                     <  1:321894/140‑1 (MQ=255)
                                                                                    aaCTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGt                                                  >  1:137213/1‑140 (MQ=255)
                                                                                     aCTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAATAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTa                                                 <  2:289261/140‑1 (MQ=255)
                                                                                        ttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTAAACGTTAATTCCAAAAAACGTAACt                                              <  2:148749‑M1/115‑1 (MQ=255)
                                                                                          aaaaaataaaatggacaagataaAAAAATTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACtat                                            >  2:136967‑M1/24‑140 (MQ=255)
                                                                                          aaaaaaaaaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACtat                                            >  2:329584‑M1/24‑140 (MQ=255)
                                                                                           aaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACtata                                           >  1:166323‑M1/23‑140 (MQ=255)
                                                                                           aaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACtata                                           >  1:195517‑M1/23‑140 (MQ=255)
                                                                                             tACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACGGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAg                                         >  1:340338/1‑140 (MQ=255)
                                                                                              acgaaCGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGt                                        <  1:150344/137‑1 (MQ=255)
                                                                                              aCTAACGAAAATTAATTTTAAAAAGTTATTGGCGGAAATGTTAATGAGATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCACAACCCGTACCTCTACGt                                        <  2:161001/140‑1 (MQ=25)
                                                                                               ataaaatggaaaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGtt                                       >  2:19830‑M1/19‑140 (MQ=255)
                                                                                                 aaaatggacaagataaAAAAAGTTATTGAATTAAATGTTGATAAAATGGATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAAAGTAAATACCAAAAAACGTAACTATAAGTTAc                                     >  1:20925‑M1/17‑140 (MQ=255)
                                                                                                  aaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa                                    >  1:17745‑M1/16‑140 (MQ=255)
                                                                                                  aaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa                                    >  1:36264‑M1/16‑140 (MQ=255)
                                                                                                   aatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTATTGTCGGTAAGAAAAATGAAAATTGAAAACTGAATGAAATTATGTCTACGTAAATTCCAAAAACAGTAACTATAAgttaaaa                                   <  2:105204‑M1/126‑4 (MQ=255)
                                                                                                    gAAAATTAATTTTAAAAAGTTATTGACTTAAATTTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa                                  <  2:218907/140‑1 (MQ=255)
                                                                                                    atggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAGTCTGAACTATAAGTTACaaa                                  <  1:257110‑M1/127‑1 (MQ=255)
                                                                                                    atggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa                                  <  2:256050‑M1/127‑1 (MQ=255)
                                                                                                               aaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGt                       >  1:103806‑M1/3‑140 (MQ=255)
                                                                                                               aaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATAACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGt                       <  1:286400‑M1/138‑1 (MQ=255)
                                                                                                               aaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAAGTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGt                       >  1:19336‑M1/3‑140 (MQ=255)
                                                                                                                        tATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTCCAAAAATTATTACGTATTTCTGaa              >  2:138279/1‑139 (MQ=21)
                                                                                                                        tATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAg              <  1:79459/140‑1 (MQ=255)
                                                                                                                        tATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAg              >  2:316855/1‑140 (MQ=255)
                                                                                                                          ttGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCt            <  1:329467/140‑1 (MQ=255)
                                                                                                                            gACTTAAATGTTGATAAAATGTATAATTGATTCTTGTCGGTAAGAAAAATGAACGTTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACTAACAATATCTAGAATTTATGAGCTaa          >  1:230277/1‑140 (MQ=18)
                                                                                                                             tcTTAAATGTTGATAAAATGTATAATTTATACTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAAt         <  1:209011/139‑1 (MQ=255)
                                                                                                                                    tGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcat  >  2:305616/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                       
GTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACAT  >  CP000730/512507‑512778

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