Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 512,643 A→G 60.1% intergenic (+711/‑6052) recR → / → speA_2 Recombination protein RecR/Arginine decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE512,6430AG60.1% 48.1 / 82.6 81intergenic (+711/‑6052)recR/speA_2Recombination protein RecR/Arginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (18/14);  new base G (30/18);  total (48/33)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

ATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATA  >  USA300TCH1516_ALE/512510‑512782
                                                                                                                                     |                                                                                                                                            
gttaagctccttagcgtcgatggtagtcgaacttacgttccggtagaagagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATaaaa                                                                                                                                        >  1:372240‑M1/111‑140 (MQ=255)
   aagctccttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGt                                                                                                                                     >  1:519505‑M1/108‑140 (MQ=255)
          gtagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggataagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaatta                                                                                                                              >  1:490458‑M1/101‑140 (MQ=255)
           tagcatcgatggtagtcgaacttacgttccgctagagtaaaacgtagccaggcagttttttaatcaaattttgattaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaattaa                                                                                                                             >  2:155725‑M1/100‑140 (MQ=255)
            cgcgtccatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaattaat                                                                                                                            >  1:219850‑M1/99‑140 (MQ=255)
                     ggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCgg                                                                                                                   <  2:442056‑M1/51‑1 (MQ=255)
                      ttGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGt                                                                                                                  >  1:403036/1‑140 (MQ=255)
                      gtagtcgaaattacgttccgcaagagtagaacgaagccaggcaattatttaatcaaattttgattaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGt                                                                                                                  >  2:206874‑M1/89‑140 (MQ=255)
                         gtcgaacttccgttccgctagagtagaacgttgccaggcagttttttaaccaaattttgcgttaaaaataaaatggacaagataaAAAAAGTTCTTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAg                                                                                                               >  1:326474‑M1/86‑140 (MQ=255)
                         gtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAg                                                                                                               <  2:52922‑M1/55‑1 (MQ=255)
                         aaaGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATTTATAATTAATTCTTGACGGTAAg                                                                                                               >  2:2758/1‑140 (MQ=255)
                         aaaGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAg                                                                                                               <  1:180695/140‑1 (MQ=255)
                              cTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAAt                                                                                                         >  1:476410/1‑140 (MQ=255)
                                  acgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGaa                                                                                                      <  1:433910‑M1/64‑1 (MQ=255)
                                  acgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATAATTTCTTCACGGTAAGAAAAATGaa                                                                                                      <  2:45697‑M1/64‑1 (MQ=255)
                                  acggtccgctagagtagaacgttgcccggccgctttttaatcaaattttggttaaaaattaaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGaa                                                                                                      <  1:332054‑M1/64‑1 (MQ=255)
                                   taGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATATTAATAAAATGTATAATTACATCTTGTAGGTAAGAAAAATAAAc                                                                                                     >  2:444971/1‑140 (MQ=255)
                                   cgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagatacAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGGAAGAAAAATGAAc                                                                                                     <  1:60577‑M1/65‑1 (MQ=255)
                                    aGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACa                                                                                                    >  1:134404/1‑140 (MQ=255)
                                      tccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAtt                                                                                                  >  1:428966‑M1/73‑140 (MQ=255)
                                      tatttattTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACAtt                                                                                                  >  2:76797/1‑140 (MQ=255)
                                             gagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACt                                                                                           >  1:374779‑M1/66‑140 (MQ=255)
                                              tgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTTTTGGCTGAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTg                                                                                          >  1:417934/1‑140 (MQ=255)
                                               gttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGa                                                                                         >  1:557546/1‑140 (MQ=255)
                                                tagaacgtttccaggcagttttttaatcaaattttggttaaaaaataaaatggacaatataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGaa                                                                                        >  1:256073‑M1/63‑140 (MQ=255)
                                                tagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGaa                                                                                        >  2:537309‑M1/63‑140 (MQ=255)
                                                 tgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAGCATTGAAAACTGAAt                                                                                       >  1:274328/1‑140 (MQ=255)
                                                    acgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGAc                                                                                    >  2:399703‑M1/59‑140 (MQ=255)
                                                     gCAAACAAAACAACTCAACTTAAGAAATACCTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACa                                                                                   >  2:16989/1‑140 (MQ=255)
                                                          acatgcagttttttaaaaaaaatttgattaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATg                                                                              >  2:533956‑M1/53‑140 (MQ=255)
                                                            aggcagttttgtaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTc                                                                            >  1:320764‑M1/51‑140 (MQ=255)
                                                             tgcagttttttaatcaaattttagttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa                                                                          >  2:187984‑M1/50‑140 (MQ=255)
                                                             ggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCa                                                                          >  1:77153‑M1/50‑140 (MQ=255)
                                                              gcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTGTTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCaa                                                                         <  1:358817‑M1/92‑1 (MQ=255)
                                                              acaacTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCAGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTaaa                                                                         >  2:356616/1‑140 (MQ=255)
                                                                aactcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACg                                                                       >  1:95949/1‑140 (MQ=255)
                                                                 gttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAACCTGACCCTTGAAAACTTAATGACAATATGTCAACGt                                                                      >  1:244988‑M1/46‑140 (MQ=255)
                                                                 gttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGt                                                                      >  1:183792‑M1/46‑140 (MQ=255)
                                                                    ttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTaa                                                                   <  2:445411‑M1/98‑1 (MQ=255)
                                                                          aaGAAAAAATTATAATTACTAACGAAAATTAATTATAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaa                                                             <  2:214852/140‑1 (MQ=38)
                                                                             aatttttgttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaa                                                          >  1:406971‑M1/34‑140 (MQ=255)
                                                                             aattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaa                                                          >  1:200674‑M1/34‑140 (MQ=255)
                                                                              attttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaaa                                                         >  2:445732‑M1/33‑140 (MQ=255)
                                                                              attttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCaaaaaa                                                         >  1:141360‑M1/33‑140 (MQ=255)
                                                                                tttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACg                                                       <  2:358032‑M1/110‑1 (MQ=255)
                                                                                tttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATCAACATTGAAAACTGAATGACAATCTGTCAACGATACTTCCGACAAACg                                                       <  1:204844‑M1/110‑1 (MQ=255)
                                                                                 ttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGt                                                      >  1:158950‑M1/30‑140 (MQ=255)
                                                                                   ggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTaa                                                    <  1:561502‑M1/113‑1 (MQ=255)
                                                                                    ttGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc                                                   >  1:340992/1‑140 (MQ=255)
                                                                                    gtGAATAACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTCATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAAc                                                   <  1:28236/139‑1 (MQ=255)
                                                                                     ttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACt                                                  <  1:568906‑M1/115‑1 (MQ=255)
                                                                                     aaaaTTAATAACGAAAATTAATTTGAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGAAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACt                                                  <  2:268152/138‑1 (MQ=255)
                                                                                        aTTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACtata                                               <  2:134404/140‑1 (MQ=255)
                                                                                         ttgcTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAAGGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATaa                                              <  2:147958/137‑1 (MQ=255)
                                                                                         ttACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATaa                                              >  1:134752/1‑140 (MQ=255)
                                                                                           aCTAACGAAAATTAATTTTAAAAAGTTATTTACTAAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGt                                            <  2:345744/140‑1 (MQ=255)
                                                                                            ataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGtt                                           >  1:50516‑M1/19‑140 (MQ=255)
                                                                                            aaaaaatggaaaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGtt                                           >  1:354678‑M1/19‑140 (MQ=255)
                                                                                             taaaatggacaagagaaAAAAAGTTATTGACTGAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTa                                          >  2:426037‑M1/18‑140 (MQ=255)
                                                                                             tAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTAAACGTTAATTCCAAAAAACGTAACTATAAGTTa                                          >  2:525866/1‑140 (MQ=255)
                                                                                              aaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTAc                                         <  2:451202‑M1/124‑1 (MQ=255)
                                                                                              aaaaaggacaagctaaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTAc                                         >  1:60362‑M1/17‑140 (MQ=255)
                                                                                               aaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATCACAATATGTCAACGTTAATGCCGAAAAACGGAACTCTAACTTACa                                        <  2:274530‑M1/125‑1 (MQ=255)
                                                                                               aCGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCACCGTTAATTCCAAAAAACGAAACTATAAGTTACa                                        >  1:152589/1‑140 (MQ=255)
                                                                                               aCGAAAAATAAATTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACa                                        <  2:166979/140‑1 (MQ=255)
                                                                                                cGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTGCGCTTCCTTCCCAGGGAAAAATGACTGTGGTATTTCTGTCTTTAATTCCAAAAAACGTAACTATAAGTTACaa                                       <  1:570284/140‑1 (MQ=255)
                                                                                                aatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaa                                       >  1:52922‑M1/15‑140 (MQ=255)
                                                                                                 gAAAATTAATTTTAAAAAGTTAGTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATAGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACaaa                                      <  1:541356/140‑1 (MQ=255)
                                                                                                       taattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAATTTACAAACattat                                >  2:431041/1‑140 (MQ=255)
                                                                                                         gataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTAttt                              <  1:149098‑M1/135‑1 (MQ=255)
                                                                                                           tatAAAAAGTTATTGAATTAAATGTTAAAAAAATGTATAATTAATTCTTGTCAGTAAGAAAAATGAACAAAGAAAACTGAATGACAATATGTCAAAAAAAAAAAAAAAAAACGTAACTATAAGTTACAAACATTATTTAg                            <  2:484650/138‑1 (MQ=21)
                                                                                                            ttAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGt                           <  2:225710/140‑1 (MQ=255)
                                                                                                            aaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGt                           <  1:235873‑M1/138‑1 (MQ=255)
                                                                                                            aaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATGCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTagc                           <  2:521024‑M1/138‑2 (MQ=255)
                                                                                                               aaaaGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAtt                        >  2:390117/1‑140 (MQ=255)
                                                                                                                 aaGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCAAAAACACGTACCTATCCGCTACAAACATACTATAGTATCTa                      >  1:251113/1‑140 (MQ=25)
                                                                                                                  aGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTAt                     <  2:425025/140‑1 (MQ=255)
                                                                                                                        tGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTa               <  1:415379/140‑1 (MQ=255)
                                                                                                                           cTTAAATGTAAATAAAATGTATAATTAATTCTTGTCGCTCAGGCAGCGGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATc            <  2:553499/140‑1 (MQ=17)
                                                                                                                           ataaaaTGTTGAAAAAACGGATAATTAAGTCTTGTCACTAAGAAAAATAAACATAAAAAACTGAACGATAATATGTCAAAAAAAAAAAAAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATc            <  2:490050/137‑1 (MQ=12)
                                                                                                                             caaaTGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATCTCAACGTTAATTCCCAAAAACGTAACTACAAGCTACAAACATTATTTACTACTCATGAGCTAATCaa          >  2:547837/2‑140 (MQ=18)
                                                                                                                                     aataaAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATa  >  2:296519/1‑140 (MQ=255)
                                                                                                                                     aataaAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATa  >  1:329286/1‑140 (MQ=255)
                                                                                                                                     aataaAATGTATAAGTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATa  <  2:332449/140‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                            
ATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGACAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATA  >  USA300TCH1516_ALE/512510‑512782

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