Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 842,690 A→G 7.4% intergenic (+76/‑586) uvrA → / → hprK UvrABC system protein A/HPr kinase/phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE842,6900AG 7.4% 175.7 / 4.1 65intergenic (+76/‑586)uvrA/hprKUvrABC system protein A/HPr kinase/phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (29/31);  new base G (3/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

ACAGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGATAAATAAGATTAAAAGAAGTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACTAATGATTTATTATGTAGTGGTTCTTTGTCATTAGCCACAGCT  >  USA300TCH1516_ALE/842555‑842829
                                                                                                                                       |                                                                                                                                           
aCAGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAAAAATTTCCGGAAAACAAATAAGGTTCAAAGAAGTGAAGGATGTTATAAATTTATCCTTCGCTTCTGTTTATTAATTTAGTAATGAATAATAGaaagaaaag                                                                                                                                         >  2:782339/1‑140 (MQ=255)
      aactagttAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGATAAATAAGATTAAAAGAAGTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGt                                                                                                                                   <  2:378914/134‑1 (MQ=255)
      aaGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGATAAATAAGATTAAAAGAAGTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGt                                                                                                                                   <  1:741598/140‑1 (MQ=255)
        gTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGATAAATAAGATTAAAAGAAGTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAAGTTAGTAATGAATAGTAGAAAGAAAAGATGCGTaa                                                                                                                                 <  1:237247/140‑1 (MQ=255)
        gTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGATAAATAAGATTAAAAGAAGTGAAGGATGGTATAAATTTATCCTTCGCTTATTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTaa                                                                                                                                 >  2:284538/1‑140 (MQ=255)
         tATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGATAAATAAGATTAAAAGAAGTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGTGGCGTaaa                                                                                                                                >  2:453080/1‑140 (MQ=255)
           tttAAAAGAAGTAGTTGAACGAGATAAACAAAATACTGAAGATAACTAAGATCAAAAGAAGTGAAGTATGTTAGAAATTCATCCTTCGCTTCTTTTTATTAATTGAGTAATGAATAGTAGAAAGAAAAGATGCGTaaaaa                                                                                                                              <  2:135000/140‑1 (MQ=255)
           tttAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGATAAATAAGATTAAAAGAAGTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTaaaaa                                                                                                                              <  2:717849/140‑1 (MQ=255)
           tttAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGATAAATAAGATTAAAAGAAGTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTaaaaa                                                                                                                              >  1:85055/1‑140 (MQ=255)
                  gaagTACTTGAACGAGATAAACAAAATACTGAAGATAAATAAGATTAAAAGAAGTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTAt                                                                                                                       <  2:274190/140‑1 (MQ=255)
                    agTACTTGAACGAGATAAACAAAATACTGAAGATAAATAAGATTAAAAGAAGTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGt                                                                                                                     >  1:669049/1‑140 (MQ=255)
                    agTACTTGAACGAGATAAACAAAATACTGAAGATAAATAAGATTAAAAGAAGTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGt                                                                                                                     >  1:264217/1‑140 (MQ=255)
                     gTACTTGAACGAGATAAACAAAATACTGAAGATAAATAAGATTAAAAGAAGTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAACGATGCGTAAAAAGAATTATGtt                                                                                                                    >  2:72377/1‑140 (MQ=255)
                       aCTTGAACGAGATAAACAAAATACTGAAGATAAATAAGATTAAAAGAAGTGAAGGATGTTATAAATTTATCCTTCACTTCTATTTATTAGTTTAGTAAAGAATAGTAGAAAGAAAAGATGCGTAAAAAGAAATAGGATaa                                                                                                                  >  2:483394/1‑140 (MQ=255)
                            aaCGAGATAAACAAAATACTGAAGATAAATAAGATTAAAAGAAGTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAg                                                                                                             <  1:140888/140‑1 (MQ=255)
                                aaaTAAACAAAATACTGAAGATAAATAAGATTAAACGAAGTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTc                                                                                                         <  2:475800/138‑1 (MQ=255)
                                  aTAAACAAAATACTGAAGATAAATAAGATTAAAAGAAGTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATGTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCaa                                                                                                       <  2:253263/140‑1 (MQ=255)
                                    aaaCAAAATACTGAAGATAAATAAGATTAAAAGAAGTGAAGGATGTTATAAATTTATCCTTCTCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATc                                                                                                     >  2:519784/1‑140 (MQ=255)
                                     aaCAAAATACTGAAGATAAATAAGATTAAAAGAAGTGAATGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCt                                                                                                    <  2:434575/140‑1 (MQ=255)
                                     aaCAAAATACTGAAGATAAATAAGATTAAAAGAAGTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCt                                                                                                    >  1:779813/1‑140 (MQ=255)
                                        aaaaTACTGAAGATAAATAAGATTTAAAGAAGTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTGAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTaga                                                                                                 <  2:208150/140‑1 (MQ=255)
                                        aaaaCACTGAAGATAAATAAGATTAAAAGAAGTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTaga                                                                                                 <  1:229292/140‑1 (MQ=255)
                                         aaaTACTGAAGATAAATAAGATTAAAAGAAGTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTagag                                                                                                <  2:378899/140‑1 (MQ=255)
                                          aagaCTGAAGATAAATAAGATTAAAAGAAGTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAAATATGTTAAGATAGGGTCAATCTAGagt                                                                                               <  1:521067/137‑1 (MQ=255)
                                                  aGATAAATAAGATTAAAAGAAGTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTCAGTAATGAATAGTAGAAAGGAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAAc                                                                                       >  2:698619/1‑140 (MQ=255)
                                                  aGATAAATAAGATTAAAAGAAGTGAAGGATGTAATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGGAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAAc                                                                                       >  2:702156/1‑140 (MQ=255)
                                                    ataaagaaGATTAAAAGAAGTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGATTTGCTGCTTAGAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACAt                                                                                     <  1:151340/140‑1 (MQ=255)
                                                                   agaagTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAg                                                                      >  2:11430/1‑140 (MQ=255)
                                                                   agaagTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAg                                                                      >  2:417840/1‑140 (MQ=255)
                                                                   agaagTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAg                                                                      >  2:577799/1‑140 (MQ=255)
                                                                    gaagTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGGAAAGATGCGTAAAAAGACTTACATTACTTAAGGGTGAATCTAGAGTAGTGAAACATAAATCGAACTGGGCGt                                                                     >  2:167538/1‑140 (MQ=255)
                                                                     aagTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTg                                                                    >  1:570389/1‑140 (MQ=255)
                                                                     aagTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATATAGAGTAGTTAAACATAAATAGAAATGGGAggg                                                                    >  2:403104/1‑138 (MQ=255)
                                                                         gAAGGATGTTATAAATTTTTCCTTCGCTTCTTTTTATTAATTTAGTAACGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGAc                                                                <  2:118118/140‑1 (MQ=255)
                                                                         gAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAACAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGAc                                                                >  1:772943/1‑140 (MQ=255)
                                                                           aGGATGTTATAAAGTTATCCTTCGCTTCTGTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGAAGCGTAAAGGTAATAGCCTAAAGGTCATACTATGCTAGAGTAGTTAAACATAAATCGAACGGGGAGTGGGACAg                                                              <  2:360978/140‑1 (MQ=255)
                                                                                 ttATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGCAATGa                                                        >  2:372900/1‑140 (MQ=255)
                                                                                 ttATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGa                                                        >  1:409047/1‑140 (MQ=255)
                                                                                  tataAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGAt                                                       <  1:151633/140‑1 (MQ=255)
                                                                                    taAATTTCTCTTGTGTTTTTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGCAGAGTAAAAAGAAGTATGGTAAGAGAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATaa                                                     <  1:157261/140‑1 (MQ=255)
                                                                                    taAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATaa                                                     <  1:526503/140‑1 (MQ=255)
                                                                                      aaTTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAg                                                   >  1:141571/1‑140 (MQ=255)
                                                                                            tCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATATAGAGCAGTTAAACATAAATCGAACTGGGAGTGGGAAAGAAATGCTAAAGAAtaca                                             >  2:441881/1‑137 (MQ=255)
                                                                                            gccTTCGCTGCTTTTTATTAATATAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGCACGATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCAc                                             <  1:113746/139‑1 (MQ=255)
                                                                                             ccTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACt                                            >  1:120439/1‑140 (MQ=255)
                                                                                                tCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAACAGATGCGTAGAAAGAATTAGGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACTAAt                                         <  2:585357/140‑1 (MQ=255)
                                                                                                   cttcgtTGTATTAAGATAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAGACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACTAATGAt                                      <  1:602682/140‑1 (MQ=255)
                                                                                                        ttttATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGAGTCAATCTGGAGTAGTAATAGTGTATCCGAACTGAGAGTGGGACAGAAATGATAAAGAATCACTAATGATttata                                 >  1:221817/1‑139 (MQ=37)
                                                                                                         tttATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACTAGTGATttatta                                >  1:14961/1‑140 (MQ=255)
                                                                                                           tattaATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACTAATGATTTATTATg                              >  2:79160/1‑140 (MQ=255)
                                                                                                            attaATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACTAATGATTTATTATGt                             <  2:728446/140‑1 (MQ=255)
                                                                                                             ttaATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTAGGACAGAAATGATAAAGAATCACTAATGATTTATTATGTa                            <  2:779813/140‑1 (MQ=255)
                                                                                                              taATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACTAATGATTTATTATGTAg                           <  1:541251/140‑1 (MQ=255)
                                                                                                                aTTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGACTCACTAATGATTTATTATGTAGTg                         >  1:549002/1‑140 (MQ=255)
                                                                                                                aTTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGCAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACTAATGACTTATTATGTAGTg                         >  2:647578/1‑140 (MQ=255)
                                                                                                                aTTTAGTAATGAAGAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACTAATGATTTATTATGTAGTg                         <  1:451131/140‑1 (MQ=255)
                                                                                                                  ttattaATGAATAGTATAAAGAAGAGATGGGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACTAATGATTTATTATGTAGTGGt                       <  2:274666/140‑1 (MQ=255)
                                                                                                                   tAGTAATGAATAGTAGAAAGGAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACTAATGATTTATTATGTAGTGGtt                      <  1:698619/140‑1 (MQ=255)
                                                                                                                   tAGTAATGAATAGTAGAAAGGAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACTAATGATTTATTATGTAGTGGtt                      <  1:702156/140‑1 (MQ=255)
                                                                                                                     gTAATGAATAGTAGAAAGAAAAGATGGGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACTAATGATTTATTATGTAGTGGTTCt                    <  1:111288/140‑1 (MQ=255)
                                                                                                                      tAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATGAAGAATCACTAATGATTTATTATGTAGTGGTTCtt                   >  1:217737/1‑140 (MQ=255)
                                                                                                                         tGAATGGGAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACTAATGATTTATTATGTAGTGGTTCTTTGt                <  2:85055/140‑1 (MQ=255)
                                                                                                                          gAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACTAATGATTTATTATGTAGTGGTTCTTTGTc               >  1:595154/1‑140 (MQ=255)
                                                                                                                             tagtagAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGGTAAACATAAATCGAACTGGGAGCGGGACAGAAATGATAAAGAATCACTAATGATTTATTATGTAGTAGTTCTTTGTCAtt            >  2:395269/1‑140 (MQ=255)
                                                                                                                               gtagAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACTAATGATTTATTATGTAGTGGTTCTTTGTCATTAg          >  1:137593/1‑140 (MQ=255)
                                                                                                                               gtagAAAGAAAAGATGCGGAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACTAATGATTTATTATGTAGTGGTTCTTTGTCATTAg          <  2:555327/140‑1 (MQ=255)
                                                                                                                                 agaaagaaAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACTAATGATTTATTATGTAGTGGTTCGTTGTCATTAGcc        >  2:2523/1‑140 (MQ=255)
                                                                                                                                  gaaagaaaAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACTAATGATTTATTATGTAGTGGTTCTTTGTCATTACCca       >  2:471106/1‑140 (MQ=255)
                                                                                                                                   aaagaaaagaTGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACTAATGATTTATTATGTAGTGGTTCTTTGTCATTAGCcac      <  1:417840/140‑1 (MQ=255)
                                                                                                                                   aaagaaaagaTGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACTAATGATTTATTATGTAGTGGTTCTTTGTCATTAGCcac      <  1:11430/140‑1 (MQ=255)
                                                                                                                                     agaaaagaTGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACTAATGATTTATTATGTAGTGGGTCTTTGTCATTAGCCACAg    >  2:428033/1‑140 (MQ=255)
                                                                                                                                      gaaaagaTGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACTAATGATTTATTATGTAGTGGTTCTTTGTCATTAGCCACAGc   <  1:132848/140‑1 (MQ=255)
                                                                                                                                       aaaagaTGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACTAATGATTTATTATGTAGTGGTTCTTTGTCATTAGCCACAGCt  >  2:470249/1‑140 (MQ=255)
                                                                                                                                       aaaagaTGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACTAATGATTTATTATGTAGTGGTTCTTTGTCATTAGCCACAGCt  <  2:224364/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                           
ACAGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGATAAATAAGATTAAAAGAAGTGAAGGATGTTATAAATTTATCCTTCGCTTCTTTTTATTAATTTAGTAATGAATAGTAGAAAGAAAAGATGCGTAAAAAGAATTATGTTAAGATAGGGTCAATCTAGAGTAGTTAAACATAAATCGAACTGGGAGTGGGACAGAAATGATAAAGAATCACTAATGATTTATTATGTAGTGGTTCTTTGTCATTAGCCACAGCT  >  USA300TCH1516_ALE/842555‑842829

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: