Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 1,963,768 C→G V365L (GTA→CTA)  hemY ← Protoporphyrinogen oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,963,7680CG98.6% 233.9 / NA 70V365L (GTA→CTA) hemYProtoporphyrinogen oxidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (33/36);  minor base A (1/0);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCT  >  USA300TCH1516_ALE/1963634‑1963904
                                                                                                                                      |                                                                                                                                        
tttCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCCCCTGGTTTCCCTAGAtcgc                                                                                                                                     >  2:142282/1‑137 (MQ=255)
 ttCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGCTACGc                                                                                                                                    >  1:963471/1‑140 (MQ=255)
 ttCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGc                                                                                                                                    >  1:574002/1‑140 (MQ=255)
 ttCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATAATCAACCACCGCAGCACATGGTTTACCTAGATACGc                                                                                                                                    >  2:206706/1‑140 (MQ=255)
 ttCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGAACATCAACCCCAGTGGCACCTGGTTTACCTCAATACGc                                                                                                                                    >  2:206460/1‑140 (MQ=255)
     gCAAAGGATTGACAATTGTAAATTCAGGATCACCTTTCAATGTCATCATTTGACTTAAATCTCTCCGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCAGCTGGTTTACCTAGATACGCACGa                                                                                                                                <  2:241845/140‑1 (MQ=255)
     gCAAACGATTTACAATTGTAAATTCAGCAGCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATGCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGa                                                                                                                                <  2:467473/140‑1 (MQ=255)
      cAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGaa                                                                                                                               <  1:722079/140‑1 (MQ=255)
      cAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGaa                                                                                                                               <  1:462466/140‑1 (MQ=255)
        aaCGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCCTCATTTGACTTAAATCTCTACGTACACTCGTTACTAATTCATTATCTGTATGATCATCAACCACAGTAGCACCTGGTTTACCTCGCTACGCACGAATc                                                                                                                             <  2:506342/140‑1 (MQ=255)
        aaCGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTCACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATc                                                                                                                             <  2:522830/140‑1 (MQ=255)
               gACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAacct                                                                                                                      <  2:664034/140‑1 (MQ=255)
                aCAATTGTAGATTCAGGACCACCTTTAAATGTCAGCATTTTACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAacctt                                                                                                                     <  2:668625/140‑1 (MQ=255)
                  aaTTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAaccttac                                                                                                                   <  2:438541/140‑1 (MQ=255)
                    ttGTAAATTCACTATCACCTTTAAATGCCATCATTTGACTTAAATCCCTACGTACAATCGATACTCATTCATTATCTGTATGATCATCCACCACAGTATCACCTGGTTTACCTAGATACCCACGAATCAAAaccttacct                                                                                                                 <  2:189863/140‑1 (MQ=255)
                         aaTTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCgg                                                                                                            >  2:266274/1‑140 (MQ=255)
                           ttCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCGCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGgtg                                                                                                          <  2:989063/140‑1 (MQ=255)
                               ggATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTgtagt                                                                                                      >  2:60470/1‑140 (MQ=255)
                                      cTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGc                                                                                               >  2:143142/1‑140 (MQ=255)
                                         tataggTCATCATTTGACTTACATCCCTACGGACAACCGATACTACTTCATGATCTGTATGATCATCAACCACAGTATCACCTGTTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCAt                                                                                            <  1:276134/135‑1 (MQ=255)
                                         tAACTGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCAt                                                                                            <  1:186958/140‑1 (MQ=255)
                                          aaaTGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATCATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTTTAGCAAATGGCCaat                                                                                           >  2:561179/1‑138 (MQ=255)
                                             tGTCATCATTTCACTTAAATCTCCACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCAttttt                                                                                        <  2:513788/140‑1 (MQ=255)
                                             tGTCACCGTGTGACGGAAATCTCTACGCACATCAGATACGAATTCACTAGCGGCACAAACAGCAACCACAGGGTCGCCGGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGCGTAGTAAACGGCCAtcttt                                                                                        >  2:708564/1‑140 (MQ=255)
                                                   catTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATg                                                                                  >  1:613641/1‑140 (MQ=255)
                                                   catTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATAATCACACACAGTATCACCTGGTTTACCTAGATACTAAAGAATAAAAACCTTACCTACCGGAGTAGTACATGGCCATATTTTCGATg                                                                                  >  1:462819/1‑140 (MQ=255)
                                                     tttGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTc                                                                                >  2:779857/1‑140 (MQ=255)
                                                      ttGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTcc                                                                               <  2:825038/140‑1 (MQ=255)
                                                        gACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGCTGTCCaa                                                                             >  1:433792/1‑140 (MQ=255)
                                                            tcatctcctACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTAc                                                                         <  2:423409/133‑1 (MQ=255)
                                                            tAACTCTCTACGTACAAATGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGCTACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTAc                                                                         <  1:350978/140‑1 (MQ=255)
                                                            tAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCGGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTAc                                                                         >  2:834633/1‑140 (MQ=255)
                                                            tAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCCTTTTTTCGATTTCCAAGTAc                                                                         >  1:575396/1‑140 (MQ=255)
                                                            tAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTAc                                                                         <  1:966486/140‑1 (MQ=255)
                                                            tAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTAc                                                                         >  1:809754/1‑140 (MQ=255)
                                                             aaaTCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACa                                                                        <  1:266274/140‑1 (MQ=255)
                                                                   ctACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGt                                                                  <  2:618276/140‑1 (MQ=255)
                                                                   ctACGTACAACCGCTACCAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGt                                                                  <  1:356770/140‑1 (MQ=255)
                                                                       gTACAATCGATACTAATTCATTATCTGTATTATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAATCCTTACCTTCCGGTGTAGTGAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATg                                                              <  2:544374/140‑1 (MQ=255)
                                                                          cAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTTTCTGTCCTTGTACATGCTGTAATgagt                                                           >  2:699819/1‑137 (MQ=255)
                                                                          cAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATAACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGCAGTAGTACATGGCCATTTTTTCGATGTACAAGTAACTGCGGGAgtgtat                                                           >  1:116274/1‑138 (MQ=255)
                                                                           aaTCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTCCCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTg                                                          <  2:884211/140‑1 (MQ=255)
                                                                            aTCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGt                                                         >  2:293524/1‑140 (MQ=255)
                                                                                  aCTAATTCATTATCTGTATGATCATCAGCCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACt                                                   <  2:697344/140‑1 (MQ=255)
                                                                                  aCTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACt                                                   <  2:372648/140‑1 (MQ=255)
                                                                                  aCTAATTCATTATCCGTATGATCATCAACCACAGTATCACCTGGTTTACCCAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACt                                                   <  2:759225/140‑1 (MQ=255)
                                                                                    tAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCg                                                 <  2:844632/140‑1 (MQ=255)
                                                                                             aTCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCaa                                        <  1:975795/140‑1 (MQ=255)
                                                                                              tCTGTATGATCATCGACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAAt                                       >  1:554734/1‑140 (MQ=255)
                                                                                                  tATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACg                                   >  2:780508/1‑140 (MQ=255)
                                                                                                   aTGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGa                                  >  1:804912/1‑140 (MQ=255)
                                                                                                    tGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGaa                                 >  2:507072/1‑140 (MQ=255)
                                                                                                         atcaACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAg                            <  1:293524/140‑1 (MQ=255)
                                                                                                           caACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTa                          >  1:490311/1‑140 (MQ=255)
                                                                                                           caACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGACAGTa                          >  1:99661/1‑140 (MQ=255)
                                                                                                             aCCACAGTATCACCTGGTTTACCTAGATACGCACGAACCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTGCTATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTAcc                        <  1:384376/140‑1 (MQ=255)
                                                                                                               cacaGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTTTTTTTTTTTGTACATGCGGTAATGTCTGAATCACTCCTCTTCGCTTATACGATACCAGTACtat                      >  2:296030/1‑140 (MQ=255)
                                                                                                               cacaGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCCGTACcat                      >  1:456240/1‑140 (MQ=255)
                                                                                                               cacaGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGGCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCACTTACGAAGCCAATACgat                      >  1:13907/1‑140 (MQ=255)
                                                                                                                 caGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCACTTACGAAGCCAGTACcgtca                    >  2:674374/1‑140 (MQ=255)
                                                                                                                   tcaTCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTCCCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATg                  <  2:528280/138‑1 (MQ=255)
                                                                                                                      tCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGt               >  2:985094/1‑140 (MQ=255)
                                                                                                                         ccTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTAtt            <  2:580743/140‑1 (MQ=255)
                                                                                                                          cTGTTTTACTTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAAGGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTAttt           <  2:789326/140‑1 (MQ=255)
                                                                                                                               ttACCTAGATACGCACGAAACAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTGTATTTTTTTTTACATGCGGTAAGGTCAGGATCACTCGTTCTCGCAATAAAGAAGCAAGTTCCATCATGGAGATTTTCaag      >  2:560364/1‑139 (MQ=255)
                                                                                                                                tACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATGACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATg     >  1:882634/1‑140 (MQ=255)
                                                                                                                                 aCCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAAGGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGt    <  2:78720/140‑1 (MQ=255)
                                                                                                                                  ccTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTc   >  1:661552/1‑140 (MQ=255)
                                                                                                                                   cTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCCAGTACATGCGGTAATGTCTGTATCACTCGTTCTACCAAGTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCt  <  2:917401/140‑1 (MQ=255)
                                                                                                                                   cTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGAACGCTGCATCATGGGTATTTTCAATGTCt  <  2:712932/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                        
TTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCT  >  USA300TCH1516_ALE/1963634‑1963904

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