Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 512,643 A→G 48.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,6430AG48.1% 5.1 / 105.5 77intergenic (+711/‑6052)recR/speA_2Recombination protein RecR/Arginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (24/16);  new base G (19/18);  total (43/34)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATA  >  USA300TCH1516_ALE/512504‑512782
                                                                                                                                           |                                                                                                                                            
ggTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTtaa                                                                                                                                             >  1:250888/1‑141 (MQ=255)
      aTTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAAt                                                                                                                                       >  1:202375/1‑141 (MQ=255)
        aaagctacgtagagtcaatggtagtcgaaattaagtaccgctagagtagaaagttgccaggaagtttaataatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGt                                                                                                                                     >  2:2781‑M1/109‑141 (MQ=255)
         tttACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGta                                                                                                                                    <  1:355481/141‑1 (MQ=255)
         tttACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGta                                                                                                                                    <  1:32198/141‑1 (MQ=255)
              tcttagcgtcgatagtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaatt                                                                                                                               >  2:412185‑M1/103‑141 (MQ=255)
               gttgttAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGACATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTAtaatta                                                                                                                              <  1:229604/141‑1 (MQ=255)
               cttagcgtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTAtaatta                                                                                                                              <  1:206389‑M1/40‑1 (MQ=255)
                    tAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCt                                                                                                                         <  1:399884/141‑1 (MQ=255)
                     gtcgatggtagtcgaacttacgttcctctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCtt                                                                                                                        <  1:83979‑M1/46‑1 (MQ=255)
                     gtcgatggtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCtt                                                                                                                        <  2:121759‑M1/46‑1 (MQ=255)
                      aGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTg                                                                                                                       >  1:235668/1‑141 (MQ=255)
                           tttGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAGTTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGt                                                                                                                  <  1:180295/141‑1 (MQ=255)
                           tttGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGt                                                                                                                  >  1:374386/1‑141 (MQ=255)
                            gtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTa                                                                                                                 <  1:317034‑M1/53‑1 (MQ=255)
                            gtagtcgaacttacgttccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTa                                                                                                                 <  1:405677‑M1/53‑1 (MQ=255)
                               aaaGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGa                                                                                                              >  2:347388/1‑141 (MQ=255)
                                    cTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAACTGTTACTAAAATGTATAATTTATTCTTGTCGGTAAGAAAAAt                                                                                                         >  1:26626/1‑141 (MQ=255)
                                            tccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg                                                                                                 >  1:278657‑M1/73‑141 (MQ=255)
                                            tccgctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg                                                                                                 >  1:137253‑M1/73‑141 (MQ=255)
                                            tccgagagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTg                                                                                                 >  1:125795‑M1/73‑141 (MQ=255)
                                                tatttGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGaaaa                                                                                             >  2:380547/1‑141 (MQ=255)
                                                ctagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACACTGaaaa                                                                                             >  2:252519‑M1/69‑141 (MQ=255)
                                                 tagagtagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAAc                                                                                            >  1:391432‑M1/68‑141 (MQ=255)
                                                   ttgttgAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTg                                                                                          <  2:13542/141‑1 (MQ=255)
                                                      tagaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAAt                                                                                       >  1:298483‑M1/63‑141 (MQ=255)
                                                        gaacgttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTCGTCGGTAAGAAAAATGAACATTACAAActgaacga                                                                                     <  2:34640‑M1/81‑4 (MQ=255)
                                                        gaacgttgccaggcagttttttaatcaaattttggttaaaaaagaaaatcgacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTCTCGGAAAGAAAAATTAACATTGAAaaatgaatgc                                                                                     <  2:225726‑M1/81‑2 (MQ=255)
                                                          ggCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATCTGAAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑Ca                                                                                  >  2:90851/1‑141 (MQ=255)
                                                           gCAAACAAAACAACTCAACTCAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑Caa                                                                                 <  1:12624/141‑1 (MQ=255)
                                                            gttgccaggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAat                                                                                <  1:354205‑M1/85‑1 (MQ=255)
                                                                caaaacaaCTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAACATGAACATTGAAAACTGAATGA‑CAATATGt                                                                            >  1:218934/1‑141 (MQ=255)
                                                                 caggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTc                                                                           <  2:312027‑M1/90‑1 (MQ=255)
                                                                 caggcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTc                                                                           <  1:257957‑M1/90‑1 (MQ=255)
                                                                    gcagttttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAAc                                                                        <  1:166989‑M1/93‑1 (MQ=255)
                                                                      aactcaactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGt                                                                      >  1:398865/1‑141 (MQ=255)
                                                                      aaatcaactTAAGAAATAACTTGAATTAATAAAAAAGAATGGTTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGt                                                                      <  2:4126/138‑1 (MQ=38)
                                                                          ttttaatcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGCATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAATTGAATGA‑CAATATGTAAACGTTAAt                                                                  <  2:406340‑M1/99‑1 (MQ=255)
                                                                          caactTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAAt                                                                  >  1:284633/1‑141 (MQ=255)
                                                                           aactTAAGAAATAACCTGAATAGCTAACGAAAAGTAATTTTAAAAAGCTATTGACTTAGATGTTAATAAAATGTATAATTAATACTTGTCGGTAAGAAAAATGAACATTGAAAACACAATGA‑CAATGTGTCAACGATagtg                                                                 >  1:114469/1‑138 (MQ=255)
                                                                            actTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTc                                                                >  1:278308/1‑141 (MQ=255)
                                                                            actTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTc                                                                >  1:325606/1‑141 (MQ=255)
                                                                            actTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTc                                                                >  2:230180/1‑141 (MQ=255)
                                                                               tAAGAAATAACGTGAATTAATAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCaa                                                             <  2:253784/141‑1 (MQ=255)
                                                                               atcaaattttggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTAAACGTTAATTCCaa                                                             <  2:162351‑M1/104‑1 (MQ=255)
                                                                               ataaaattttgaataaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCaa                                                             >  2:257079‑M1/38‑141 (MQ=255)
                                                                                 aGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCaaaa                                                           >  1:195184/1‑141 (MQ=255)
                                                                                    aaTAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACATTAATTCCAAAAAAc                                                        >  2:360278/1‑141 (MQ=255)
                                                                                       ttagttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCAAAAAACGTa                                                     >  2:51013‑M1/30‑141 (MQ=255)
                                                                                         ggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCAAAAAACGTAAc                                                   >  2:330639‑M1/28‑141 (MQ=255)
                                                                                         ggttaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACATTAATTCCAAAAAACGtaaa                                                   <  2:201401‑M1/114‑2 (MQ=255)
                                                                                         cTTGAATTACTAACGAAAATTAATTTTAAAATGTTATTGACTGAAATGTCAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAATACTGCATGA‑CAATATGTCAACGTTCATTCCaaaaaaaataaa                                                   >  2:198110/1‑135 (MQ=255)
                                                                                          gaaaaaaaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCAAAAAACGTAACt                                                  >  2:368077‑M1/27‑141 (MQ=255)
                                                                                                aaataaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGt                                            >  2:405677‑M1/21‑141 (MQ=255)
                                                                                                 aCTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGtt                                           <  1:294501/141‑1 (MQ=255)
                                                                                                  atAACGAAAATTAATTATAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTa                                          <  2:223806/140‑1 (MQ=255)
                                                                                                   taaaatggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTAc                                         <  1:190643‑M1/124‑1 (MQ=255)
                                                                                                       gAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAAc                                     <  2:235668/141‑1 (MQ=255)
                                                                                                        tggacaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACa                                    <  2:398157‑M1/129‑1 (MQ=255)
                                                                                                          aattaattTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACatt                                  >  1:254848/1‑141 (MQ=255)
                                                                                                          aattaattATAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACatt                                  <  1:95427/141‑1 (MQ=255)
                                                                                                           acaagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACatta                                 <  1:418808‑M1/132‑1 (MQ=255)
                                                                                                           acaagataaAAAAAGTTATTGACGTAAATGTTGATAAAATGTATAATTAATTCTTGTCAGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACatta                                 >  2:218175‑M1/10‑141 (MQ=255)
                                                                                                            caagataaAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACattat                                >  1:153856‑M1/9‑141 (MQ=255)
                                                                                                               attTTAAAAAGTTATTGACTTAAGTGTTAATAAAATGTATAATTAATTCTTGTCGGAAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATCATTTa                             >  2:114845/1‑141 (MQ=255)
                                                                                                                   tAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAACTGAACATTGAAAACTGGATGATCTCGAAG‑CAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAt                         >  1:81851/1‑141 (MQ=25)
                                                                                                                   tAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCC‑AAAAACGTAACTATAAGTTACAAACATTATTTAGTAtt                        <  2:111083/141‑1 (MQ=255)
                                                                                                                   aAAAAAGTTATTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATGGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAt                         <  1:419213‑M1/140‑1 (MQ=255)
                                                                                                                   aAAAAAGTTATTGAATAAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTAt                         >  2:314525‑M1/2‑141 (MQ=255)
                                                                                                                         gTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGa                   >  1:356414/1‑141 (MQ=255)
                                                                                                                         gTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCAAAAAAACGTAACTATAAATTAAAAAAATTCTTTAttatttatta                   >  2:35804/1‑139 (MQ=25)
                                                                                                                         gTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGCAAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGa                   >  1:81546/1‑141 (MQ=255)
                                                                                                                          ttATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTTAACGTTAAGTCCAAAAAGCGTAACTATAAGTTACAAACATTATTTAGTATTTATGAg                  <  1:327693/141‑1 (MQ=255)
                                                                                                                          ttATTGAATTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAg                  <  2:250888/141‑1 (MQ=255)
                                                                                                                            aTTGACTTAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCt                >  2:122566/1‑141 (MQ=255)
                                                                                                                               gACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAAt             >  1:39784/1‑141 (MQ=255)
                                                                                                                                   tAAATGTTGATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCaaa         >  1:63196/1‑141 (MQ=255)
                                                                                                                                      aTGTTGATAAAATGTATAATTAATTCTTGTCGATAAGAAAAATGAACATTGAAAACTAAATGA‑CAATATGTCAGCGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAAcat      >  2:379047/1‑141 (MQ=255)
                                                                                                                                          taataaAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATa  <  2:230348/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                            
GGTGTAATTTTTACTGTTGTTAAGCAGTTTGAAAGCCTGTATAGTATTTATTTGTTGAGGCAAACAAAACAACTCAACTTAAGAAATAACTTGAATTACTAACGAAAATTAATTTTAAAAAGTTATTGACTTAAATGTTAATAAAATGTATAATTAATTCTTGTCGGTAAGAAAAATGAACATTGAAAACTGAATGA‑CAATATGTCAACGTTAATTCCAAAAAACGTAACTATAAGTTACAAACATTATTTAGTATTTATGAGCTAATCAAACATCATA  >  USA300TCH1516_ALE/512504‑512782

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: ‑