Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,929,573 C→A 5.9% R31L (CGC→CTC)  USA300TCH1516_01793 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,929,5730CA 5.9% 1004.3 / 26.4 289R31L (CGC→CTC) USA300TCH1516_01793hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (130/140);  new base A (14/4);  total (144/145)
Fisher's exact test for biased strand distribution p-value = 2.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 306 total aligned reads displayed.

GAATGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTTCTCTTGTCA  >  USA300TCH1516_ALE/1929483‑1929663
                                                                                          |                                                                                          
gAATGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTcgc                                                                                            <  1:2367542/91‑1 (MQ=255)
gAATGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTcgc                                                                                            >  1:2088053/1‑91 (MQ=255)
  aTGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACCCGCGt                                                                                          >  1:1625342/1‑91 (MQ=255)
   tGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGta                                                                                         <  2:2609997/91‑1 (MQ=255)
    gTATCCGTGCTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATTATACGCGAGtat                                                                                        <  1:3553495/91‑1 (MQ=255)
      aCCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATaa                                                                                      <  2:753362/91‑1 (MQ=255)
      aCCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATaa                                                                                      >  1:2108181/1‑91 (MQ=255)
       cccGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAAt                                                                                     <  1:3713162/91‑1 (MQ=255)
       cccGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAAt                                                                                     <  2:2948308/91‑1 (MQ=255)
        ccGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAAtt                                                                                    <  1:717098/91‑1 (MQ=255)
         cGTGTTTCGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGGGTATAAttt                                                                                   <  1:1971275/91‑1 (MQ=255)
          gtgtTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTc                                                                                  >  1:1836514/1‑91 (MQ=255)
          gtgtTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTc                                                                                  >  1:1930698/1‑91 (MQ=255)
           tgtTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCa                                                                                 >  2:1420230/1‑91 (MQ=255)
            gtTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCAt                                                                                <  2:3413214/91‑1 (MQ=255)
                 gATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTa                                                                           >  2:454261/1‑91 (MQ=255)
                   tGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAgg                                                                         >  2:2423419/1‑91 (MQ=255)
                    gTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGc                                                                        >  1:2554562/1‑91 (MQ=255)
                      tGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCtt                                                                      <  1:1444963/91‑1 (MQ=255)
                      tGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCtt                                                                      <  1:855938/91‑1 (MQ=255)
                      tGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCtt                                                                      <  1:520558/91‑1 (MQ=255)
                       gCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTa                                                                     >  2:3012737/1‑91 (MQ=255)
                       gCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTa                                                                     >  2:1036233/1‑91 (MQ=255)
                        cTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTAc                                                                    <  1:216042/91‑1 (MQ=255)
                        cTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTAc                                                                    <  2:1579438/91‑1 (MQ=255)
                        cTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTAc                                                                    >  2:2843594/1‑91 (MQ=255)
                        cTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTAc                                                                    <  2:2864589/91‑1 (MQ=255)
                          ttATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCa                                                                  >  1:1908874/1‑91 (MQ=255)
                           tATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCAt                                                                 >  2:1159836/1‑91 (MQ=255)
                           tATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCAt                                                                 >  1:340571/1‑91 (MQ=255)
                             tCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCAttt                                                               <  2:3100252/91‑1 (MQ=255)
                              ccATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCAtttt                                                              >  2:902521/1‑91 (MQ=255)
                               cATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCAttttt                                                             <  2:12969/91‑1 (MQ=255)
                                 ttttCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTtat                                                           <  2:1864437/91‑1 (MQ=255)
                                   ttCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTtatat                                                         <  1:133372/91‑1 (MQ=255)
                                   ttCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTtatat                                                         <  2:1095681/91‑1 (MQ=255)
                                      cTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAAt                                                      <  2:2187439/91‑1 (MQ=255)
                                      cTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAAt                                                      <  1:454261/91‑1 (MQ=255)
                                        aCATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCt                                                    >  1:1347080/1‑91 (MQ=255)
                                        aCATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGAATAATTTCATTCCTAGGCTTACCATTTTTATATAATCt                                                    <  2:2277653/91‑1 (MQ=255)
                                         cATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTa                                                   >  1:2905200/1‑91 (MQ=255)
                                         cATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTa                                                   >  1:1680537/1‑91 (MQ=255)
                                           ttttATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAAc                                                 <  1:41022/91‑1 (MQ=255)
                                              tataACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCAt                                              <  1:3483734/91‑1 (MQ=255)
                                               ataACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCAtt                                             >  2:3602654/1‑91 (MQ=255)
                                                taACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCAttt                                            >  1:437706/1‑91 (MQ=255)
                                                taACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCAttt                                            >  1:59714/1‑91 (MQ=255)
                                                 aaCCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTg                                           <  2:590331/91‑1 (MQ=255)
                                                  aCCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAACTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGt                                          >  2:2566031/1‑91 (MQ=255)
                                                   ccGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTa                                         <  2:139045/91‑1 (MQ=255)
                                                     gCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAAc                                       >  2:385057/1‑91 (MQ=255)
                                                     gCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAAc                                       >  2:1622889/1‑91 (MQ=255)
                                                     gCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAAc                                       <  2:2528639/91‑1 (MQ=255)
                                                     gCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAAc                                       >  1:3327480/1‑91 (MQ=255)
                                                       cATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGCAACtt                                     >  2:3647791/1‑91 (MQ=255)
                                                        aTTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTa                                    <  2:2488987/91‑1 (MQ=255)
                                                        aTTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTa                                    >  1:2411544/1‑91 (MQ=255)
                                                          ttACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTaaa                                  <  1:3281854/91‑1 (MQ=255)
                                                          ttACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTaaa                                  >  1:2439721/1‑91 (MQ=255)
                                                          ttACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTaaa                                  <  2:3221001/91‑1 (MQ=255)
                                                            aCAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAAct                                >  2:726830/1‑91 (MQ=255)
                                                            aCAAACATCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAAct                                >  1:3074059/1‑91 (MQ=255)
                                                             cAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAActc                               >  1:2952470/1‑91 (MQ=255)
                                                              aaaCGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAActct                              <  1:1591910/91‑1 (MQ=255)
                                                              aaaCGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAActct                              <  1:3760238/91‑1 (MQ=255)
                                                              aaaCGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAActct                              <  2:584769/91‑1 (MQ=255)
                                                               aaCGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTg                             >  1:103081/1‑91 (MQ=255)
                                                               aaCGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTg                             <  2:444229/91‑1 (MQ=255)
                                                                 cGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATGTAATCTAACCATTTGTAACTTAAACTCTGaa                           >  1:3370760/1‑91 (MQ=255)
                                                                  gTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAAc                          <  2:538551/91‑1 (MQ=255)
                                                                    cGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTa                        <  1:3212056/91‑1 (MQ=255)
                                                                    cGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTa                        <  1:3840216/91‑1 (MQ=255)
                                                                     gAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTaa                       >  2:2215549/1‑91 (MQ=255)
                                                                     gAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTaa                       <  1:2733334/91‑1 (MQ=255)
                                                                       aGGTGTGAAATCATTCTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAAt                     <  2:2815314/91‑1 (MQ=255)
                                                                       aGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAAt                     >  2:3283331/1‑91 (MQ=255)
                                                                        ggTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATg                    <  2:2021357/91‑1 (MQ=255)
                                                                        ggTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATg                    <  1:1291925/91‑1 (MQ=255)
                                                                        ggTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATg                    <  1:3609804/91‑1 (MQ=255)
                                                                         gtgtGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGa                   <  1:635543/91‑1 (MQ=255)
                                                                         gtgtGAAATCATACCCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGa                   >  2:3359629/1‑91 (MQ=255)
                                                                           gtgAAATCATACGCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtc                 >  1:3029118/1‑91 (MQ=255)
                                                                             gAAATCATACTCGCGTATAATTTCATTCGTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtctt               <  2:1539652/91‑1 (MQ=255)
                                                                             gAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtctt               <  1:2191648/91‑1 (MQ=255)
                                                                              aaaTTATACGCGAGTATGATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtcttc              >  2:2704978/1‑91 (MQ=14)
                                                                               aaTCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtcttct             <  1:1713736/91‑1 (MQ=255)
                                                                                attaTACGCGAGTATGATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATcttctt            >  1:1375245/4‑91 (MQ=14)
                                                                                aTCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATcttctt            <  1:2825327/91‑1 (MQ=255)
                                                                                aTCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATcttctt            <  1:2571581/91‑1 (MQ=255)
                                                                                aTCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATcttctt            <  1:2664349/91‑1 (MQ=255)
                                                                                 tCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCttcttt           <  2:414332/91‑1 (MQ=255)
                                                                                   atatTCGCTTAGAATGTTATTCCTATGGTTACCATTGTTATATAATCTAACCATTTGTAACTTAAACGCTGAACTAAATGATCtgctttct         <  2:182089/91‑1 (MQ=14)
                                                                                   aTACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCttctttct         <  2:1731484/91‑1 (MQ=255)
                                                                                   aTACGCGAGTATGATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCttctttct         >  1:2341937/1‑91 (MQ=14)
                                                                                    tACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTtctc        >  1:2120728/1‑91 (MQ=255)
                                                                                    tACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTtctc        <  1:1185949/91‑1 (MQ=255)
                                                                                       tCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTTCTCTTg     >  1:2515307/1‑91 (MQ=255)
                                                                                         gcgTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTTCTCTTGTc   >  1:1689235/1‑91 (MQ=255)
                                                                                         gcgTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTTCTCTTGTc   >  1:2476472/1‑91 (MQ=255)
                                                                                          cgTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTTCTCTTGTCa  <  1:79358/91‑1 (MQ=255)
                                                                                          |                                                                                          
GAATGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTTCTCTTGTCA  >  USA300TCH1516_ALE/1929483‑1929663

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