Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,653,340 C→T 8.0% intergenic (‑45/+114) gcvT ← / ← aroK Aminomethyltransferase/Shikimate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,653,3400CT 8.0% 689.7 / 30.6 221intergenic (‑45/+114)gcvT/aroKAminomethyltransferase/Shikimate kinase
Reads supporting (aligned to +/- strand):  ref base C (87/115);  new base T (13/5);  total (101/120)
Fisher's exact test for biased strand distribution p-value = 2.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 231 total aligned reads displayed.

CATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACCTTCGGCGAT  >  USA300TCH1516_ALE/1653253‑1653430
                                                                                       |                                                                                          
cATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCt                                                                                         <  1:2855836/91‑1 (MQ=255)
   ttttttTATATAAAGGTGTTTGTTTTAAATCACTTGACATCGATATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGtt                                                                                      >  2:442847/3‑91 (MQ=255)
    aTTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTATCCTGttt                                                                                     <  1:2392659/91‑1 (MQ=38)
    aTTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGttt                                                                                     <  1:2343081/91‑1 (MQ=255)
    aTTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGttt                                                                                     <  1:3665012/91‑1 (MQ=255)
     ttttGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTATCCTGtttt                                                                                    <  2:416652/91‑2 (MQ=37)
     ttttGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTta                                                                                    >  1:2873006/1‑91 (MQ=255)
       ttGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTtata                                                                                  >  1:3058110/1‑91 (MQ=255)
        tGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATaa                                                                                 <  2:2906846/91‑1 (MQ=255)
        tGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATaa                                                                                 >  1:2990009/1‑91 (MQ=255)
        tGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATaa                                                                                 >  1:1250370/1‑91 (MQ=255)
        tGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATaa                                                                                 >  2:1919316/1‑91 (MQ=255)
        tGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATaa                                                                                 >  1:1784886/1‑91 (MQ=255)
        tGATATAAAGGTGTTTGTTATAAATCACTTGGCATTGAAATACCCCCTTAAAATTATAAAAACAGGGTAGCGATTGCTATCCTGTTtttat                                                                                 <  2:2600321/91‑2 (MQ=21)
          atatAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAAt                                                                               <  1:2814374/91‑1 (MQ=255)
          atatAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAAt                                                                               <  2:158894/91‑1 (MQ=255)
           tataAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATc                                                                              <  2:3039809/91‑1 (MQ=255)
             taAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCAc                                                                            >  2:184749/1‑91 (MQ=255)
             taAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCAc                                                                            <  1:2568850/91‑1 (MQ=255)
                 ggTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAAt                                                                        <  1:3537128/91‑1 (MQ=255)
                    gtttgtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTcg                                                                     <  1:2138002/91‑2 (MQ=255)
                       tgtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTAGTCTGTTTATAAATCACCAATTCAGGa                                                                  >  2:3354732/1‑91 (MQ=255)
                       tgtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGa                                                                  <  1:1842948/91‑1 (MQ=255)
                        gtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGaa                                                                 <  2:320320/91‑1 (MQ=255)
                        gtttTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGaa                                                                 >  1:2781441/1‑91 (MQ=255)
                            tAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAAcgcg                                                             >  1:3299578/1‑91 (MQ=255)
                             aaaTCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGt                                                            >  2:3018843/1‑91 (MQ=255)
                             aaaTCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGt                                                            <  1:16156/91‑1 (MQ=255)
                              aaTCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGtt                                                           <  1:209502/91‑1 (MQ=255)
                                tCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTAc                                                         <  2:1857069/91‑1 (MQ=255)
                                tCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTAc                                                         >  2:692250/1‑91 (MQ=255)
                                tCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTAc                                                         <  1:1470432/91‑1 (MQ=255)
                                tCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTAc                                                         <  1:2010971/91‑1 (MQ=255)
                                  aCTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAg                                                       >  2:2113629/1‑91 (MQ=255)
                                   cTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGt                                                      <  1:3401468/91‑1 (MQ=255)
                                     tGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtat                                                    <  2:143650/91‑1 (MQ=255)
                                     tGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtat                                                    <  1:2423866/91‑1 (MQ=255)
                                      gACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtatt                                                   >  1:434129/1‑91 (MQ=255)
                                       aCATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtatta                                                  <  2:1726527/91‑1 (MQ=255)
                                       aCATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtatta                                                  <  1:1229351/91‑1 (MQ=255)
                                        cATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtattat                                                 >  2:657800/1‑91 (MQ=255)
                                        cATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGtattat                                                 <  2:2936064/91‑1 (MQ=255)
                                         aTTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATc                                                <  1:1308702/91‑1 (MQ=255)
                                             aaaTACCGCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCttt                                            <  2:3451226/91‑1 (MQ=255)
                                              aaTACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCtttt                                           >  2:3320798/1‑91 (MQ=255)
                                              aaTACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCtttt                                           >  1:3140892/1‑91 (MQ=255)
                                                 aCCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGcc                                        >  1:379489/1‑91 (MQ=255)
                                                  cccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCt                                       <  2:3669365/91‑1 (MQ=255)
                                                  cccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCt                                       <  2:1076796/91‑1 (MQ=255)
                                                   ccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTg                                      <  1:2184463/91‑1 (MQ=255)
                                                   ccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTg                                      <  1:1298946/91‑1 (MQ=255)
                                                   ccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTg                                      <  1:1298476/91‑1 (MQ=255)
                                                   ccccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTg                                      <  1:808405/91‑1 (MQ=255)
                                                    cccTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTga                                     <  2:2317239/91‑1 (MQ=255)
                                                      cTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTgaga                                   <  2:2367052/91‑1 (MQ=255)
                                                        tAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTACCCTTTTGCCTGAGAGt                                 >  2:1445004/1‑91 (MQ=255)
                                                          aaaTTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGttt                               <  1:1733762/91‑1 (MQ=255)
                                                          aaaTTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGttt                               <  2:386545/91‑1 (MQ=255)
                                                          aaaTTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGttt                               >  1:2780774/1‑91 (MQ=255)
                                                           aaTTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTc                              <  1:843805/91‑1 (MQ=255)
                                                           aaTTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTc                              >  1:3571860/1‑91 (MQ=255)
                                                            aTTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCg                             >  1:2555886/1‑91 (MQ=255)
                                                             ttATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGt                            >  1:1453535/1‑91 (MQ=255)
                                                             ttATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGt                            <  1:2521778/91‑1 (MQ=255)
                                                              tataAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGta                           >  1:1133919/1‑91 (MQ=255)
                                                              tataAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGta                           <  1:3170204/91‑1 (MQ=255)
                                                               ataAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGtaa                          >  1:1698498/1‑91 (MQ=255)
                                                                taAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGtaat                         >  1:1009392/1‑91 (MQ=255)
                                                                  taaaCAGGGTAGCGATTGCTATCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGtaataa                       >  2:2516290/2‑91 (MQ=37)
                                                                  aaaaCAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGtaataa                       >  2:2157385/1‑91 (MQ=255)
                                                                   aaaCAGGGTAGCGATTGCTATCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATaaa                      >  2:734171/1‑91 (MQ=38)
                                                                    aaCAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAAt                     <  1:1359989/91‑1 (MQ=255)
                                                                     aCAGGGTAGCGATTGCTATCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAAtt                    <  2:216851/91‑1 (MQ=38)
                                                                     aCAGGGTAGCGATTGCTATCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAAtt                    >  1:777805/1‑91 (MQ=38)
                                                                     aCAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAAtt                    <  2:795115/91‑1 (MQ=255)
                                                                     aCAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAAtt                    <  1:56530/91‑1 (MQ=255)
                                                                      cAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTa                   >  1:1912169/1‑91 (MQ=255)
                                                                       aGGGTAGCGATTGCTATCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTAc                  >  1:702428/1‑91 (MQ=39)
                                                                          aTAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTg               <  2:900672/91‑1 (MQ=255)
                                                                          aTAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTg               <  2:1250370/91‑1 (MQ=255)
                                                                           tAGCGATTGCTATCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGc              >  2:1206470/1‑91 (MQ=255)
                                                                           tAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGc              >  2:2777913/1‑91 (MQ=255)
                                                                           tAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGc              <  1:991289/91‑1 (MQ=255)
                                                                           tAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGc              <  2:1333074/91‑1 (MQ=255)
                                                                            aGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCa             >  2:1354336/1‑91 (MQ=255)
                                                                            aGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCa             >  1:3394266/1‑91 (MQ=255)
                                                                             gCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCAc            >  2:342903/1‑91 (MQ=255)
                                                                              cAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCAcc           >  1:3819373/1‑91 (MQ=255)
                                                                              cAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCAcc           <  2:867168/91‑1 (MQ=255)
                                                                                aTCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACCtt         >  2:3522453/1‑91 (MQ=255)
                                                                                aTCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACCtt         >  2:3145819/1‑91 (MQ=255)
                                                                                aTCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACCtt         >  2:2395018/1‑91 (MQ=255)
                                                                                   gCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACCTTCgg      >  1:950929/1‑91 (MQ=255)
                                                                                    cTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACCTTCGGc     <  1:3145355/91‑1 (MQ=255)
                                                                                     tACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACCTTCGGCg    >  1:526880/1‑91 (MQ=255)
                                                                                       cccTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACCTTCGGCGAt  >  1:3576331/1‑91 (MQ=255)
                                                                                       cccTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACCTTCGGCGAt  >  2:3508976/1‑91 (MQ=255)
                                                                                       cccTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACCTTCGGCGAt  >  2:1249049/1‑91 (MQ=255)
                                                                                       cccTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACCTTCGGCGAt  >  2:2052340/1‑91 (MQ=255)
                                                                                       cccTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACCTTCGGCGAt  >  2:1540759/1‑91 (MQ=255)
                                                                                       |                                                                                          
CATAATTTTGATATAAAGGTGTTTGTTTTAAATCACTTGACATTGAAATACCCCCTTAAAATTATAAAAACAGGATAGCAATCGCTACCCTGTTTATAAATCACCAATTCAGGAACGCGTTACAGTATTATCCTTTTGCCTGAGAGTTTCGTAATAAATTACTTGCACCTTCGGCGAT  >  USA300TCH1516_ALE/1653253‑1653430

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: