Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 950,098 T→A 16.2% intergenic (+92/‑270) yugI_2 → / → USA300HOU_RS04740 General stress protein 13/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE950,0980TA16.2% 629.3 / 101.9 228intergenic (+92/‑270)yugI_2/USA300HOU_RS04740General stress protein 13/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (83/108);  new base A (24/13);  total (107/121)
Fisher's exact test for biased strand distribution p-value = 1.96e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 241 total aligned reads displayed.

ACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATCT  >  USA300TCH1516_ALE/950011‑950188
                                                                                       |                                                                                          
aCAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAg                                                                                         <  2:2820965/91‑1 (MQ=38)
    aTAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATAc                                                                                     >  1:1695555/1‑91 (MQ=255)
    aTAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATAc                                                                                     <  1:2373653/91‑1 (MQ=37)
     tAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACt                                                                                    <  2:3752731/91‑1 (MQ=255)
     tAAAGCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATACt                                                                                    <  1:3230695/91‑1 (MQ=25)
      aaaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATACTa                                                                                   <  2:1666655/91‑1 (MQ=37)
      aaaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATACTa                                                                                   <  2:1381469/91‑1 (MQ=37)
       aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                                  <  2:1958711/91‑1 (MQ=255)
       aaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                                  <  2:181695/91‑1 (MQ=255)
        aTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAAt                                                                                 <  2:140971/91‑1 (MQ=255)
         tCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAAtc                                                                                <  1:3446993/91‑1 (MQ=255)
          cGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGTCTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTCAtct                                                                               <  1:3069923/91‑1 (MQ=255)
             aCCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCtt                                                                            >  2:3276213/1‑91 (MQ=255)
               cGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCtttt                                                                          >  1:2746067/1‑91 (MQ=255)
                gAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCttttt                                                                         <  2:905941/91‑1 (MQ=255)
                gAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCttttt                                                                         <  1:2512430/91‑1 (MQ=255)
                   aaTCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTc                                                                      <  1:1121494/91‑1 (MQ=255)
                    aTCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCt                                                                     >  1:3873859/1‑91 (MQ=255)
                      cATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTAc                                                                   <  1:2342577/91‑1 (MQ=255)
                       aTACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACg                                                                  >  1:772711/1‑91 (MQ=255)
                       aTACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACg                                                                  <  2:2484589/91‑1 (MQ=255)
                        tACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGa                                                                 >  1:3756644/1‑91 (MQ=255)
                          cAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATa                                                               >  1:903049/1‑91 (MQ=255)
                          cAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATa                                                               >  2:3473744/1‑91 (MQ=255)
                           aaaGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATaa                                                              >  1:1909151/1‑91 (MQ=255)
                              gggTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGt                                                           <  1:1442180/91‑1 (MQ=255)
                              gggTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTACTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGt                                                           >  2:3194610/1‑91 (MQ=255)
                                 tCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAat                                                        >  2:1403897/1‑91 (MQ=255)
                                 tCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAat                                                        >  2:1442164/1‑91 (MQ=255)
                                 tCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAat                                                        >  1:3389538/1‑91 (MQ=255)
                                  cTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAata                                                       <  1:826289/91‑1 (MQ=255)
                                  cTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAata                                                       >  2:2074456/1‑91 (MQ=255)
                                    gaaatgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtg                                                     >  2:3391594/1‑91 (MQ=255)
                                      aatgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGCTACTAATCTCTCTTAGTCTACTATAACGTTATAagat                                                   >  2:2529874/1‑91 (MQ=37)
                                      aatgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgat                                                   >  1:2373758/1‑91 (MQ=255)
                                       atgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgatt                                                  >  1:689678/1‑91 (MQ=255)
                                         gaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattga                                                >  1:2260573/1‑91 (MQ=255)
                                         gaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattga                                                <  1:377702/91‑1 (MQ=255)
                                          aaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattgat                                               >  2:3861849/1‑91 (MQ=255)
                                          aaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattgat                                               <  2:3510886/91‑1 (MQ=255)
                                          aaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattgat                                               >  1:582449/1‑91 (MQ=255)
                                          aaaGTTTCTTAGACTATAAAAGAGATAAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgattgat                                               >  2:2211787/1‑91 (MQ=255)
                                           aaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATgattgatt                                              <  1:2649181/91‑1 (MQ=255)
                                            aGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTc                                             >  2:2367459/1‑91 (MQ=255)
                                            aGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTc                                             >  2:1144353/1‑91 (MQ=255)
                                             gTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCt                                            <  1:1107879/91‑1 (MQ=255)
                                              tttCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTa                                           >  1:462743/1‑91 (MQ=255)
                                              tttCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTa                                           >  1:345534/1‑91 (MQ=255)
                                              tttCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTa                                           >  2:1888792/1‑91 (MQ=255)
                                              tttCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTa                                           <  1:2210/91‑1 (MQ=255)
                                                 cTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTAttt                                        >  2:378627/1‑91 (MQ=255)
                                                 cTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTAttt                                        >  2:569462/1‑91 (MQ=255)
                                                   tAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTac                                      >  1:2810474/1‑91 (MQ=255)
                                                    aGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTaca                                     <  2:587527/91‑1 (MQ=255)
                                                     gACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTacac                                    <  2:1276198/91‑1 (MQ=255)
                                                     gACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTacac                                    <  1:2023861/91‑1 (MQ=255)
                                                      aCTATAAAAGATATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACg                                   <  1:3523625/91‑1 (MQ=255)
                                                       cTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGt                                  <  2:2388831/91‑1 (MQ=255)
                                                          taAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACa                               <  2:3649684/91‑1 (MQ=255)
                                                          aaaaaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACa                               <  1:2559084/90‑1 (MQ=255)
                                                           aaaaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACaa                              >  2:2021279/1‑91 (MQ=255)
                                                             aaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTCCAAAt                            >  2:3581037/1‑91 (MQ=255)
                                                              agagATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATg                           <  2:1706747/91‑1 (MQ=255)
                                                                agaTTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGt                         <  2:1518035/91‑1 (MQ=255)
                                                                agaTTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGt                         <  2:1767973/91‑1 (MQ=255)
                                                                agaTTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGt                         <  2:2444357/91‑1 (MQ=255)
                                                                 gaTTATTATCTATTAAATTTTATTATATACTAATCTCTTTTTGTCTACGATAATGTAATATGATTGATTCTATTTACACGTACAAATGGtt                        <  1:3608705/91‑1 (MQ=255)
                                                                 gaTTAGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACATAATATGATTGATTCTATTTACACGTACAAATGGtt                        >  2:3063457/1‑91 (MQ=25)
                                                                  aTTAGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGttt                       >  2:259194/1‑91 (MQ=25)
                                                                  aTTAGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGttt                       >  2:1301229/1‑91 (MQ=25)
                                                                  aTTAGTATCTAATAAAATTTAATAGATACTAATCTCTCTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGttt                       >  1:3276058/1‑91 (MQ=25)
                                                                   ttAGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTa                      >  1:1063653/1‑91 (MQ=25)
                                                                    tAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTaa                     >  1:954078/1‑91 (MQ=255)
                                                                     aGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAg                    <  2:1046604/91‑1 (MQ=255)
                                                                     aGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAg                    >  2:2783253/1‑91 (MQ=37)
                                                                     aGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAg                    >  1:57298/1‑91 (MQ=37)
                                                                      gTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAgg                   >  2:2851357/1‑91 (MQ=37)
                                                                        aTCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTg                 >  1:499666/1‑91 (MQ=37)
                                                                         tCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGa                <  1:1412575/91‑1 (MQ=255)
                                                                         tCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGa                >  1:2025186/1‑91 (MQ=255)
                                                                           tATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACa              >  2:271810/1‑91 (MQ=255)
                                                                            aTTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGatat             <  2:943960/91‑4 (MQ=255)
                                                                             atAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACata            >  1:1969883/2‑91 (MQ=255)
                                                                              tAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACatat           >  1:957918/1‑91 (MQ=255)
                                                                                aaTTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTGATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATcc         <  1:3165656/91‑1 (MQ=255)
                                                                                aaTTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATcc         <  1:2874073/91‑1 (MQ=255)
                                                                                 aTTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCa        >  1:2538534/1‑91 (MQ=255)
                                                                                 aTTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCa        <  1:1525600/91‑1 (MQ=255)
                                                                                   tttATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATTTGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACTTATCCAtt      <  2:81851/91‑1 (MQ=255)
                                                                                   tttATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCAtt      <  2:146512/91‑1 (MQ=255)
                                                                                   tttATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCAtt      <  1:1465851/91‑1 (MQ=255)
                                                                                   tttATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCAtt      <  1:1730076/91‑1 (MQ=255)
                                                                                   tttATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCAtt      <  2:3821035/91‑1 (MQ=255)
                                                                                   tttATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCAta      <  2:70081/91‑2 (MQ=255)
                                                                                    ttattaGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACGAATGGTTTAAGGTGACATATCCAtaa     <  1:1885054/91‑3 (MQ=255)
                                                                                    ttattaGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTa     <  2:419409/91‑1 (MQ=255)
                                                                                    ttattaGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTa     <  2:865016/91‑1 (MQ=255)
                                                                                      attaGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATc   >  1:1774291/1‑91 (MQ=255)
                                                                                      attaGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATc   <  1:2620939/91‑1 (MQ=255)
                                                                                       ttaGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATCt  >  2:301039/1‑91 (MQ=255)
                                                                                       |                                                                                          
ACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATCT  >  USA300TCH1516_ALE/950011‑950188

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