Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,172,078 T→A 5.3% F43I (TTT→ATT)  leuB → 3‑isopropylmalate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,172,0780TA 5.3% 487.8 / 16.4 152F43I (TTT→ATT) leuB3‑isopropylmalate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (55/88);  new base A (5/3);  total (60/92)
Fisher's exact test for biased strand distribution p-value = 2.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.69e-01

Only 100 of 158 total aligned reads displayed.

TAGGTAATTAATTATGACTTATAACATTGTTGCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCTAAATGGACAGATCCTAACAATCGACCAGAACAAGGA  >  USA300TCH1516_ALE/2171939‑2172218
                                                                                                                                            |                                                                                                                                             
taggtaATTAATTATGACTTATAACATTGTTGCCCTACCTGGTGCTGGAATCGGTCAAGAAATTTTGACCGGATGTCTACCATGGCTCAAAATGATGTGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAAtt                                                                                                                                              >  1:535552/1‑141 (MQ=255)
taggtaATTAATTATGACTTATAACATTGTTGCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAAtt                                                                                                                                              <  4:407591/141‑1 (MQ=255)
taggtaATTAATTATGACTTATAACATTGTTGCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAAtt                                                                                                                                              <  2:202375/141‑1 (MQ=255)
      attaattaTGACTTATAACATTGTTGCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtggtgg                                                                                                                                        <  4:310311/141‑1 (MQ=255)
      attaattaTGACTTATAACATTGTTGCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtggtgg                                                                                                                                        <  1:131371/141‑1 (MQ=255)
      attaattaTGACTTATAACATTGTTGCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtggtgg                                                                                                                                        <  1:968128/141‑1 (MQ=255)
        taattaTGACTTATAACATTGTTGCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTACTTATTAAATAGAGCACTACGAATTTggttgtt                                                                                                                                      >  4:43007/1‑140 (MQ=255)
         aattaTGACTTATAACATTGTTGCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACCAAATTCGTGGTGc                                                                                                                                     <  1:360602/141‑1 (MQ=255)
          attattACTTATAACATTGTTGCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCCCCACGAATTTGGTGGTGcc                                                                                                                                    <  1:470839/141‑1 (MQ=255)
             aTGACTTATAACATTGTTGCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAATGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCctct                                                                                                                                 <  2:229865/141‑1 (MQ=255)
              tGACTTATAACATTGTTGCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTa                                                                                                                                >  1:809787/1‑141 (MQ=255)
               gACTTATAACATTGTTGCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACCAAATTCGTGGTGCTCTAtt                                                                                                                             <  1:238714/141‑1 (MQ=255)
                aCTTATAACATTGTTGCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTAtt                                                                                                                             <  3:356538/141‑1 (MQ=255)
                  ttATAACATTGTTGCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCCCGAATTTGGTGGTGCATCTATTGa                                                                                                                           >  1:516130/1‑141 (MQ=255)
                    ataACATTGTTGCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATa                                                                                                                         >  1:630269/1‑141 (MQ=255)
                     taACATTGTTGCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATAc                                                                                                                        >  2:530993/1‑141 (MQ=255)
                      aaCATTGTTGCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACAACGAATTTGGTGGTGCCTCTATTGATACa                                                                                                                       >  3:65535/1‑141 (MQ=255)
                       aCATTGTTGCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACAt                                                                                                                      <  3:152123/141‑1 (MQ=255)
                       aCATTGTTGCCCTACCTGGTGATAGAATCGGTCAAGAAATTTTGAACGGATTTTTATCATTTTTTTAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATActt                                                                                                                      <  2:352202/141‑3 (MQ=255)
                            gttgCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCCCTATTGATACATTCGGc                                                                                                                 >  1:755829/1‑141 (MQ=255)
                             ttgCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCg                                                                                                                <  1:41657/141‑1 (MQ=255)
                               gCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAg                                                                                                              <  4:53320/141‑1 (MQ=255)
                                 ccTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGCGcc                                                                                                            <  1:769861/141‑1 (MQ=255)
                                 ccTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGcc                                                                                                            <  2:1092721/141‑1 (MQ=255)
                                 ccTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGcc                                                                                                            <  1:530993/141‑1 (MQ=255)
                                  cTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAACAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCt                                                                                                           >  1:466193/1‑141 (MQ=255)
                                   tACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCATCTATTGATACATTCGGCGCGCCtt                                                                                                          >  2:776157/1‑141 (MQ=255)
                                      cTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTaa                                                                                                       >  2:1151013/1‑141 (MQ=255)
                                      cTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTaa                                                                                                       <  2:306768/141‑1 (MQ=255)
                                      cTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTaa                                                                                                       >  1:802693/1‑141 (MQ=255)
                                       tGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAAc                                                                                                      <  3:73958/141‑1 (MQ=255)
                                        ggTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACt                                                                                                     <  2:23578/141‑1 (MQ=255)
                                        ggTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACt                                                                                                     <  1:653443/141‑1 (MQ=255)
                                         gTGATGGAATCGGTCCAGAAATTTTGAACCGTTTTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAATTAGAGCATCATGTATTTTGTGGTTCCTCTATTAATATATTCTTCGAGCCTTTAACtc                                                                                                    >  4:66881/1‑140 (MQ=255)
                                          tGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTga                                                                                                   <  1:1067075/141‑1 (MQ=255)
                                           gATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTgag                                                                                                  >  1:996115/1‑141 (MQ=255)
                                           gATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTgag                                                                                                  >  4:364734/1‑141 (MQ=255)
                                           gATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTgag                                                                                                  >  2:885465/1‑141 (MQ=255)
                                           gATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTCTAACTgag                                                                                                  >  2:311665/1‑141 (MQ=255)
                                            aTGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTgaga                                                                                                 >  4:360598/1‑141 (MQ=255)
                                                aaTCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAAc                                                                                             <  1:701685/141‑1 (MQ=255)
                                                  tCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCt                                                                                           <  2:483361/141‑1 (MQ=255)
                                                  tCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCt                                                                                           <  4:282362/141‑1 (MQ=255)
                                                   cGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACt                                                                                                     <  3:247526/130‑1 (MQ=255)
                                                   cGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACt                                                                                                     >  4:247526/1‑130 (MQ=255)
                                                        cAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATg                                                                                     <  1:846273/141‑1 (MQ=255)
                                                          gAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAGGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCg                                                                                   <  4:440179/141‑1 (MQ=255)
                                                          gAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCg                                                                                   >  2:182561/1‑141 (MQ=255)
                                                          aaaaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCg                                                                                   <  3:104319/140‑1 (MQ=255)
                                                            aaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCgtg                                                                                 <  1:1157812/141‑1 (MQ=255)
                                                             aTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAGACCTTAAATGCgtgt                                                                                <  3:383567/141‑1 (MQ=255)
                                                             aTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCgtgt                                                                                <  2:809787/141‑1 (MQ=255)
                                                                 tGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTaaaa                                                                            <  1:1177276/141‑1 (MQ=255)
                                                                 tGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTaaaa                                                                            >  2:1119729/1‑141 (MQ=255)
                                                                  gAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATAc                                                                                                                        >  2:70275/1‑96 (MQ=255)
                                                                  gAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATAc                                                                                                                        <  1:70275/96‑1 (MQ=255)
                                                                  gAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAag                                                                           >  1:1105081/1‑141 (MQ=255)
                                                                           tctATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGc                                                                  <  1:293633/141‑1 (MQ=255)
                                                                           tctATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGc                                                                  <  1:935056/141‑1 (MQ=255)
                                                                            ctATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAATAGCatctgct                                                                 <  3:334837/141‑7 (MQ=255)
                                                                              aTCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAt                                                               <  2:369509/141‑1 (MQ=255)
                                                                               tCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtt                                                              <  2:516688/141‑1 (MQ=255)
                                                                               tCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtt                                                              <  3:220900/141‑1 (MQ=255)
                                                                                cATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGTTACATTCGGCGAGCCTTTAACTGAGAATACCTTATATGCGTGTAAAAGAGCAGATGCTAttt                                                             <  2:441251/141‑1 (MQ=255)
                                                                                 aTTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtttt                                                            <  3:243999/141‑1 (MQ=255)
                                                                                   tGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTAc                                                          <  2:656364/141‑1 (MQ=255)
                                                                                      ttGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTgg                                                       <  1:485100/141‑1 (MQ=255)
                                                                                      ttGAAATGATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTgg                                                       >  3:399278/1‑141 (MQ=255)
                                                                                         aaaTTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTg                                                    <  1:225207/141‑1 (MQ=255)
                                                                                           aTTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCa                                                  >  1:1065679/1‑141 (MQ=255)
                                                                                           aTTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCa                                                  <  3:434483/141‑1 (MQ=255)
                                                                                            ttATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCaa                                                 <  4:279481/141‑1 (MQ=255)
                                                                                              ataAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATc                                               <  1:182561/141‑1 (MQ=255)
                                                                                               taAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCg                                              >  1:942968/1‑141 (MQ=255)
                                                                                                      taaataTAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGAcc                                       <  1:139774/141‑1 (MQ=255)
                                                                                                       aaataTAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCt                                      >  1:799458/1‑141 (MQ=255)
                                                                                                       aaataTAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCt                                      <  1:522645/141‑1 (MQ=255)
                                                                                                          tataACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCTaaa                                   >  1:398968/1‑141 (MQ=255)
                                                                                                           ataACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCTAAAt                                  >  4:346234/1‑141 (MQ=255)
                                                                                                           ataACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCTAAAt                                  >  1:544029/1‑141 (MQ=255)
                                                                                                             aaCTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCTAAATgg                                <  4:134373/141‑1 (MQ=255)
                                                                                                              aCTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGAc                                        >  1:250106/1‑132 (MQ=255)
                                                                                                              aCTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGAc                                        <  2:250106/132‑1 (MQ=255)
                                                                                                              aCTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCTAAATGGa                               <  2:473661/141‑1 (MQ=255)
                                                                                                               cTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCTAAATGGAc                              <  2:621834/141‑1 (MQ=255)
                                                                                                               cTTTAATTATCAAATAGAGCACCACGAATTTGGTGATGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCTAAATGGAc                              >  4:479803/1‑141 (MQ=255)
                                                                                                                   aaTTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCTAAATGGACAGAt                          >  2:570713/1‑141 (MQ=255)
                                                                                                                   aaTTATCAAATAGAGCACCACGAATTTGATGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCTAAATGGACAGAt                          <  2:867024/141‑1 (MQ=255)
                                                                                                                    aTTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCTAAATGGACAGATc                         <  1:923133/141‑1 (MQ=255)
                                                                                                                      tATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCTAAATGGACAGATCCt                       >  1:336078/1‑141 (MQ=255)
                                                                                                                            aTAGAGGCACCACCAAATTCGTGGTGCTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCTAAATGGACAGAGCCTAACaa                  >  1:1028001/1‑141 (MQ=37)
                                                                                                                              agagCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCTAAATGGACAGATCCTAACAATCg               >  1:367063/1‑141 (MQ=255)
                                                                                                                               gagCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCTAAATGGACAGATCCTAACAATCgc              >  4:209371/1‑140 (MQ=255)
                                                                                                                                aggCACCACCAAATTCGTGGTGCTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCTAAATGGACAGATCCTAACAATCGa              >  2:739663/3‑141 (MQ=38)
                                                                                                                                agCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCTAAATGGACAGATCCTAACAATCGAc             >  1:1022166/1‑141 (MQ=255)
                                                                                                                                  gCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCTAAATGGACAGATCCTAACAATCGAcc            >  1:458759/1‑141 (MQ=255)
                                                                                                                                  gCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCTAAATGGACAGATCCTAACAATCGAcc            >  3:161004/1‑141 (MQ=255)
                                                                                                                                  gCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCTAAATGGACAGATCCTAACAATCGAcc            >  2:982337/1‑141 (MQ=255)
                                                                                                                                    accacGAATTTGGTGGTGCCTCTATTTATACATGCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCTAAATGGACAGATCCTAACAATCGACCAg          <  2:535552/141‑1 (MQ=255)
                                                                                                                                            tttGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCTAAATGGACAGATCCTAACAATCGACCAGAACAAGGa  <  2:1006999/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                             
TAGGTAATTAATTATGACTTATAACATTGTTGCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCTAAATGGACAGATCCTAACAATCGACCAGAACAAGGA  >  USA300TCH1516_ALE/2171939‑2172218

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