Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,781,307 T→C 5.1% K471E (AAG→GAG)  lysP_2 ← Lysine‑specific permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,781,3070TC 5.1% 218.3 / 4.1 79K471E (AAG→GAG) lysP_2Lysine‑specific permease
Reads supporting (aligned to +/- strand):  ref base T (39/34);  new base C (2/2);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.92e-02

TTAATCATATCAACAATTCACATTTAAGTCAACACTATAACTGTAATTATTTTTATTTTAACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACATCTTGACCAATAATTACTACAATACATAAAAAGCCAGCAA  >  USA300TCH1516_ALE/1781167‑1781441
                                                                                                                                            |                                                                                                                                      
ttaATCATATCAACAATTCACATTTAAGTCAACACTATAACTGTAATTATTTTTATTTTAACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCtt                                                                                                                                        >  2:158037/1‑141 (MQ=255)
  aaTCATATCAACAATTCACATTTAAGTCAACACTATAACTGTAATTATTTTTATTTTAACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTAt                                                                                                                                      >  2:659437/1‑141 (MQ=255)
   aTCATATCAACAATTCACATTTAAGTCAACACTATAACTGTAATTATTTTTATTTTAACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATg                                                                                                                                     <  2:448672/141‑1 (MQ=255)
   aTCATATCAACAATTCACATTTAAGTCAACACTATAACTGTAATTATTTTTATTTTAACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATg                                                                                                                                     <  1:805662/141‑1 (MQ=255)
   aTCATATCAACAATTCACATTTAAGTCAACACTATAACTGTAATTATTTTTATTTGAACCTTTTATTTCATCCATTGACACGTCTTGACTTAAATCTACTCGTTCTATTGGTATTTTTTTCGTTTTCTATCGAGGCTTATg                                                                                                                                     >  1:577089/1‑141 (MQ=255)
     cATATCAACAATTCACATTTAAGTCAACACTATAACTGTAATTATTTTTATTTTAACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGat                                                                                                                                   <  2:135192/141‑1 (MQ=255)
     cATATCAACAATTCACATTTAAGTCAACACTATAACTGTAATTATTTTTATTTTAACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGaa                                                                                                                                   >  2:655898/1‑140 (MQ=255)
             caaTTCACATTTAAGTCAACACTATAACTGTAATTATTTTTATTTTAACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAg                                                                                                                           >  1:795011/1‑141 (MQ=255)
               aTTCACATTTAAGTCACCACTATAACTGTAATTATTTTTATTTTAACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGaa                                                                                                                         <  2:500092/141‑1 (MQ=255)
               aTTCACATTTAAGTCAACACTATAACTGTAATTATTTTTATTTTAACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGaa                                                                                                                         >  2:528723/1‑141 (MQ=255)
                  cacaTTTAAGTCAACACTATAACTGTAATTATTTTTATTTTAACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGc                                                                                                                      <  2:472823/141‑1 (MQ=255)
                   acaTTTAAGTCAACACTATAACTGTAATTATTTTTATTTTAACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTCATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCt                                                                                                                     <  2:694889/141‑1 (MQ=255)
                   acaTTTAAGTCAACACTATAACTGTAATTATTTTTATTTTAACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGACTTTTTTTCGTTTTATATCGAAGCTGATGACAAATAAAGAATGCt                                                                                                                     <  1:515845/141‑1 (MQ=255)
                    caTTTAAGTCAACACTATAACTGTAATTATTTTTATTTTAACCTTTTGTTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAGGAATGCTa                                                                                                                    >  2:326184/1‑141 (MQ=255)
                    caTTTAAGTCAACACTATAACTGTAATTATTTTTATTTTAACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTa                                                                                                                    >  2:712630/1‑141 (MQ=255)
                    caTTTAAGTCAACACTATAACTGTAATTATTTTTATTTTAACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTa                                                                                                                    >  2:416076/1‑141 (MQ=255)
                            tccacacTATAACTGTAATTATTTTTACTTTAACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACt                                                                                                            <  1:191321/138‑1 (MQ=255)
                            tCAATATTATAATTGTAATTATTTTTTTTTTAACCTTTTATTTCATCCATTGACACGTCTTGACGTAAAACTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACt                                                                                                            <  2:332973/141‑1 (MQ=255)
                                       aCTGTAATTATTTTTATTTTAACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATTTCGACGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCat                                                                                                 >  2:37128/1‑141 (MQ=255)
                                         tGTAATTATTTTTATTTTAACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGATATTCCCatat                                                                                               >  2:598634/1‑141 (MQ=255)
                                          gTAATTATTTTTATTTTAACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCAtata                                                                                              >  2:305333/1‑141 (MQ=255)
                                              ttatttttatttTAACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTa                                                                                          >  1:397056/1‑141 (MQ=255)
                                              ttatttttatttTAACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTCATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTa                                                                                          >  1:185882/1‑141 (MQ=255)
                                              ttatttttatttTAACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTCATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTa                                                                                          >  1:186007/1‑141 (MQ=255)
                                               tatttttattttAACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTCATTCCGCTTATTTAATTTTTTTCGTTTTATATCGAAGCTTATGCTAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTaa                                                                                         <  1:926993/141‑1 (MQ=255)
                                                  ttttattttaACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATAc                                                                                                             <  2:194464/118‑1 (MQ=255)
                                                  ttttattttaACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATAc                                                                                                             >  1:194464/1‑118 (MQ=255)
                                                       ttttaGCCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCCACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTaagaag                                                                                 <  2:351997/141‑1 (MQ=255)
                                                       ttttaACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATAAACGTACTTaaaaaa                                                                                 >  2:534202/1‑140 (MQ=255)
                                                                tttATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACtaaaa                                                                        <  2:467906/141‑1 (MQ=255)
                                                                     ttCATCCATTGACACGTCTTGACGTAAATCTACTTATTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACtaaaattaaa                                                                   <  2:46296/141‑1 (MQ=255)
                                                                      tCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTTGTTTAAGAGGGGAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACtaaaattaaaa                                                                  <  1:927278/141‑1 (MQ=255)
                                                                       cATCCATTGACACGTCTTGACATAAATCTACTTGTTCTAATAGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACtaaaattaaaat                                                                 <  1:787181/141‑1 (MQ=255)
                                                                        aTCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTTTTTTTTTTTCTAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTATGACGCGCCtaaaatttaaatt                                                                >  2:588816/1‑140 (MQ=255)
                                                                        aTCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAAtc                                                                >  2:311462/1‑141 (MQ=255)
                                                                             ttGACACGCCATGACAAAAATCTACTTGTTCTAGTGGAATTGTTTTCGTTCTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCACCTg                                                           >  2:692608/1‑141 (MQ=255)
                                                                                acacGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGttt                                                        >  2:424511/1‑141 (MQ=255)
                                                                                acacGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAACTTATGATAAATAAAGAATGCTAAAAATAATGGAATTACCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGttt                                                        >  2:183890/1‑141 (MQ=255)
                                                                                       tttacctAAATCTATTTGTTATAATTGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATCCCTACTTAAGATGCGACTAACATTAAAATCTCCTGTTTTAATaaa                                                 <  1:597544/135‑1 (MQ=255)
                                                                                       ttGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATaaa                                                 <  2:900960/141‑1 (MQ=255)
                                                                                       ttGACGTAAATCTACTTCTTCTAAATGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATaaa                                                 <  1:401745/141‑1 (MQ=255)
                                                                                           cGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCa                                             >  1:780862/1‑141 (MQ=255)
                                                                                            gTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCaa                                            <  2:397056/141‑1 (MQ=255)
                                                                                            gTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTCATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCaa                                            <  2:185882/141‑1 (MQ=255)
                                                                                            gTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTCATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCaa                                            <  2:186007/141‑1 (MQ=255)
                                                                                             tAAATCTTCTTGTTTTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAAc                                           <  2:68101/141‑1 (MQ=255)
                                                                                             tAAATCTACTTTTTCTAATGGAATTTTTTTCGTTTTAAATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAAc                                           <  1:222669/141‑1 (MQ=255)
                                                                                             tAAATCGTCTTGTTCTAATGGAATTCTTTTCGTTTTATATCGAAGCTGATGATAAATAAAGAATGCTAAAAATACTGGAATGCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAAc                                           <  2:68021/141‑1 (MQ=255)
                                                                                             aaaaCCTACTTGCTCTCATGGAATTTTCTTCGTTTTACAACGAAGCTTATGACAAATAAAGAATGCTAAAAATACTGGAATTCCCATAGACGTAATCAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAAc                                           <  1:399498/140‑1 (MQ=255)
                                                                                              aaaTCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACa                                          <  1:161897/141‑1 (MQ=255)
                                                                                                aTCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACATc                                        >  1:862124/1‑141 (MQ=255)
                                                                                                  cTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACATCtt                                      >  1:114099/1‑141 (MQ=255)
                                                                                                      ttGTTCTAATGGATTTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACATCTTGAcc                                  <  2:23830/141‑1 (MQ=255)
                                                                                                      ttGTTCTAAAGGAATTTTTTTCGTTATATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACATCTTGAcc                                  <  2:941109/141‑1 (MQ=255)
                                                                                                       tGTTCTAATGGAATTTTTTTCGCCATCGTTGGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTCAAATTAAAATCTCGTGTTTTAATAAAGTCAACATCTTGACCa                                 >  2:668428/1‑141 (MQ=255)
                                                                                                            taaatgAATTTTTTTCGTTTTATATCCAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGTAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACATCTTGACCaataat                            <  1:223830/136‑1 (MQ=255)
                                                                                                              aTGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACATCTTGACCAATAATta                          >  1:96945/1‑141 (MQ=255)
                                                                                                                 gAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAAAATCTTGACAAATAATtacta                       >  2:103744/1‑141 (MQ=255)
                                                                                                                  aaTTTTTTTCGTTTTATTTCGAAGCTTATGATAAATAAAAAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACATCTTGACCAATAATtactac                      <  2:345570/141‑1 (MQ=255)
                                                                                                                  aaTTGTTTTCGTGTTATATCGAAGCTTATGATAAATAAAGAATGCTTAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAAGCTCCTCTTTTAATAAAGTCAACATCTTGACCAATAATtactac                      >  1:664842/1‑141 (MQ=255)
                                                                                                                    tttttttCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACATCTTGACCAATAATTACTACaa                    >  2:70241/1‑141 (MQ=255)
                                                                                                                        tttCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACATCTTGACCAATAATTACTACAATACa                <  1:659437/141‑1 (MQ=255)
                                                                                                                         ttCGTTTTATATCGAAGCTTATGATAAATAGGGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACATCTTGACAAATAATTACTACAATACAt               >  1:224890/1‑141 (MQ=255)
                                                                                                                         ttCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACATCTTGACCAATAATTCCTACATCACAt               >  2:944217/1‑141 (MQ=255)
                                                                                                                         ttCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACATCTTGACCAATAATTACTACAATACAt               >  1:313198/1‑141 (MQ=255)
                                                                                                                         ttCGTTTTATATCGAAGCTTATGATAAATAAAAAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACATCTTGACCAATAATTACTACAATACAt               >  1:597689/1‑141 (MQ=255)
                                                                                                                         ttCGTCTTATATCAAAGCTTATAATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAAGTAAAATCTCCTGTTTTAATAAAGTCAACATCTTGACCAATAATAACTACAATACag               <  2:296029/141‑2 (MQ=255)
                                                                                                                          tCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACATCTTGACCAATAATTACTACAATACATa              <  1:796335/141‑1 (MQ=255)
                                                                                                                           cGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACATCTTGACCAATAATTACTACAATACATaa             >  2:755624/1‑141 (MQ=255)
                                                                                                                             ttttATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACATCTTGACCAATAATTACTACAATACATaaaa           >  2:851155/1‑141 (MQ=255)
                                                                                                                                 atatCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAGTTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACATCTTGACCACTAATTCCAACAATAAATTAAAAgaa       >  2:492800/1‑139 (MQ=255)
                                                                                                                                 atatCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTAGTAAGAAAGTCAACATCTTGACCAATAATTACTACAATACATAAAAAGcc       >  1:306888/1‑141 (MQ=255)
                                                                                                                                 atatCAAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATGCGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACATCTTGACCAATAATTACTACAATACCAAAAAAGcc       >  2:217597/1‑141 (MQ=255)
                                                                                                                                   atggAAGTTTATGATAACCACAGAATGATAAACATACTGGAATCCGCATATACGTCATTAAGAAGTGACTAAAATTACAATCTCCTGTTTTAATAAAGTCAACATCTTGACCAATAATTACTACAATACATAAAAAGCCAg     <  1:902826/138‑1 (MQ=255)
                                                                                                                                   atCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACATCTTGACCAATAATTACTACAATACATAAAAAGCCAg     <  2:19083/141‑1 (MQ=255)
                                                                                                                                   atCGAAGCTTAAGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACATCTTGACCAATAATTACTACAATACATAAAAAGCCAg     <  2:800018/141‑1 (MQ=255)
                                                                                                                                   atCAAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACATCTTGACCAATAATTACTACAATACATAAAAAGCCAg     <  2:503612/141‑1 (MQ=255)
                                                                                                                                     cGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACATCTTGACCAATAATTACTACAATACATAAAAAGCCAGCa   >  2:787664/1‑141 (MQ=255)
                                                                                                                                     cGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACATCTTGACCAATAATTACTACAATACATAAAAAGCCAGCa   >  2:872014/1‑141 (MQ=255)
                                                                                                                                      gAAGCTTATGATAAATAAAGAATGGAATCCAAACTAGAATTCCCATATACGTAATTAAGAAGAGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACATCTTGACCAATAATTACTACAATACATACAAAGCCAGCaa  >  1:693003/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                      
TTAATCATATCAACAATTCACATTTAAGTCAACACTATAACTGTAATTATTTTTATTTTAACCTTTTATTTCATCCATTGACACGTCTTGACGTAAATCTACTTGTTCTAATGGAATTTTTTTCGTTTTATATCGAAGCTTATGATAAATAAAGAATGCTAAAAATACTGGAATTCCCATATACGTAATTAAGAAGCGACTAAAATTAAAATCTCCTGTTTTAATAAAGTCAACATCTTGACCAATAATTACTACAATACATAAAAAGCCAGCAA  >  USA300TCH1516_ALE/1781167‑1781441

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: ‑