Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 690,101 A→T 10.5% intergenic (‑30/+5) USA300HOU_RS03320 ← / ← USA300HOU_RS03325 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE690,1010AT10.5% 855.1 / 67.2 282intergenic (‑30/+5)USA300HOU_RS03320/USA300HOU_RS03325hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (122/128);  new base T (11/19);  total (134/148)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 287 total aligned reads displayed.

CGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAA  >  USA300TCH1516_ALE/690011‑690191
                                                                                          |                                                                                          
cGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATa                                                                                            >  2:2236333/1‑91 (MQ=255)
cGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATa                                                                                            >  1:3039440/1‑91 (MQ=255)
cGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATa                                                                                            >  2:3524391/1‑91 (MQ=255)
 gCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAg                                                                                           <  2:2330602/91‑1 (MQ=255)
 gCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAg                                                                                           <  1:3471609/91‑1 (MQ=255)
  ccAGTTTTCGTTAATTACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCAAGTCATCCCTCCAATAAAAATAACCCCAATAgg                                                                                          >  1:3365378/1‑91 (MQ=255)
   cAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAggg                                                                                         <  2:563460/91‑1 (MQ=255)
   cAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAggg                                                                                         <  1:3776523/91‑1 (MQ=255)
     gTTTTCGTTAATAACAAATTATCTCTTAAACCTGCGTACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGc                                                                                       <  2:622849/91‑1 (MQ=255)
     gTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGc                                                                                       <  1:3199897/91‑1 (MQ=255)
      ttttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTggggtt                                                                                      <  1:2781827/91‑3 (MQ=255)
      ttttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTggggtt                                                                                      <  1:3500096/91‑3 (MQ=255)
      ttttCGTCAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTggggtt                                                                                      <  1:1974929/91‑3 (MQ=255)
       tttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCtt                                                                                     >  2:120054/1‑90 (MQ=255)
        ttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTAt                                                                                    <  2:2754581/91‑1 (MQ=255)
          cGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTAttt                                                                                  >  2:2853952/1‑91 (MQ=255)
          cGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTAttt                                                                                  >  2:3002503/1‑91 (MQ=255)
          cGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTAttt                                                                                  >  1:3254562/1‑91 (MQ=255)
            ttAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTaa                                                                                >  2:1916993/1‑91 (MQ=255)
            ttAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTaa                                                                                >  2:1570419/1‑91 (MQ=255)
               ataaCAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATg                                                                             <  1:163103/91‑1 (MQ=255)
               ataaCAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATg                                                                             >  2:226036/1‑91 (MQ=255)
                taaCAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGt                                                                            >  2:637359/1‑91 (MQ=255)
                    aaaTTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAtt                                                                        <  1:738028/91‑1 (MQ=255)
                     aaTTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAttt                                                                       >  1:3135007/1‑91 (MQ=255)
                        tATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAtttttt                                                                    <  1:1561169/91‑1 (MQ=255)
                          tctcTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTa                                                                  >  1:1016683/1‑91 (MQ=255)
                           ctctTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAGTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTAt                                                                 <  2:1005228/91‑1 (MQ=255)
                           ctctTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTAt                                                                 <  1:114743/91‑1 (MQ=255)
                           ctctTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTAt                                                                 <  1:370470/91‑1 (MQ=255)
                           ctctTAAACCTGCATACGGTGCTTTTAAAGTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTAt                                                                 <  2:1386274/91‑1 (MQ=255)
                            tctTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATc                                                                <  1:3306417/91‑1 (MQ=255)
                            tctTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATc                                                                <  1:966325/91‑1 (MQ=255)
                            tctTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATc                                                                <  1:509005/91‑1 (MQ=255)
                            tctTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATc                                                                <  1:441295/91‑1 (MQ=255)
                            tctTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATc                                                                <  1:3939449/91‑1 (MQ=255)
                             ctTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCa                                                               >  2:3116862/1‑91 (MQ=255)
                              ttAAACCTGCATACGGTGCTTTTAAATTCATGTCATTCCTCCAGTAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATcgt                                                              <  2:2446301/91‑3 (MQ=255)
                              ttAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCAt                                                              >  2:1831387/1‑91 (MQ=255)
                                aaaCCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGa                                                            >  1:3022546/1‑91 (MQ=255)
                                   ccTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGatat                                                         <  1:1479183/91‑1 (MQ=255)
                                   ccTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGatat                                                         <  2:3704544/91‑1 (MQ=255)
                                     tGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTc                                                       <  2:3244201/91‑1 (MQ=255)
                                     tGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTc                                                       >  2:2799798/1‑91 (MQ=255)
                                      gCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCt                                                      >  1:1321591/1‑91 (MQ=255)
                                      gCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCt                                                      >  1:937622/1‑91 (MQ=255)
                                      gCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCt                                                      >  1:1900173/1‑91 (MQ=255)
                                       cATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTa                                                     >  2:380048/1‑91 (MQ=255)
                                         tACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATCGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTAct                                                   >  1:3339696/1‑91 (MQ=255)
                                           cGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTActtc                                                 <  2:2535915/91‑1 (MQ=255)
                                           cGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTActtc                                                 >  2:3875357/1‑91 (MQ=255)
                                            ggTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTActtct                                                <  1:3549177/91‑1 (MQ=255)
                                               gCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTtca                                             >  2:3879841/1‑91 (MQ=255)
                                                  tttAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTtcatca                                          <  2:1288646/91‑1 (MQ=255)
                                                   ttAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTcatcat                                         >  1:2053666/1‑91 (MQ=255)
                                                   ttAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTcatcat                                         >  2:2156898/1‑91 (MQ=255)
                                                    tAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATCCTACTTCTTCATCATa                                        >  2:771505/1‑91 (MQ=255)
                                                     aaaTTCATGTCATCCCTCCAATAAAAATACCCCCAATAGGGGCTATTTAAATGTTATTTTTTTATAATGATATTCTACTTCTGCATCATaa                                       >  2:3433627/1‑91 (MQ=255)
                                                     aaaTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATaa                                       <  2:403610/91‑1 (MQ=255)
                                                      aaTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAAc                                      <  1:1751850/91‑1 (MQ=255)
                                                      aaTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAAc                                      <  2:3222435/91‑1 (MQ=255)
                                                       aTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACt                                     <  1:2910849/91‑1 (MQ=255)
                                                           aTGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAAt                                 >  1:2615146/1‑91 (MQ=255)
                                                           aTGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAAt                                 >  1:3866910/1‑91 (MQ=255)
                                                            tGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAAtt                                >  1:276444/1‑91 (MQ=255)
                                                             gTCATCCCTCCAATAAAAATAACTCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAAttt                               <  2:3024757/91‑1 (MQ=255)
                                                              tCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTa                              <  1:2939536/91‑1 (MQ=255)
                                                                aTCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAAc                            >  2:678080/1‑91 (MQ=255)
                                                                aTCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAAc                            >  1:1030914/1‑91 (MQ=255)
                                                                    cTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACtttt                        <  1:574484/91‑1 (MQ=255)
                                                                    cTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACtttt                        <  1:2638562/91‑1 (MQ=255)
                                                                    cTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACtttt                        >  1:839781/1‑91 (MQ=255)
                                                                     tCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTc                       >  1:216277/1‑91 (MQ=255)
                                                                      ccAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCt                      >  2:2605908/1‑91 (MQ=255)
                                                                         atttaaaTAGCCCCTATTGGGGTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTCTTCTTaa                   >  2:232773/5‑91 (MQ=37)
                                                                         aTAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTaa                   <  1:1376515/91‑1 (MQ=255)
                                                                           aaaaaTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACt                 <  2:2638324/91‑1 (MQ=255)
                                                                            aaaaTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTg                <  1:2881457/91‑1 (MQ=255)
                                                                            aaaaTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTg                >  1:523182/1‑91 (MQ=255)
                                                                             aaaTAGCCCCTATTGGGGTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGc               >  2:1789461/1‑91 (MQ=38)
                                                                              aaTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATTATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGca              >  2:3273643/1‑91 (MQ=255)
                                                                                tAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGcaca            >  1:2704856/1‑91 (MQ=255)
                                                                                 aaCCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAt           >  1:794402/1‑91 (MQ=255)
                                                                                 aaCCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAt           <  1:2214185/91‑1 (MQ=255)
                                                                                  gccccTATTGGGGTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAtt          >  1:2092180/2‑91 (MQ=255)
                                                                                  gccccTATTGGGGTTATTTAAATGTTATTTTTTTATAATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAtt          >  1:1294516/2‑91 (MQ=255)
                                                                                  aCCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAtt          <  1:2493614/91‑1 (MQ=255)
                                                                                  aCCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAtt          <  2:1699883/91‑1 (MQ=255)
                                                                                  aCCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGAACAtt          >  2:2603058/1‑91 (MQ=255)
                                                                                  aCCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACCTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAtt          <  1:3420454/91‑1 (MQ=255)
                                                                                   ccccTATTGGGGTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTa         >  2:1426488/1‑91 (MQ=255)
                                                                                   ccccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTa         >  1:1413441/1‑91 (MQ=255)
                                                                                    cccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaa        >  2:3337748/1‑91 (MQ=255)
                                                                                    cccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaa        >  2:1593587/1‑91 (MQ=255)
                                                                                     ccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaaa       <  1:362887/91‑1 (MQ=255)
                                                                                      cAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaaaa      <  1:3834377/91‑1 (MQ=255)
                                                                                      cAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaaaa      <  2:2014749/91‑1 (MQ=255)
                                                                                       aaTAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaaaaa     >  2:1045267/1‑91 (MQ=255)
                                                                                       aaTAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaaaaa     <  1:168126/91‑1 (MQ=255)
                                                                                          aGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACaa  >  2:1498881/1‑91 (MQ=255)
                                                                                          |                                                                                          
CGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAA  >  USA300TCH1516_ALE/690011‑690191

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