Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 512,643 A→G 57.4% 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.4% 35.0 / 81.3 75intergenic (+711/‑150)recR/rrs1recombination protein RecR/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (18/14);  new base G (19/24);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

TAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCAT  >  CP000730/512508‑512781
                                                                                                                                        |                                                                                                                                          
tAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACATGAATTACTAACGACAATTAATAATAAAAAGCAATTGACTTAAATGTTAATaaa                                                                                                                                       >  1:148849/1‑141 (MQ=255)
gaagttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcaggtttttaatcaaattttggtgaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTtgataa                                                                                                                                        <  1:254418‑M1/28‑1 (MQ=255)
 aagttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATaaa                                                                                                                                       <  2:71911‑M1/29‑1 (MQ=255)
  agttaagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATaaaa                                                                                                                                      <  2:213460‑M1/30‑1 (MQ=255)
      aagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaattatattttggtgaaaaaataacatggacaagataaAAAAAGATATTGACTTAAATGTTGCTAAAATGta                                                                                                                                  <  2:58631‑M1/34‑1 (MQ=255)
       tACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGtata                                                                                                                                >  2:269585/1‑141 (MQ=255)
        gctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGtata                                                                                                                                <  1:236468‑M1/36‑1 (MQ=255)
          tgttgtTAAGCAGTTTTAAAGCCTTTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACGTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaatt                                                                                                                             <  2:30073/141‑1 (MQ=255)
           ccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATaaaa                                                                                                                                      <  1:87294‑M1/30‑1 (MQ=255)
           ccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATaaaa                                                                                                                                      >  2:87294‑M1/103‑132 (MQ=255)
            cttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaatta                                                                                                                            >  1:326679‑M1/102‑141 (MQ=255)
                   gCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGt                                                                                                                    >  1:263433/1‑141 (MQ=255)
                     aGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCg                                                                                                                  >  2:105872/1‑141 (MQ=255)
                          gAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTCAATAAAATGTATAATTAATTCTTGTCGGTAAg                                                                                                             >  2:131654/1‑141 (MQ=255)
                            gtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGa                                                                                                            <  2:322515‑M1/56‑1 (MQ=255)
                             tcgaacttacgttccgctagagtataacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGaa                                                                                                           <  2:60170‑M1/57‑1 (MQ=255)
                                ttttataTTATTTATTTGTTGAGAAAAAAAAAACAAATCAACATAAAAAATAACTTGAAATACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTATAAAAAAGAAAAAt                                                                                                       <  2:43241/138‑1 (MQ=255)
                                   tataGTATTTATTTGTTGAGGCAAAGAAAACAAATCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTATTGTAGATAAGAAAAATGaa                                                                                                    >  1:139966/1‑141 (MQ=255)
                                     acgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAAc                                                                                                   <  1:173918‑M1/65‑1 (MQ=255)
                                      cgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACa                                                                                                  <  2:23841‑M1/66‑1 (MQ=255)
                                       gttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAt                                                                                                 >  1:249922‑M1/75‑141 (MQ=255)
                                       gTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTATTGTCGGTAAGAAAAATGAACAt                                                                                                 >  2:313842/1‑141 (MQ=255)
                                       gTAGTTATTTTTTAAGGTAAACAAAACAACTAAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTAACTTAAATGTTAATAAAATGTATAATTAATTCTTGTAGATAAGAAAAATGAACAtt                                                                                                <  2:163471/141‑1 (MQ=255)
                                        tatttattTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg                                                                                               <  1:151120/141‑1 (MQ=255)
                                        tatttattTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg                                                                                               <  2:307815/141‑1 (MQ=255)
                                                    gAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGa                                                                                   >  1:309347/1‑141 (MQ=255)
                                                     gaactttgccaggcagttttgtaatcaaattttggtaaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCTGTAAGAAAAATGAACATTAAAAACTAAATGa                                                                                   <  2:270365‑M1/81‑1 (MQ=255)
                                                     gaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGa                                                                                   >  1:306878‑M1/61‑141 (MQ=255)
                                                     gaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGa                                                                                   <  1:44221‑M1/81‑1 (MQ=255)
                                                       acgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACa                                                                                 <  1:75374‑M1/83‑1 (MQ=255)
                                                        agttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAAAAAAATGAACATTGAAAACTGAATGACaa                                                                                >  2:229297‑M1/58‑141 (MQ=255)
                                                          aacaaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAatat                                                                             >  1:115328/1‑141 (MQ=255)
                                                            gccaggcagtttttttatcaaattttggttaaaaagtaagatggagaggataaAAAGGGTTATTGACTTAAATGTGGATAAAATGTGTAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACaatatt                                                                            <  2:249240‑M1/88‑2 (MQ=255)
                                                               aggcagttttttaatcaaattttggttaaaaaataaaatgggcaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa                                                                         >  1:148639‑M1/51‑141 (MQ=255)
                                                                  cagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACg                                                                      <  1:243533‑M1/94‑1 (MQ=255)
                                                                  aactcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGt                                                                     >  1:49958/1‑141 (MQ=255)
                                                                  aactcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGt                                                                     >  1:241300/1‑141 (MQ=255)
                                                                    ataaactTAAGAAATAAATTGAATTACTAACGAAAAATAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTAAAAAATGAATAAAAATATGTaaaaaaaa                                                                   <  2:109809/138‑6 (MQ=255)
                                                                        tttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAAtt                                                                >  1:213460‑M1/42‑141 (MQ=255)
                                                                              gAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaa                                                         >  1:60927/1‑141 (MQ=255)
                                                                              caaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaa                                                          <  1:17507‑M1/106‑1 (MQ=255)
                                                                              caaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaa                                                          <  1:122913‑M1/106‑1 (MQ=255)
                                                                                 aTAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACg                                                      <  1:86547/141‑1 (MQ=255)
                                                                                   aaCTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAACGTaa                                                   >  1:75005/1‑141 (MQ=255)
                                                                                    ttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTa                                                    <  2:72912‑M1/112‑1 (MQ=255)
                                                                                      ggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAACTGAACATTGAAAACTGAATGACAATATGTCAATGTTATTTACAAAAAACGTAAc                                                  <  2:135804‑M1/114‑1 (MQ=255)
                                                                                          aTTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATaa                                             >  1:284760/1‑141 (MQ=255)
                                                                                           ttACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAg                                            >  1:22280/1‑141 (MQ=255)
                                                                                            tACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGt                                           <  1:311274/141‑1 (MQ=255)
                                                                                              atAACGAAAATTAATTTTAAAAATTTATTTACTTAAATGTTAATAAGATGTATAATGAATTCTTGTAGGTAAGAAAAATTAACATTGAAAACTGAATGACAACATGTCAACGTTAATTCCAAAAAACGTAACttatagtta                                         <  2:229035/140‑9 (MQ=25)
                                                                                               ataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAaacgtaacttatagtta                                         >  2:89333‑M1/19‑133 (MQ=255)
                                                                                                 aaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa                                       >  1:60170‑M1/17‑141 (MQ=255)
                                                                                                   aatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTAGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATACCAAAAAACGTAACTATAAGTTACaaa                                     >  2:112388‑M1/15‑141 (MQ=255)
                                                                                                    atggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAAc                                    <  1:214275‑M1/128‑1 (MQ=255)
                                                                                                        ttaattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACattct                               >  1:173322/1‑139 (MQ=255)
                                                                                                           attTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTa                            >  2:88006/1‑141 (MQ=255)
                                                                                                            gataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTa                            <  2:168604‑M1/136‑1 (MQ=255)
                                                                                                             ataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAg                           >  1:229735‑M1/5‑141 (MQ=255)
                                                                                                             ataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAg                           >  2:135598‑M1/5‑141 (MQ=255)
                                                                                                              taaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAGACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGt                          >  2:158954‑M1/4‑141 (MQ=255)
                                                                                                              taaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGt                          >  1:121886‑M1/4‑141 (MQ=255)
                                                                                                                     aGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATg                   >  2:210487/1‑141 (MQ=255)
                                                                                                                      gTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAAAGAACATTGAAAACTGAATGACAATATGTCAAAGTAAAATAAAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGa                  <  2:236468/141‑1 (MQ=21)
                                                                                                                       ttATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAg                 <  1:72912/141‑1 (MQ=255)
                                                                                                                       ttATTGAACTCAATGTTAATAAAATGTACAATAAATACTTGTCAGTCAGAAAAATGAAAATAGAACACTGAATCACAATATGAAAAAAAAAAAAAAAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAg                 <  2:76719/141‑1 (MQ=12)
                                                                                                                        tATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGc                >  1:191167/1‑141 (MQ=255)
                                                                                                                        tATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGc                >  1:145250/1‑141 (MQ=255)
                                                                                                                        tATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGc                >  2:73658/1‑141 (MQ=255)
                                                                                                                         aTTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCt               >  1:62775/1‑141 (MQ=255)
                                                                                                                         aTTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAACAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCt               <  2:46366/141‑1 (MQ=255)
                                                                                                                         aTTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCt               <  1:295333/141‑1 (MQ=255)
                                                                                                                            gACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAAt            <  1:322515/141‑1 (MQ=255)
                                                                                                                                 aaaTGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAc       >  1:193531/1‑141 (MQ=255)
                                                                                                                                 aaaTGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGGAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAc       <  1:16509/141‑1 (MQ=255)
                                                                                                                                 aaaTGGTGATAAAAATTATAATTAAGTCTTGTCGGTAAGAAAAATGAACGTGGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCaatc       <  2:61238/141‑3 (MQ=17)
                                                                                                                                      ttAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATACAAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcatcat  <  2:90138/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                          
TAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCAT  >  CP000730/512508‑512781

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