Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 629,681 C→A L84I (CTA→ATA)  USA300HOU_RS02990 → 3‑hexulose‑6‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE629,6810CA98.7% 227.4 / NA 74L84I (CTA→ATA) USA300HOU_RS029903‑hexulose‑6‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (33/40);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGAAGGTTTACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGG  >  USA300TCH1516_ALE/629543‑629810
                                                                                                                                          |                                                                                                                                 
aaCCAATGTTTACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGTTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTataataata                                                                                                                                 <  1:187623/141‑1 (MQ=255)
    aaGGTTTACAAGCAAGTAAACATATGGCAGACAACATTAGTAAGGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGgtg                                                                                                                             <  2:709963/141‑1 (MQ=255)
      ggTTTACCAGCAGTTAAACATATGGCAGGCAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGATGGTTATGAAGTTAGCCAAGCAATTAAATTTGCCGCGGATGTAATGACAATAATAGGTGtt                                                                                                                           >  2:3684/1‑141 (MQ=255)
      ggTTTACCAGCAGTTAAACATATGGCAGACCACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATAATGAAGTTAGCCAAGACATTAAATTTGGCGCGGAGGAAATCACAATAATAgatctt                                                                                                                           >  2:176337/1‑136 (MQ=255)
      ggTTTACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGtt                                                                                                                           >  2:713738/1‑141 (MQ=255)
      ggTTTACCAGCAGTGAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGtt                                                                                                                           >  1:623834/1‑141 (MQ=255)
       gTTTACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTg                                                                                                                          >  1:1130152/1‑141 (MQ=255)
       gTTTACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTg                                                                                                                          >  1:1133844/1‑141 (MQ=255)
       gTTTACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTg                                                                                                                          >  2:610167/1‑141 (MQ=255)
        tttACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGAGGATGGAATTACAATAATAGGTGTTGc                                                                                                                         >  1:214674/1‑141 (MQ=255)
         ttACCAGCAGTTAACCATATGGCAGACAGCATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCa                                                                                                                        >  1:779868/1‑141 (MQ=255)
          tACCAGCAGTTAAACATATGGCAGACAACATAAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCag                                                                                                                       <  1:496920/141‑1 (MQ=255)
          tACCAGAAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCag                                                                                                                       <  1:32607/141‑1 (MQ=255)
             cagcagTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCagaag                                                                                                                    <  1:1044882/141‑1 (MQ=255)
                   taaaaCATATGGCAGAAAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATATAATTACAATAATAGGTGTTGCAGAAGATGCAt                                                                                                              <  1:366853/139‑1 (MQ=255)
                       aCATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAAt                                                                                                          <  2:677225/141‑1 (MQ=255)
                        cgtatGGCAGACAACATTAGTAATGTAAAAGTATTAGTAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTCGCCAAGCCATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAAtt                                                                                                         <  2:697921/139‑1 (MQ=255)
                        cATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAAtt                                                                                                         <  1:1149696/141‑1 (MQ=255)
                        aatagGGCAGAAAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAAtt                                                                                                         <  1:666119/140‑1 (MQ=255)
                            tGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAag                                                                                                     <  1:473950/141‑1 (MQ=255)
                                  acaacaTTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTa                                                                                               >  2:406531/1‑141 (MQ=255)
                                    aacaTTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCCATTAAAGCAGCTAtt                                                                                             >  2:186048/1‑141 (MQ=255)
                                    aacaTTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTAtt                                                                                             >  1:893442/1‑141 (MQ=255)
                                     acaTTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTg                                                                                            >  2:170683/1‑141 (MQ=255)
                                     acaTTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTg                                                                                            >  1:364639/1‑141 (MQ=255)
                                     acaTTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGCTAGCCAAGCAAGTAAATTTGGCGCGGATGTAATTACAACAAGCGGTGTGGCAGAGGATGCATCAGTTAAAGAAGCtttgg                                                                                            >  2:330320/1‑137 (MQ=255)
                                       aTTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTgaa                                                                                          >  2:230655/1‑141 (MQ=255)
                                        ttAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGTTGCAGCTGATTATGAAGTTAGCTAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTgaag                                                                                         <  2:927236/141‑1 (MQ=255)
                                         tAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAATCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTgaaga                                                                                        <  2:864834/141‑1 (MQ=255)
                                         tAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGGAGAAGATGCATCAATTAAAGCAGCTATTgaaga                                                                                        <  2:430497/141‑1 (MQ=255)
                                         tAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTgaaga                                                                                        >  1:1148480/1‑141 (MQ=255)
                                         tAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTgaaga                                                                                        <  2:140215/141‑1 (MQ=255)
                                          aGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCGTCAATTAAAGCAGCTATTgaagaa                                                                                       >  2:156191/1‑141 (MQ=255)
                                           gTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGaagaag                                                                                      <  2:393468/141‑1 (MQ=255)
                                                gTAAAAGTATTAGCAGACATGAAGATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCAt                                                                                 <  2:904970/141‑1 (MQ=255)
                                                gTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCAt                                                                                 <  2:469078/141‑1 (MQ=255)
                                                   aaaGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATaaa                                                                              <  1:704804/141‑1 (MQ=255)
                                                       tATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGCAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAaata                                                                          >  2:1025065/1‑141 (MQ=255)
                                                         ttAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAAGTTGGCGCGGATGTCATTACAATAATAGGTGTTCCAGAAGATGCATCAATTAAATAGTCTATCGGAGAAGCTCATAAAaataat                                                                        >  1:396952/1‑141 (MQ=255)
                                                          tAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAataata                                                                       >  1:532012/1‑141 (MQ=255)
                                                          tAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAataata                                                                       <  1:901960/141‑1 (MQ=255)
                                                          tAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAataata                                                                       >  1:921744/1‑141 (MQ=255)
                                                          tAGCAGACATGAAAATTATAGATGCAGCGGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAataata                                                                       <  2:1039782/141‑1 (MQ=255)
                                                            gCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATaaa                                                                     >  1:99402/1‑141 (MQ=255)
                                                              aGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACa                                                                   <  1:406531/141‑1 (MQ=255)
                                                                  aTGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCACCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTa                                                               <  2:779868/141‑1 (MQ=255)
                                                                       aaTTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGt                                                          <  2:272643/141‑1 (MQ=255)
                                                                          tataGTTGAAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGa                                                       <  2:315487/141‑1 (MQ=255)
                                                                              gATGCAGCTGATTATGAAGTTAGCCAAGAAATTAAATTTGGCGCGTATTGAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATg                                                   <  1:514864/141‑1 (MQ=255)
                                                                                    gCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCt                                             >  1:778706/1‑141 (MQ=255)
                                                                                       gATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGtt                                          >  2:1189188/1‑141 (MQ=255)
                                                                                              aatttAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATCACTAGTTGATATGATTGCTGTTCAAGAtt                                   <  1:355349/138‑1 (MQ=255)
                                                                                              aaGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAATAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGAtt                                   <  1:1056694/141‑1 (MQ=255)
                                                                                              aaGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGAtt                                   <  1:894554/141‑1 (MQ=255)
                                                                                                gTTAGCCAAGCAATTAAAATTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTa                                 <  1:1162653/141‑1 (MQ=255)
                                                                                                 ttAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAg                                <  1:966782/141‑1 (MQ=255)
                                                                                                 ttAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAg                                <  1:32507/141‑1 (MQ=255)
                                                                                                  tAGCCAAGCAATTAAATTTTGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGTTATGATTGCTGTTCAAGATTTAGa                               <  1:762787/141‑1 (MQ=255)
                                                                                                  tAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGa                               <  2:379709/141‑1 (MQ=255)
                                                                                                  tAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGa                               <  2:385329/141‑1 (MQ=255)
                                                                                                   aGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGGTCATAAAAATAATAAACAATTAATAGATGATATGATTGCTGTTCAAGATTTAGaa                              >  1:834333/1‑141 (MQ=255)
                                                                                                      cAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGaaaaa                           >  2:387785/1‑141 (MQ=255)
                                                                                                        aGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACg                         <  1:417284/141‑1 (MQ=255)
                                                                                                        aGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACg                         >  2:511461/1‑141 (MQ=255)
                                                                                                             ttAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCaa                    >  2:687944/1‑141 (MQ=255)
                                                                                                              tAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCaaa                   <  1:950436/141‑1 (MQ=255)
                                                                                                                aaTTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAg                 >  1:712210/1‑141 (MQ=255)
                                                                                                                aaTTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAAAAATTAATAGTTGATGTGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAg                 >  1:1112136/1‑141 (MQ=255)
                                                                                                                          ggATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGa       >  2:165067/1‑141 (MQ=255)
                                                                                                                           gATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGaa      >  1:792695/1‑141 (MQ=255)
                                                                                                                            aTGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGaaa     <  1:1083901/141‑1 (MQ=255)
                                                                                                                            aTGTAATTACAATAATAGGTCTTGCAGAAGATGCATCAATCAAAGCCGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGaaa     <  1:102249/141‑1 (MQ=255)
                                                                                                                              gTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATg   <  2:1130152/141‑1 (MQ=255)
                                                                                                                               ttaTTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATgg  <  2:571741/139‑1 (MQ=255)
                                                                                                                               tAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATAGGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATgg  >  2:616115/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                 
AACGAAGGTTTACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGG  >  USA300TCH1516_ALE/629543‑629810

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