Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,172,075 G→C 5.1% E42Q (GAA→CAA)  leuB → 3‑isopropylmalate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,172,0750GC 5.1% 526.2 / 18.9 156E42Q (GAA→CAA) leuB3‑isopropylmalate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (56/91);  new base C (5/3);  total (61/95)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.26e-01

Only 100 of 158 total aligned reads displayed.

AAGGTAGGTAATTAATTATGACTTATAACATTGTTGCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCTAAATGGACAGATCCTAACAATCGACCAGAACA  >  USA300TCH1516_ALE/2171935‑2172214
                                                                                                                                             |                                                                                                                                            
aaGGTAGGTAATTAATTATGACTTATAACATTGTTGCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACg                                                                                                                                              <  4:179380/141‑1 (MQ=255)
aaGGTAGGTAATTAATTATGACTTATAACATTGTTGCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACg                                                                                                                                              <  2:122538/141‑1 (MQ=255)
    taggtaATTAATTATGACTTATAACATTGTTGCCCTACCTGGTGCTGGAATCGGTCAAGAAATTTTGACCGGATGTCTACCATGGCTCAAAATGATGTGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAAtt                                                                                                                                          >  1:535552/1‑141 (MQ=255)
    taggtaATTAATTATGACTTATAACATTGTTGCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAAtt                                                                                                                                          <  4:407591/141‑1 (MQ=255)
       gtaATTAATTATGACTTATAACATTGTTGCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtgg                                                                                                                                       >  2:209910/1‑141 (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)
                  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)
                           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)
                                  tgCCCTACCTGGTGCTGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGCAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGCTGACTCTATTGATACATTCCGCGa                                                                                                           >  2:109707/1‑141 (MQ=255)
                                   gCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAg                                                                                                          <  4:53320/141‑1 (MQ=255)
                                     ccTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGCGcc                                                                                                        <  1:769861/141‑1 (MQ=255)
                                     ccTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGcc                                                                                                        <  1:530993/141‑1 (MQ=255)
                                      cTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCt                                                                                                       >  2:1172061/1‑141 (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                                                                                              >  2:885465/1‑141 (MQ=255)
                                               gATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTgag                                                                                              >  4:364734/1‑141 (MQ=255)
                                               gATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTgag                                                                                              >  1:996115/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                                                                                       <  4:282362/141‑1 (MQ=255)
                                                      tCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCt                                                                                       <  2:483361/141‑1 (MQ=255)
                                                       cGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACt                                                                                                 >  4:247526/1‑130 (MQ=255)
                                                       cGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACt                                                                                                 <  3:247526/130‑1 (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)
                                                                 aTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCgtgt                                                                            <  2:809787/141‑1 (MQ=255)
                                                                     tGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTaaaa                                                                        <  3:84796/141‑1 (MQ=255)
                                                                     tGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTaaaa                                                                        <  1:1177276/141‑1 (MQ=255)
                                                                      gAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATAc                                                                                                                    <  1:70275/96‑1 (MQ=255)
                                                                      gAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATAc                                                                                                                    >  2:70275/1‑96 (MQ=255)
                                                                      gAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAag                                                                       >  1:1105081/1‑141 (MQ=255)
                                                                          gcaTCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAATAGCa                                                                   <  2:1022549/139‑1 (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)
                                                                                           tGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTggg                                                  <  3:52972/141‑1 (MQ=255)
                                                                                             aaaTTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTg                                                <  1:225207/141‑1 (MQ=255)
                                                                                               aTTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCa                                              >  1:1065679/1‑141 (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                              >  1:544029/1‑141 (MQ=255)
                                                                                                               ataACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCTAAAt                              >  4:346234/1‑141 (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)
                                                                                                                             cAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGTTATTTTACTGGGTGCAATCGGTGGACCTAAATGGACAGATCCTAAc                <  3:416695/141‑1 (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)
                                                                                                                                             |                                                                                                                                            
AAGGTAGGTAATTAATTATGACTTATAACATTGTTGCCCTACCTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTGCAATCGGTGGACCTAAATGGACAGATCCTAACAATCGACCAGAACA  >  USA300TCH1516_ALE/2171935‑2172214

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: