Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,386,969 A→T 6.9% K117* (AAA→TAA)  USA300HOU_RS12255 → putative HTH‑type transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,386,9690AT 6.9% 221.5 / 5.7 68K117* (AAA→TAA) USA300HOU_RS12255putative HTH‑type transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (24/39);  new base T (3/2);  total (27/41)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAAATGACCTTATAACAGTAGAAAACAGTTTAGAGGATAAACGACAAAAGTTTTTAACTTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAATATTATTAAGGGAGGCACACAATGACTAATCAAGGCTCATTAAATAAATCACCTATTTTCACCAAAAGTTTTA  >  USA300TCH1516_ALE/2386834‑2387105
                                                                                                                                       |                                                                                                                                        
gAAAATGACCTTATAACAGTAGAAAACAGTTTAGAGGATAAACGACAAAAGTTTTTAACTTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAaaaca                                                                                                                                     <  2:714024/141‑1 (MQ=255)
  aaaTGACCTTATAACAGTAGAAAACAGTTTAGAGGATAAACGACAAAAGTTTTTAACTTTAACACCTAAAGGGCATAAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAaaacaaa                                                                                                                                   >  2:494611/1‑141 (MQ=255)
   aaTGACCTTATAACAGTAGAAAACAGTTTAGAGGATAAACGACAAAAGTTTTTAACTTTAACACCTAAAGGGCATGAATTATATGAGATTATTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAaaacaaac                                                                                                                                  <  1:508528/141‑1 (MQ=255)
         cTTATAAAAGGAGAAAACAGTTGAGAGGATAAACCACAAGAGTTTTTAACTTTAACACCCAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCg                                                                                                                            <  2:458084/141‑1 (MQ=255)
          ttATACAAGTAGAAAACAGTTTAGAGGATAAACGACAAAAGATTTTCCATTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCACCAAGCAGTAGTTGCCAAAACAAACATTTCGc                                                                                                                           <  1:59316/141‑1 (MQ=255)
          ttATAACAGTAGAAAACAGTTTAGAGGATAAACGACAAAAGTTTTTAACTTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGc                                                                                                                           <  2:713530/141‑1 (MQ=255)
          ttATAACAGTAGAAAACAGTTTAGAGGATAAACGACAAAAGTTTTTAACTTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGc                                                                                                                           <  1:436833/141‑1 (MQ=255)
            ataACAGTAGAAAACAGTTTAGAGGATAAACGACAAAAGTTTTTAACTTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCaa                                                                                                                         <  2:416845/141‑1 (MQ=255)
             taACAGTAGAAAACAGTTTAGAGGATAAACGACAAAAGTTTTTAACTTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAg                                                                                                                        <  1:35575/141‑1 (MQ=255)
                  gTAGAAAACAGTTTAGAGGATAAACGACAAAAGTTTTTAACTTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCaa                                                                                                                   >  1:331241/1‑141 (MQ=255)
                      aaaaCAGTTTAGAGGATAAACGACAAAAGTTTTTAACTTTAACAGCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGc                                                                                                               >  1:435284/1‑141 (MQ=255)
                      aaaaCAGTTTAGAGGATAAACGACAAAAGTTTTTAACTTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGc                                                                                                               <  2:103516/141‑1 (MQ=255)
                        aaCAGTTTAGAGGATAAACGACAAAAGTTTTTAACTTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCaa                                                                                                             <  2:389340/141‑1 (MQ=255)
                         actgTTTAGAGGATAAACTACTATAGTTTTTAACTTTAACACCTAAAGGGCATCAATTATATTAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAg                                                                                                            <  2:421076/138‑1 (MQ=255)
                          cAGTTTAGAGGATAAACGACAAAAGTTTTTAACTTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGa                                                                                                           <  1:456447/141‑1 (MQ=255)
                               tAGAGGATAAACGACAAAAGTTTTTAACTTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCa                                                                                                      >  2:581966/1‑141 (MQ=255)
                               tAGAGGATAAACGACAAAAGCAATGCACGTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCa                                                                                                      >  1:564162/1‑141 (MQ=255)
                                agagGATAAACGACAAAAGTTTTTAACTTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCAt                                                                                                     >  1:99169/1‑141 (MQ=255)
                                agagGATAAACGACAAAAGTTTTTAACTTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATAAAATGCAAGAAACCAt                                                                                                     >  1:694588/1‑141 (MQ=255)
                                   ggATAAACGGCAAAAGTTTTTAACTTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATaaa                                                                                                  <  2:697261/141‑1 (MQ=255)
                                   ggATAAACGACAAAAGTTTTTAACTTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCCATAATTACTATTTGGCACCGCAAACAGTTCGCAAGATCAAATGCAAGAAACCATaaa                                                                                                  <  1:369425/141‑1 (MQ=255)
                                    gATAAACGACAAAAGTTTTTAACTTTAACACCTAACGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAACGCTCCAACAATCAGTAGTGACCAAAACAATCATTTCGCAAGATCAAATGCAAGAAACCATAAAt                                                                                                 <  2:409189/141‑1 (MQ=255)
                                       aaaCGACAAAAGTTTTTAACTTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAAtgtg                                                                                              <  1:278841/141‑1 (MQ=255)
                                         aCGACAAAAGTTTTTAACTTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGACACCATAAATGtgat                                                                                            >  2:322840/1‑141 (MQ=255)
                                           gACAAAAGTTTTTAACTTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGtgatga                                                                                          >  2:339326/1‑141 (MQ=255)
                                                 aGTTTTTAACTTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCaaa                                                                                    >  1:321337/1‑141 (MQ=255)
                                                  gTTTTTAACTTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCTAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAAt                                                                                   >  1:536812/1‑141 (MQ=255)
                                                  gTTTTTAACTTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCTAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAAt                                                                                   >  1:646523/1‑141 (MQ=255)
                                                  gTTTTTAACTTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAAt                                                                                   >  2:649733/1‑141 (MQ=255)
                                                  gTTTTTAACTTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGGGATGAATCAAAt                                                                                   >  1:316696/1‑141 (MQ=255)
                                                  gTTTTTAACTTTAACACCTAAAGGGCATGAATTATATGAACTTGCTTCTCTTGCTGTACAAAACCTCCAACAAGCAGTAGTTCCCTAGACAAACATTTCGCAAGATCGAATGCAAGCAACCATCAAACTGATGAACCAAAt                                                                                   >  2:422300/1‑141 (MQ=255)
                                                    ttttAACTTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTc                                                                                 >  2:472672/1‑141 (MQ=255)
                                                       taactttaacACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATg                                                                              <  2:704105/141‑1 (MQ=255)
                                                       taactttaacACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATg                                                                              <  2:628982/141‑1 (MQ=255)
                                                       taactttaacACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATg                                                                              <  2:704086/141‑1 (MQ=255)
                                                           tttaacACCTAAAGGGCATGAATTATATTAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAat                                                                          <  1:709577/141‑1 (MQ=255)
                                                            ttaacACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAata                                                                         >  2:74064/1‑141 (MQ=255)
                                                               acacCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAAtatta                                                                      <  1:60573/141‑1 (MQ=255)
                                                                cacCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAAtattat                                                                     <  1:20361/141‑1 (MQ=255)
                                                                 acCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAATattatt                                                                    <  2:435284/141‑1 (MQ=255)
                                                                        ggtcATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAATATTATTAAGGGAg                                                             <  1:709299/138‑1 (MQ=255)
                                                                           cATGAATTATATGAGATTGTTTGTGTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAATCCATAAATGTGATGAATCAAATTCATGAAATATTATTAAGGGAGGca                                                          <  1:52374/141‑1 (MQ=255)
                                                                                aTTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAATATTATTAAGGAAGGCACACaa                                                     >  2:717645/1‑141 (MQ=255)
                                                                                     atGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAATATTATTAAGGGAGGCACACAATGACt                                                <  1:649733/141‑1 (MQ=255)
                                                                                       aagaTTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAATATTATTAAGGGAGGCACACAATGACTaa                                              <  1:593420/140‑1 (MQ=255)
                                                                                            tgtttggCTTGATGGACCAAACCTCCAACAAGCAGTCGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAATATTATTAAGGGAGGCACACAATGACTAATCAAg                                         <  1:237593/141‑1 (MQ=255)
                                                                                              tttgtCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAATATTATTAAGGGAGGCACACAATGACTAATCAAGGc                                       <  2:331241/141‑1 (MQ=255)
                                                                                               ttgtCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAATATTATTAAGGGAGGCACACAATGACTAATCAAGGCt                                      <  1:165930/141‑1 (MQ=255)
                                                                                               ttgtCTTGATGTACAAAAGATACACCAATCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAATATTATTAAGGGAGGCACACAATGACTAATCAAAGca                                      >  2:738361/1‑140 (MQ=255)
                                                                                                  tCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAATATTATTAAGGGAGGCACACAATGACTAATCAAGGCTCAt                                   >  1:263644/1‑141 (MQ=255)
                                                                                                      gATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAATATTATTAAGGGAGGCACACAATGACTAATCAAGGCTCATtaaa                               <  1:190155/141‑1 (MQ=255)
                                                                                                       aTGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAATATTATTAAGGGAGGCACACAATGACTAATCAAGGCTCATtaaaa                              >  2:707102/1‑140 (MQ=255)
                                                                                                       aTGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGAAAGAAACCATAAATGTGATGAATCAAATTCATGAAATATTATTAAGGGAGGCACACAATGACTAATCAAGGCTCATtaaat                              <  2:746573/141‑1 (MQ=255)
                                                                                                         gTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAATATTATTAAGGGAGGCACACAATGACTAATCAAGGCTCATtaaataa                            <  1:664875/141‑1 (MQ=255)
                                                                                                         gTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAATATTATTAAGGGAGGCACACAATGACTAATCAAGGCTCATtaaataa                            <  1:435190/141‑1 (MQ=255)
                                                                                                              aaaGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAATATTATTAAGGGAGGCACACAATGACTAATCAAGGCTCATTAAATAAATCAc                       <  1:735375/141‑1 (MQ=255)
                                                                                                                      aacaaGCAGTAGTTGCCAAAACAAACATTTCGCAAGCTCAAATGCAAGAAACCATAAATGAGACGAATCAAATTCATGAAATATTATTAAGGGAGGCACACAATGACTAATCAAGGCTCATTAAATAAATCACCTATTTTc               <  1:519201/141‑1 (MQ=255)
                                                                                                                       acaaGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAATATTATTAAGGGAGGCACACAATGACTAATCAAGGCTCATTAAATAAATCACCTATTTTCa              >  2:137496/1‑141 (MQ=255)
                                                                                                                        caaGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAATATTATTAAGGGAGGCACACAATGACTAATCAAGGCTCATTAAATAAATCACCTATTTTCAc             >  2:558821/1‑141 (MQ=255)
                                                                                                                          aGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGTAATATTATTAAGGGAGGCACACAATGACTAATCCTGGCTTGCGAAATAAATCACCTATTTGCACCa           <  1:563957/141‑1 (MQ=255)
                                                                                                                          aGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAATATTATTAAGGGAGGCACACAATGACTAATCAAGGCTCATTAAATAAATCACCTATTTTCACCa           <  1:732366/141‑1 (MQ=255)
                                                                                                                           gCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATTAATGTGATGAATCAAATTCATGAAATATTATTAAGGGAGGCACACAATGACTAATCAAGGCTCATTAAATAAATCACCTATTTTCACCaa          >  2:564265/1‑141 (MQ=255)
                                                                                                                           gCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAAGATTATTAAGGGAGGCACACAATGACTAATCAAGGCTCATTAAATTAATTACCTATTTTCACCaa          <  1:572445/141‑1 (MQ=255)
                                                                                                                            cAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAATATTATTAAGGGAGGCACACAATGACTAATCAAGGCTCATTAAATAAATCACCTATTTTCACCaaa         >  1:650567/1‑141 (MQ=255)
                                                                                                                            cAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAATATTATTAAGGGAGGCACACAATGACTAATCAAGGCTCATTAAATAAATCACCTATTTTCACCaaa         >  1:148746/1‑141 (MQ=255)
                                                                                                                               tagtTGCCTAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAATATTATTAAGGGAGGCACACAATGACTAATCAAGGCTCATTAAATAAATCACCTATTTTCACCAAAAGt      <  2:536812/141‑1 (MQ=255)
                                                                                                                               tacgtGCCTAAACAAACATTTCGCAAGATCAAATGAAAGAAACCATAAATGTGATGAATCAAATTCATGAAATATTATTAAGGGAGGCACACAATGACTAATCAAGGCTCATTAAATAAATCACCTATTTTCACCAAAAGt      <  2:646523/137‑1 (MQ=255)
                                                                                                                                   tGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAATATTATTAAGGGAGGCACACAATGACTAATCAAGGCTCATTAAATAAATCACCTATTTTCACCAAAAGTTTTa  >  1:102914/1‑141 (MQ=255)
                                                                                                                                       |                                                                                                                                        
GAAAATGACCTTATAACAGTAGAAAACAGTTTAGAGGATAAACGACAAAAGTTTTTAACTTTAACACCTAAAGGGCATGAATTATATGAGATTGTTTGTCTTGATGTACAAAAGCTCCAACAAGCAGTAGTTGCCAAAACAAACATTTCGCAAGATCAAATGCAAGAAACCATAAATGTGATGAATCAAATTCATGAAATATTATTAAGGGAGGCACACAATGACTAATCAAGGCTCATTAAATAAATCACCTATTTTCACCAAAAGTTTTA  >  USA300TCH1516_ALE/2386834‑2387105

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: