Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,542,910 C→A 5.4% intergenic (‑412/+19) USA300HOU_RS07525 ← / ← hup hypothetical protein/DNA‑binding protein HU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,542,9100CA 5.4% 1016.8 / 29.5 286intergenic (‑412/+19)USA300HOU_RS07525/huphypothetical protein/DNA‑binding protein HU
Reads supporting (aligned to +/- strand):  ref base C (116/154);  new base A (5/11);  total (121/165)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 291 total aligned reads displayed.

AATGAATGTTCCAAAGATTAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGCTGGAACTTTACTTGCTGGGATATCAATTT  >  USA300TCH1516_ALE/1542820‑1542999
                                                                                          |                                                                                         
aatgaatgTTCCAAAGATTAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTc                                                                                           >  2:758180/1‑91 (MQ=255)
 atgaatgTTCCAAAGATTAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCa                                                                                          <  1:1784588/91‑1 (MQ=255)
 atgaatgTTCCAAAGATTAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCa                                                                                          <  2:3752624/91‑1 (MQ=255)
 atgaatgTTCCAAAGAGTAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCa                                                                                          <  2:3113240/91‑1 (MQ=255)
  tgaatgTTCCAAAGATTAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAg                                                                                         >  1:771251/1‑91 (MQ=255)
   gaatgTTCCAAAGATTAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAgg                                                                                        <  1:2140163/91‑1 (MQ=255)
   gaatgTTCCAAAGATTAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAgg                                                                                        <  2:1247022/91‑1 (MQ=255)
   gaatgTTCCAAAGATTAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAgg                                                                                        <  1:1344585/91‑1 (MQ=255)
   gaatgTTCCAAAGATTAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAgg                                                                                        <  1:948464/91‑1 (MQ=255)
     atgTTCCAAAGATTAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGc                                                                                      >  2:1151756/1‑91 (MQ=255)
     atgTTCCAAAGATTAAAAACTAAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGc                                                                                      >  1:929907/1‑91 (MQ=255)
      tgTTCCAAAGATTAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCt                                                                                     <  1:3698592/91‑1 (MQ=255)
      tgTTCCAAAGATTAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTGAATAAGGGCt                                                                                     <  1:2200172/91‑1 (MQ=255)
      tgTTCCAAAGATTAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTGAATAAGGGCt                                                                                     <  1:2274225/91‑1 (MQ=255)
       gTTCCAAAGATTAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCtt                                                                                    >  1:3800788/1‑91 (MQ=255)
        ttCCAAAGATTAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCttt                                                                                   >  2:649245/1‑91 (MQ=255)
        ttCCAAAGATTAAAAACTTAGAAGAACATTACAATTTATGTACAGGGTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCttt                                                                                   <  2:818928/91‑1 (MQ=255)
         tCCAAAGATTAACAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCtttt                                                                                  >  1:3038242/1‑91 (MQ=255)
         tCCAAAGATTAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCtttt                                                                                  <  1:800575/91‑1 (MQ=255)
          ccAAAGATTAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTGAATAAGGGCttttt                                                                                 <  1:425186/91‑1 (MQ=38)
          ccAAAGATTAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTGAATAAGGGCttttt                                                                                 <  1:784545/91‑1 (MQ=38)
           cAAAGATTAAAATCTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTGAATAAGGGCTTTTTa                                                                                <  2:3291389/91‑1 (MQ=38)
           cAAAGATTAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTa                                                                                <  2:3599569/91‑1 (MQ=255)
           cAAAGATTAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTGAATAAGGGCTTTTTa                                                                                <  1:2011464/91‑1 (MQ=38)
            aaaGATTAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTaa                                                                               <  2:1454112/91‑1 (MQ=255)
              aGATTAACAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGt                                                                             <  2:3503805/91‑1 (MQ=255)
              aGATTAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGt                                                                             <  2:2544626/91‑1 (MQ=255)
              aGATTAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGt                                                                             <  1:1474360/91‑1 (MQ=255)
               gATTAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTa                                                                            <  2:2792783/91‑1 (MQ=255)
                 ttAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTaaa                                                                          <  1:2785474/91‑1 (MQ=255)
                 ttAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTaaa                                                                          <  2:1060799/91‑1 (MQ=255)
                  tAAAAACTTAGAAGAAGATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTaaaa                                                                         <  2:3777033/91‑1 (MQ=255)
                  tAAAAACTTAGAAGAACGTTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTaaaa                                                                         <  1:3382301/91‑1 (MQ=255)
                    aaaaCTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATACAAAGCCCTTATTCAGGGCTTTTTAAGTAAAatt                                                                       >  1:1546305/1‑91 (MQ=255)
                    aaaaCTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAatt                                                                       <  1:2404179/91‑1 (MQ=255)
                    aaaaCTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAatt                                                                       <  1:1573741/91‑1 (MQ=255)
                    aaaaCTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAatt                                                                       <  2:2913298/91‑1 (MQ=255)
                         ttAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTAtttt                                                                  <  2:2256922/91‑1 (MQ=255)
                           agaagaACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTAc                                                                <  2:1619501/91‑1 (MQ=255)
                           agaagaACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTAc                                                                <  1:1206891/91‑1 (MQ=255)
                           agaagaACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTAc                                                                <  1:2833994/91‑1 (MQ=255)
                            ggagaaCATGACAGTGTGTGTACAGGTGGTAAGTATTAAAAGCAAAATAAAAAGCCCTTAGTCAGGGCTTTTTAAGTAAAATAATTTTAct                                                               <  2:2465114/89‑2 (MQ=255)
                            gaagaaCAGTACAATTGATGTACATGTTGTAAGTATTAAAAGCAAAATAAAAAGTCCTTATTCATGGCTTTTTAAGTAAAATTATTTTACa                                                               <  2:2380389/91‑1 (MQ=255)
                             aagaaCATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAATAGACCTTATTCCGGGCATTTTAAGTAAAATTCTTTTACaa                                                              >  2:2689548/1‑90 (MQ=255)
                               gaaCATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCa                                                            >  1:3106359/1‑91 (MQ=255)
                                aaCGTTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCAt                                                           <  1:3066748/91‑1 (MQ=255)
                                aaCATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCAt                                                           <  1:1099291/91‑1 (MQ=255)
                                  cATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCt                                                         <  1:3278663/91‑1 (MQ=255)
                                   aTTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCtt                                                        <  1:608723/91‑1 (MQ=255)
                                     tACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCtttt                                                      >  2:1935801/1‑91 (MQ=255)
                                      aCAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTa                                                     <  2:2103665/91‑1 (MQ=255)
                                      aCAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTa                                                     <  1:570704/91‑1 (MQ=255)
                                      aCAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTa                                                     <  1:1751547/91‑1 (MQ=255)
                                      aCAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTa                                                     >  1:2508927/1‑91 (MQ=255)
                                        aaTTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCTGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAAt                                                   >  2:463774/1‑91 (MQ=255)
                                        aaTTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAAt                                                   <  1:714663/91‑1 (MQ=255)
                                          tttATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGc                                                 <  1:3907935/91‑1 (MQ=255)
                                           ttATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCt                                                >  1:1691882/1‑91 (MQ=255)
                                             aTGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCttt                                              <  2:874038/91‑1 (MQ=255)
                                             aTGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCttt                                              >  1:1414739/1‑91 (MQ=255)
                                               gTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTAc                                            >  1:1619080/1‑91 (MQ=255)
                                                  cAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAg                                         >  2:299614/1‑91 (MQ=255)
                                                     gTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCtt                                      <  1:658880/91‑1 (MQ=255)
                                                       tGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTg                                    >  1:3602043/1‑91 (MQ=255)
                                                        gTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGa                                   <  1:1004877/91‑1 (MQ=255)
                                                          aaGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAAt                                 <  1:1035177/91‑1 (MQ=255)
                                                          aaGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAAt                                 >  1:232745/1‑91 (MQ=255)
                                                          aaGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAAt                                 >  2:1723678/1‑91 (MQ=255)
                                                          aaGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAAt                                 >  1:3482262/1‑91 (MQ=255)
                                                           aGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATg                                <  1:2624796/91‑1 (MQ=255)
                                                            gTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGc                               <  2:550562/91‑1 (MQ=255)
                                                            gTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGc                               >  2:2508449/1‑91 (MQ=255)
                                                            gTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGc                               <  2:266401/91‑1 (MQ=255)
                                                              aTTAAAAGCAAAATCAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGCTg                             >  2:3357846/1‑91 (MQ=255)
                                                              aTTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGCTg                             >  2:3860778/1‑91 (MQ=255)
                                                              aTTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGCTg                             >  1:910781/1‑91 (MQ=255)
                                                               ttAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGCTgg                            >  2:1863877/1‑91 (MQ=255)
                                                               ttAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGCTgg                            <  2:2340383/91‑1 (MQ=255)
                                                                tAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCATTCAAAGCTggt                           <  2:1780170/91‑2 (MQ=255)
                                                                  aaaGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCATTCAAAGCTggtaa                         <  1:2073206/91‑4 (MQ=255)
                                                                     gCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGCTGGAACttt                      <  2:3164475/91‑1 (MQ=255)
                                                                     gCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGCTGGAACttt                      >  2:1285050/1‑91 (MQ=255)
                                                                      cAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGCTGGCACTTTa                     >  2:3735521/1‑91 (MQ=255)
                                                                            aaaaaGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGCTGGAACTTTACTTGCt               >  2:2178280/1‑91 (MQ=255)
                                                                            aaaaaGCCCTGAATAAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGCTGGAACTTTACTTGCt               >  2:139147/1‑91 (MQ=38)
                                                                             aaaaGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGCTGGAACTTTACTTGCTg              <  2:3055927/91‑1 (MQ=255)
                                                                             aaaaGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGCTGGAACTTTACTTGCTg              >  1:2705224/1‑91 (MQ=255)
                                                                               aaGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTTCATGCTGGAACTTTACTTGCTggg            >  2:2835600/1‑91 (MQ=255)
                                                                                aGCCCTTATTCAGGGCTTTTTACGTAAAATTAATTTACAGCATCTCTTAATGCTTTACCCGCTTTGAATGCTGGAACTTTCCTTGCTGGGa           >  2:536388/1‑91 (MQ=255)
                                                                                aGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGCTGGAACTTTACTTGCTGGGa           >  2:930922/1‑91 (MQ=255)
                                                                                aGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGCTGGAACTTTACTTGCTGGGa           <  1:1044206/91‑1 (MQ=255)
                                                                                aGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGCTGGAACTTTACTTGCTGGGa           <  1:894344/91‑1 (MQ=255)
                                                                                 gCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGCTGGAACTTTACTTGCTGGGat          >  2:1336638/1‑91 (MQ=255)
                                                                                   ccTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGCTGGAACTTTACTTGCTGGGatat        <  1:2081712/91‑1 (MQ=255)
                                                                                   ccTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGCTGGAACTTTACTTGCTGGGatat        <  2:709400/91‑1 (MQ=255)
                                                                                    cTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGCTGGAACTTTACTTGCTGGGATATc       >  1:209645/1‑91 (MQ=255)
                                                                                      tATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGCTGGAACTTTACTTGCTGGGATATCaa     <  2:148547/91‑1 (MQ=255)
                                                                                       aTTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGCTGGAACTTTACTTGCTGGGATATCAAt    <  1:2862579/91‑1 (MQ=255)
                                                                                       aTTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGCTGGAACTTTACTTGCTGGGATATCAAt    >  1:2074649/1‑91 (MQ=255)
                                                                                        ttCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGCTGGAACTTTACTTGCTGGGATATCAAtt   <  1:400542/91‑1 (MQ=255)
                                                                                          |                                                                                         
AATGAATGTTCCAAAGATTAAAAACTTAGAAGAACATTACAATTTATGTACAGGTTGTAAGTATTAAAAGCAAAATAAAAAGCCCTTATTCAGGGCTTTTTAAGTAAAATTATTTTACAGCATCTTTTAATGCTTTACCAGCTTTGAATGCTGGAACTTTACTTGCTGGGATATCAATTT  >  USA300TCH1516_ALE/1542820‑1542999

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