Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,512,101 C→T 6.7% T96T (ACG→ACA USA300HOU_RS07370 ← Bacterial dynamin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,512,1010CT 6.7% 419.5 / 11.3 136T96T (ACG→ACAUSA300HOU_RS07370Bacterial dynamin‑like protein
Reads supporting (aligned to +/- strand):  ref base C (61/64);  new base T (7/2);  total (69/67)
Fisher's exact test for biased strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 139 total aligned reads displayed.

TTTGAATCAACACCTGGTGTATCTTGCAACGTAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATTAATAAATTTATCAGTGTCGATTTACCTGCAGA  >  USA300TCH1516_ALE/1511961‑1512241
                                                                                                                                            |                                                                                                                                            
  tGAATCAACACCTGGTGTATCTTGCAACGTAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGt                                                                                                                                            <  3:427691/141‑1 (MQ=255)
   gAATCAACACCTGGTGTATCTTGCAACGTAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATATGtt                                                                                                                                           >  1:418571/1‑141 (MQ=255)
    aaTCAACACCTGGTGTATCTTGCAACGTAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGttt                                                                                                                                          >  2:400048/1‑141 (MQ=255)
     aTCAACACCTGGTGTATCTTGCAACGTAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATATGTTTg                                                                                                                                         >  1:711556/1‑141 (MQ=255)
       cAACACCTGGTGTATCTTGCAACGTAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGAt                                                                                                                                       >  2:108124/1‑141 (MQ=255)
         acacCTGGTGTATCTTGCAACGTAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATATGTTTGATTc                                                                                                                                     >  1:719600/1‑141 (MQ=255)
             cTGGTGTATCTTGCAACGTAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCa                                                                                                                                 >  1:594914/1‑141 (MQ=255)
              tGGTGTATCTTGCAACGTAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCg                                                                                                                                    >  2:813993/1‑137 (MQ=255)
              tGGTGTATCTTGCAACGTAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCg                                                                                                                                    <  1:813993/137‑1 (MQ=255)
              tGGTGTATCTTGCAACGTAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCaa                                                                                                                                >  4:232854/1‑141 (MQ=255)
               ggTGTATCTTGCAACGTAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATAAGTTTGATTCGGCaaa                                                                                                                               >  2:32257/1‑141 (MQ=255)
                 tgtATCTTGCAACGTAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAAtt                                                                                                                             >  1:995654/1‑141 (MQ=255)
                  gtATCTTGCAACGTAAACCCATTTTCAAATTTATCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTa                                                                                                                            >  1:975923/1‑141 (MQ=255)
                  gtATCTTGCAACGTAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTa                                                                                                                            >  1:129611/1‑141 (MQ=255)
                  gtATCTTGCAACGTAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTa                                                                                                                            >  2:602898/1‑141 (MQ=255)
                  gtATCTTGCAACGTAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTa                                                                                                                            >  2:396484/1‑141 (MQ=255)
                   tATCTTGCAACGTAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAg                                                                                                                           <  4:158109/141‑1 (MQ=255)
                            aCGTAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATActat                                                                                                                  >  1:599756/1‑141 (MQ=255)
                            aCGTAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAatat                                                                                                                  >  3:191196/1‑141 (MQ=255)
                             cGTAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATc                                                                                                                 <  3:48889/141‑1 (MQ=255)
                              gTAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCg                                                                                                                <  1:640953/141‑1 (MQ=255)
                              gTAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCg                                                                                                                <  3:256166/141‑1 (MQ=255)
                               tAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATATGTTTGATTCGGCAAATTAGCAATAATATCGt                                                                                                               >  1:392641/1‑141 (MQ=255)
                               tAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGt                                                                                                               <  2:980537/141‑1 (MQ=255)
                               tAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGt                                                                                                               <  1:1035533/141‑1 (MQ=255)
                               tAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGt                                                                                                               >  1:1048899/1‑141 (MQ=255)
                                  aCCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGAt                                                                                                            <  2:946175/141‑1 (MQ=255)
                                      aTTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTc                                                                                                        <  4:445918/141‑1 (MQ=255)
                                      aTTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTc                                                                                                        <  1:114699/141‑1 (MQ=255)
                                      aTTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTc                                                                                                        <  4:445849/141‑1 (MQ=255)
                                      aTTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAAGTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTc                                                                                                        <  2:703194/141‑1 (MQ=255)
                                       ttttCAAATTTAGCTAATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTTGTCTTAAGTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATAGCGTGATTGTcg                                                                                                       <  2:1033030/141‑2 (MQ=255)
                                        tttCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATAGTGATTGTCTg                                                                                                      >  2:1078579/1‑141 (MQ=255)
                                           cAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGaaa                                                                                                   >  4:454977/1‑141 (MQ=255)
                                           cAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGaaa                                                                                                   >  2:842642/1‑141 (MQ=255)
                                            aaaTTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAAc                                                                                                  >  2:346496/1‑141 (MQ=255)
                                             aaTTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACt                                                                                                 >  1:696237/1‑141 (MQ=255)
                                              aTTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTg                                                                                                <  1:144474/141‑1 (MQ=255)
                                              aTTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTg                                                                                                <  3:263813/141‑1 (MQ=255)
                                              aTTTAGCTGATTGAAAATTAATTTCTACAGATTAATCCCCGACACTCTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTg                                                                                                <  1:209693/141‑1 (MQ=255)
                                                ttATCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATCCGTTTGATTCGGCGAATTGCCAATAATATCGTGATTGTCTGAAACTGac                                                                                              <  2:864964/141‑1 (MQ=255)
                                                    cTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACaa                                                                                          >  2:524548/1‑141 (MQ=255)
                                                      gATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATa                                                                                        <  1:626970/141‑1 (MQ=255)
                                                       aTTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAg                                                                                       <  4:385245/141‑1 (MQ=255)
                                                            aaaTTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTa                                                                                  <  2:530512/141‑1 (MQ=255)
                                                             aattaattTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTAt                                                                                 >  4:53503/1‑141 (MQ=255)
                                                              attaattTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTAtt                                                                                <  2:741394/141‑1 (MQ=255)
                                                                taattTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTaa                                                                              >  2:1047358/1‑140 (MQ=255)
                                                                taattTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTAc                                                                              <  1:717261/141‑1 (MQ=255)
                                                                 aattTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCGTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACt                                                                             <  1:556934/141‑1 (MQ=255)
                                                                   tttCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTg                                                                           <  1:1026681/141‑1 (MQ=255)
                                                                   tttCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTg                                                                           <  1:960571/141‑1 (MQ=255)
                                                                   tttCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATCACTTg                                                                           >  1:968335/1‑141 (MQ=255)
                                                                    ttCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGt                                                                          >  2:656002/1‑141 (MQ=255)
                                                                           aaTTCAACGTCGACATTTTGGCGATTTATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGt                                                                   <  2:199664/140‑1 (MQ=255)
                                                                            aTTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTa                                                                  >  1:819447/1‑141 (MQ=255)
                                                                               cAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAg                                                               <  2:288889/141‑1 (MQ=255)
                                                                               cAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAg                                                               >  2:1041392/1‑141 (MQ=255)
                                                                                aaCGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAgg                                                              <  2:685539/141‑1 (MQ=255)
                                                                                aaCGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTCGTACAgg                                                              <  3:401934/141‑1 (MQ=255)
                                                                                 aCGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATATGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAggg                                                             >  1:436602/1‑140 (MQ=255)
                                                                                 aCGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATATGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGt                                                             <  2:169348/141‑1 (MQ=255)
                                                                                 aCGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGt                                                             >  1:767320/1‑141 (MQ=255)
                                                                                  cGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATc                                                                                                                 <  3:400082/88‑1 (MQ=255)
                                                                                  cGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATc                                                                                                                 >  4:400082/1‑88 (MQ=255)
                                                                                  cGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTg                                                            >  2:990410/1‑141 (MQ=255)
                                                                                   gTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGa                                                           <  1:133338/141‑1 (MQ=255)
                                                                                            ttGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTaa                                                  <  1:337051/141‑1 (MQ=255)
                                                                                              ggCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGa                                                <  1:32450/141‑1 (MQ=255)
                                                                                                cGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGata                                              >  3:306330/1‑141 (MQ=255)
                                                                                                cGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACAAAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGata                                              >  3:312811/1‑141 (MQ=255)
                                                                                                             ttACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATta                                 >  2:821500/1‑141 (MQ=255)
                                                                                                              tACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATtaa                                <  3:71252/141‑1 (MQ=255)
                                                                                                               aCTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATtaat                               >  4:64894/1‑141 (MQ=255)
                                                                                                                cTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATNTCTTGTTCAATtaata                              <  1:1053871/141‑1 (MQ=255)
                                                                                                                cTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATtaata                              >  1:268052/1‑141 (MQ=255)
                                                                                                                  tcatcaTAATTAGATAATTTGGCATATGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATTAATaaa                            >  2:1054971/1‑141 (MQ=255)
                                                                                                                  tcatcaTAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATTAATaaa                            <  2:333130/141‑1 (MQ=255)
                                                                                                                    atcatAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATTAATAAAtt                          >  1:166274/1‑141 (MQ=255)
                                                                                                                     tcatAATTAGATCATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATTAATAAAttt                         <  3:422066/141‑1 (MQ=255)
                                                                                                                       atAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATTAATAAATTTAt                       <  2:409889/141‑1 (MQ=255)
                                                                                                                          aTTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATTAATAAATTTATCAg                    >  2:24796/1‑141 (MQ=255)
                                                                                                                          aTTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATTAATAAATTTATCAg                    >  2:833186/1‑141 (MQ=255)
                                                                                                                          aTTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATTAATAAATTTATCAg                    <  2:556197/141‑1 (MQ=255)
                                                                                                                           ttAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATTAATAAATTTATCAgt                   <  1:5528/141‑1 (MQ=255)
                                                                                                                             aGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATTAATAAATTTATCAgtgt                 >  2:934000/1‑141 (MQ=255)
                                                                                                                             aGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATTAATAAATTTATCAgtgt                 >  2:1053391/1‑141 (MQ=255)
                                                                                                                              gATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATTAATAAATTTATCAGTGTc                >  1:364959/1‑141 (MQ=255)
                                                                                                                              gATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATTAATAAATTTATCAGTGTc                <  3:126075/141‑1 (MQ=255)
                                                                                                                               aTAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATTAATAAATTTATCAGTGTCg               <  1:357344/141‑1 (MQ=255)
                                                                                                                               aTAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATTAATAAATTTATCAGTGTCg               <  1:104629/141‑1 (MQ=255)
                                                                                                                               aTAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATTAATAAATTTATCAGTGTCg               <  2:317366/141‑1 (MQ=255)
                                                                                                                                      ggCATAGGTTTGATGGGGCAGATTAGAAATAATATCGTGATTGTCTGAAACTGACACAATAGAAGTATTACGTGTTGTTGATAGTGGGGAACTGGGTAAGAGATCTTGTTAAATTAATATTTTTATCATTGTCGATTGAca        <  2:217956/141‑2 (MQ=255)
                                                                                                                                      ggCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATTAATAAATTTATCAGTGTCGATTTAcc        >  1:1039827/1‑141 (MQ=255)
                                                                                                                                       gCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATTAATAAATTTATCAGTGTCGATTTACCt       >  2:323034/1‑141 (MQ=255)
                                                                                                                                       gCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATTAATAAATTTATCAGTGTCGATTCACCt       >  1:777303/1‑141 (MQ=255)
                                                                                                                                        cATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATTAATAAATTTATCAGTGTCGATTTACCTg      <  1:951805/141‑1 (MQ=255)
                                                                                                                                          tACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATTAATAAATTTATCAGTGTCGATTTACCTGCa    <  1:524548/141‑1 (MQ=255)
                                                                                                                                           aCGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATTAATAAATTTATCAGTGTCGATTTACCTGCAg   <  1:990410/141‑1 (MQ=255)
                                                                                                                                            cGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTAGTTGGTACAGCTGCTCTTCGTAAGATATCAACTCCAATTAATAAATTTATCAGTGTCGATTTACCTGCAGa  >  1:158026/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                            
TTTGAATCAACACCTGGTGTATCTTGCAACGTAAACCCATTTTCAAATTTAGCTGATTGAAAATTAATTTCTACAGATTCAACGTCGACATTTTGGCGATTCATTTCCCTTACTTCATCATAATTAGATAATTTGGCATACGTTTGATTCGGCAAATTAGCAATAATATCGTGATTGTCTGAAACTGACACAATAGCAGTATTACTTGTCGTTGGTACAGGTGAACTTGGTAAGATATCTTGTTCAATTAATAAATTTATCAGTGTCGATTTACCTGCAGA  >  USA300TCH1516_ALE/1511961‑1512241

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: