Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,653,335 C→T 9.8% intergenic (‑40/+119) gcvT ← / ← aroK Aminomethyltransferase/Shikimate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,653,3350CT 9.8% 435.9 / 27.5 143intergenic (‑40/+119)gcvT/aroKAminomethyltransferase/Shikimate kinase
Reads supporting (aligned to +/- strand):  ref base C (64/65);  new base T (5/9);  total (69/74)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 146 total aligned reads displayed.

TATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACCTTCG  >  USA300TCH1516_ALE/1653247‑1653425
                                                                                        |                                                                                            
tATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGc                                                                                            <  1:1047517/91‑1 (MQ=255)
 aTCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCt                                                                                           <  1:2442434/91‑1 (MQ=255)
  tCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTa                                                                                          <  1:3111419/91‑1 (MQ=255)
  tCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTa                                                                                          <  1:4090326/91‑1 (MQ=255)
   cAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTAc                                                                                         <  1:873771/91‑1 (MQ=255)
    aaCATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTAcc                                                                                        >  1:2920391/1‑91 (MQ=255)
     aCATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTAccc                                                                                       <  1:2155272/91‑1 (MQ=255)
      cATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCt                                                                                      >  1:3664048/1‑91 (MQ=255)
      cATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCt                                                                                      >  1:3816958/1‑91 (MQ=255)
      cATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCt                                                                                      >  1:468537/1‑91 (MQ=255)
      cATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCt                                                                                      >  1:2788250/1‑91 (MQ=255)
      cATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCt                                                                                      >  1:587182/1‑91 (MQ=255)
      cATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCt                                                                                      >  1:2009801/1‑91 (MQ=255)
       aTAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTATCCTg                                                                                     <  1:1963208/91‑1 (MQ=38)
       aTAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTg                                                                                     <  1:354466/91‑1 (MQ=255)
         aaTTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑tt                                                                                 <  1:2606434/91‑1 (MQ=255)
         aaTTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑tt                                                                                 <  1:1912391/91‑1 (MQ=255)
          aTTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑ttt                                                                                <  1:3464345/91‑1 (MQ=255)
          aTTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑ttt                                                                                >  1:2033368/1‑91 (MQ=255)
          aTTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑ttt                                                                                >  1:2033117/1‑91 (MQ=255)
          aTTTTGATATAAAGGTGTTTGTTTTAAATCACTAGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑ttt                                                                                <  1:3402986/91‑1 (MQ=255)
           ttttGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTATCCTG‑‑tttt                                                                               <  1:20136/91‑2 (MQ=37)
            tttGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTtat                                                                              <  1:2767880/91‑1 (MQ=255)
             ttGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTtata                                                                             >  1:2839161/1‑91 (MQ=255)
              tGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTATCCTG‑‑tttttat                                                                            <  1:1713917/91‑5 (MQ=25)
              tGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATaa                                                                            >  1:1909843/1‑91 (MQ=255)
               gATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATaaa                                                                           >  1:1751045/1‑91 (MQ=255)
                atatAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAAt                                                                          <  1:654819/91‑1 (MQ=255)
                 tataAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATc                                                                         <  1:1356832/91‑1 (MQ=255)
                 tataAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATc                                                                         <  1:1361287/91‑1 (MQ=255)
                 tataAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATc                                                                         <  1:1626350/91‑1 (MQ=255)
                  ataAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATTAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCa                                                                        <  1:3507595/91‑1 (MQ=255)
                  ataAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCa                                                                        <  1:3932408/91‑1 (MQ=255)
                   taAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCAc                                                                       <  1:1796020/91‑1 (MQ=255)
                    aaaGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCAcc                                                                      <  1:2454614/91‑1 (MQ=255)
                     aaGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCa                                                                     <  1:1126763/91‑1 (MQ=255)
                     aaGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCa                                                                     >  1:1829450/1‑91 (MQ=255)
                      aGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCaa                                                                    >  1:1795288/1‑91 (MQ=255)
                       ggTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAAt                                                                   >  1:1643849/1‑91 (MQ=255)
                        gtgtTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAAtt                                                                  <  1:2697116/91‑1 (MQ=255)
                        gtgtTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAAtt                                                                  >  1:2940992/1‑91 (MQ=255)
                            ttgtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAgg                                                              <  1:3392384/91‑1 (MQ=255)
                             tgtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGa                                                             >  1:571087/1‑91 (MQ=255)
                             tgtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGa                                                             >  1:3175456/1‑91 (MQ=255)
                             tgtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGa                                                             >  1:2014224/1‑91 (MQ=255)
                              gtttTAAATCACCTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCCATTCAGGaa                                                            <  1:3626759/91‑1 (MQ=255)
                                tttAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAAcg                                                          <  1:2080474/91‑1 (MQ=255)
                                 ttAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAAcgc                                                         >  1:589107/1‑91 (MQ=255)
                                 ttAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAAcgc                                                         >  1:493278/1‑91 (MQ=255)
                                 ttAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAAcgc                                                         <  1:3760911/91‑1 (MQ=255)
                                 ttAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAAcgc                                                         <  1:3584652/91‑1 (MQ=255)
                                 ttAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAAcgc                                                         <  1:1152935/91‑1 (MQ=255)
                                 ttAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAAcgc                                                         <  1:92558/91‑1 (MQ=255)
                                        aCTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAg                                                  <  1:2614523/91‑1 (MQ=255)
                                        aCTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAg                                                  >  1:2271992/1‑91 (MQ=255)
                                         cTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGt                                                 >  1:1496139/1‑91 (MQ=255)
                                          ttGACATTGAAATACCCCCTTAAAATTATAAAACCAGAATAGCAATCGCTACCCTG‑‑TTTGTAAATCACCAATTCAGGAACGCGTTACAGta                                                >  1:3312384/1‑91 (MQ=255)
                                          ttGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGta                                                >  1:611283/1‑91 (MQ=255)
                                           tGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGtat                                               <  1:832167/91‑1 (MQ=255)
                                           tGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGtat                                               <  1:2087914/91‑1 (MQ=255)
                                           tGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGtat                                               <  1:817941/91‑1 (MQ=255)
                                                 tGAAATACCCCCTTAAAATTATAAAATCAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCt                                         <  1:3480349/91‑1 (MQ=255)
                                                  gAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCtt                                        >  1:4106708/1‑91 (MQ=255)
                                                  gAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCtt                                        >  1:2912073/1‑91 (MQ=255)
                                                  gAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCtt                                        >  1:2911654/1‑91 (MQ=255)
                                                  gAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCtt                                        >  1:2120156/1‑91 (MQ=255)
                                                     aTACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTg                                     >  1:4037653/1‑91 (MQ=255)
                                                      tACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGc                                    >  1:775744/1‑91 (MQ=255)
                                                        cccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCt                                  >  1:828427/1‑91 (MQ=255)
                                                        cccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCt                                  <  1:995491/91‑1 (MQ=255)
                                                        cccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCt                                  <  1:3051092/91‑1 (MQ=255)
                                                         ccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTg                                 >  1:3569011/1‑91 (MQ=255)
                                                         ccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTg                                 >  1:4010462/1‑91 (MQ=255)
                                                         ccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTg                                 <  1:1538075/91‑1 (MQ=255)
                                                             ttAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGagag                             >  1:1852574/1‑91 (MQ=255)
                                                             ttAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGagag                             <  1:489156/91‑1 (MQ=255)
                                                              tAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGt                            <  1:985691/91‑1 (MQ=255)
                                                                  aTTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCg                        >  1:3042725/1‑91 (MQ=255)
                                                                  aTTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCg                        >  1:623864/1‑91 (MQ=255)
                                                                   ttATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGt                       <  1:1392881/91‑1 (MQ=255)
                                                                      tataaaCAGGGTAGCGATTGCTATCCTGTTTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGta                      >  1:2477099/4‑91 (MQ=18)
                                                                       aaaaaCAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGtaata                   <  1:2744654/91‑1 (MQ=255)
                                                                            cAGGGTAGCGATTGCTATCCTG‑‑TTTTTAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTa              >  1:612036/1‑91 (MQ=37)
                                                                            cAGGGTAGCGATTGCTATCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTa              >  1:2339578/1‑91 (MQ=38)
                                                                            cAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTa              <  1:1074867/91‑1 (MQ=255)
                                                                              ggATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACt            <  1:3643245/91‑1 (MQ=255)
                                                                               gATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGTAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACtt           <  1:187604/91‑1 (MQ=255)
                                                                               gATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACtt           >  1:3356106/1‑91 (MQ=255)
                                                                               gATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACtt           >  1:1178665/1‑91 (MQ=255)
                                                                                gtAGCGATTGCTATCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTCTCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTg          >  1:2424033/2‑91 (MQ=255)
                                                                                gtAGCGATTGCTATCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTg          >  1:4030216/2‑91 (MQ=255)
                                                                                aTAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTg          >  1:469029/1‑91 (MQ=255)
                                                                                aTAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTg          >  1:1778296/1‑91 (MQ=255)
                                                                                  aGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCa        >  1:4120030/1‑91 (MQ=255)
                                                                                      aTCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACCtt    >  1:3972366/1‑91 (MQ=255)
                                                                                       tCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACCTTc   >  1:4123005/1‑91 (MQ=255)
                                                                                       tCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACCTTc   <  1:284805/91‑1 (MQ=255)
                                                                                        cGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACCTTCg  >  1:2935386/1‑91 (MQ=255)
                                                                                        cGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACCTTCg  <  1:1370375/91‑1 (MQ=255)
                                                                                        cGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACCTTCg  >  1:2383860/1‑91 (MQ=255)
                                                                                        |                                                                                            
TATCAACATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTG‑‑TTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACCTTCG  >  USA300TCH1516_ALE/1653247‑1653425

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