Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,141,228 G→A 10.4% E283K (GAA→AAA)  USA300HOU_RS05655 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,141,2280GA10.4% 187.6 / 2.2 74E283K (GAA→AAA) USA300HOU_RS05655hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/35);  new base A (2/7);  total (32/42)
Fisher's exact test for biased strand distribution p-value = 2.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGATAAAGTGCAACAAGACAGTAATAAAAAGATTGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCATAAA  >  USA300TCH1516_ALE/1141088‑1141362
                                                                                                                                            |                                                                                                                                      
aaGATAAAATGCAACAAGACAGAAATAAAAAGATTGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATg                                                                                                                                        <  2:535603/141‑1 (MQ=255)
    tAAAGTGCAACAAGACAGTAATAAAAAGATTGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGTAACATAAAGATGaaaa                                                                                                                                    <  1:467344/141‑1 (MQ=255)
      aaGTGCAACAAGACAGTAATAAAAAGATTGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGaaaaaa                                                                                                                                  >  2:545110/1‑140 (MQ=255)
      aaGTGCAACAAGACAGTAATAAAAAGATTGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAAATGAAAAAt                                                                                                                                  >  1:925285/1‑141 (MQ=255)
          gCAACCAGACAGTAATAAAAAGATTGAAAATGAGCGACCTAAAGGATCACGTACAATAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGaaaaccctcc                                                                                                                              <  2:9268/141‑7 (MQ=255)
          gCAACAAGACAGTAACAAAAAGATTGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAACCATAAAGATGAACAATCaaa                                                                                                                              <  1:589734/141‑1 (MQ=255)
           caacaaGACAGTAATAAAAAGATTGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCaaaa                                                                                                                             >  2:35971/1‑141 (MQ=255)
            aacaaTACAGCAATAAAAAGAATGAAAATAAGCAACCAAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCaaaaa                                                                                                                            <  1:351759/141‑1 (MQ=255)
               aaGACAGTAATAAAAAGATTGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAga                                                                                                                         >  1:472621/1‑140 (MQ=255)
                  aCAGTAATAAAAAGATTGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGgaaa                                                                                                                      >  1:918437/1‑141 (MQ=255)
                  aCAGTAATAAAAAGATTGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGgaaa                                                                                                                      >  2:137983/1‑141 (MQ=255)
                       aataaAAAGATTGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGgaaaagaaa                                                                                                                 <  2:143157/141‑1 (MQ=255)
                       aataaAAAGATAGAAAAGGAGCGAACAAAAGCATCAGGGACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGgaaaagaaa                                                                                                                 <  1:59597/141‑1 (MQ=255)
                        ataaAAAGATTGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGgaaaagaaaa                                                                                                                <  1:411422/141‑1 (MQ=255)
                              aaaaagAAAAAGAGCGACAAAAAAAAACAGGAACAATAAAAATTGAAAATAGCACTACAACAATAAACAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGAAAAAAAAAAAAAAAAGGAAAAGAAAAAagtag                                                                                                          <  2:524390/136‑1 (MQ=255)
                                 aagAAAAAGAGCGACCTAAAACATAAAGAAAAATAAAAGTTGAAAATAGCACACAAAAAAAAAAAAAGGTTGAAAATAATCACAAAGAGCAATAGAAACATAAAGAAGAAAAAACAAAAAAGGAAAACAAAAAagtaaaaa                                                                                                       <  2:417876/139‑5 (MQ=255)
                                   gAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTgaa                                                                                                     >  1:840339/1‑141 (MQ=255)
                                      aaTGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTgaaaaa                                                                                                  >  2:569929/1‑141 (MQ=255)
                                       aTGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTgagaaag                                                                                                 <  1:247583/141‑1 (MQ=255)
                                        agAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTgaaaaaga                                                                                                <  1:933757/140‑1 (MQ=255)
                                              aCCTAAAGCATCAGGTACATTAAAAGTTGAAAACAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGc                                                                                          <  1:537348/141‑1 (MQ=255)
                                               ccTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCg                                                                                         >  1:513899/1‑141 (MQ=255)
                                                aaaaaGAATAAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGAAAAAAAAAAAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGt                                                                                        <  2:524274/139‑1 (MQ=255)
                                                 tAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGtt                                                                                       <  1:569929/141‑1 (MQ=255)
                                                 aaaaGCATCACGTACATCAAAACTGCAAAAGAGCACAACAAAAAAAAAAAAGGTACAAAATAACCACAAAGAGCAACCGAAACATAAACACGAAAAATCAAAAAAGGAAAACAAAAAACTAGTTGAAAAAGAAAAAGCGtt                                                                                       <  2:791726/140‑1 (MQ=255)
                                                  aaaGCATCAGGTACATTAAAAGTTGAAAAAAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTa                                                                                      <  1:59497/141‑1 (MQ=255)
                                                    aGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAACGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTAcc                                                                                    >  1:709761/1‑141 (MQ=255)
                                                    aGCATCAGGTACATTAAAAGGTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTAcc                                                                                    >  2:790551/1‑141 (MQ=255)
                                                        tCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCt                                                                                <  1:591553/141‑1 (MQ=255)
                                                        tCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCt                                                                                <  1:35971/141‑1 (MQ=255)
                                                        acAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCt                                                                                <  1:483761/140‑1 (MQ=255)
                                                           ggTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCtttt                                                                             >  1:512398/1‑141 (MQ=255)
                                                            gTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTa                                                                            >  2:914013/1‑141 (MQ=255)
                                                            gTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTa                                                                            >  2:776713/1‑141 (MQ=255)
                                                            gTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTa                                                                            >  1:148206/1‑141 (MQ=255)
                                                            gTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAAAAAAAAGTAGTTGAAAAAGAAAAAGCGTGAACCGCTGTTa                                                                            >  2:660373/1‑141 (MQ=255)
                                                            gTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAAACACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTCGTAGAAAAAGAAAAAGCGATACAAGCtttaa                                                                            >  2:660391/1‑139 (MQ=255)
                                                                aTTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAATAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATag                                                                        <  1:628997/141‑1 (MQ=255)
                                                                aTTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAATAGCGTTACCAGCTTTTAATag                                                                        <  1:482637/141‑1 (MQ=255)
                                                                   aaaaGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGGAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAgcga                                                                     >  2:219328/1‑141 (MQ=255)
                                                                   aaaaGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGAGGAAAAATCAAAAAAGGAAAAGAAAAAAGGAGTTGAAAAAGAAAAAGCGATACCAGCGTTTAATAgaga                                                                     >  2:251503/1‑141 (MQ=255)
                                                                       gTTGAAAATAGCCCTCAAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAgatgat                                                                 <  2:583546/141‑1 (MQ=255)
                                                                                atcccTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAg                                                        <  1:210396/139‑1 (MQ=255)
                                                                                aGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAAACGAAACATAAAGACAAAAAATCAAAAAAAGAAAAGAAAAAAGTAGTTGAAAAAAAAAAACCGTTACAAGCTTTTAATAGAAATGATGATAGCacg                                                        >  2:373258/1‑139 (MQ=255)
                                                                                aGCCCACCAACAAAAAAAAAGGTTGAAAATAATCAAAAAGAGAAACCAAAACATAAAGATAAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAg                                                        <  1:637191/141‑1 (MQ=255)
                                                                                 gCCCTCCAACAATAAAAAAGGTTGAAAATCATCACAAAGACCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGa                                                       <  1:879702/141‑1 (MQ=255)
                                                                                 gCCCTCCAACAATAAAAAAGGTTGAAAAAAATCACAAAGAGCAACCGAAACAAAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGa                                                       <  2:110454/141‑1 (MQ=255)
                                                                                   ccTCCAACGATAAAAAAGGTAGAAAATAATCACAAAGAGCAGCGGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGGTTTTAATAGTGATGAGGATGGCaataat                                                     >  2:302616/1‑141 (MQ=255)
                                                                                    acccAACAATAAAAAAGGTAAAATAAAAACACAAAGAGCAACCGAAACATAAAAAAAAAAAAAAAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAAta                                                    <  2:731081/139‑1 (MQ=255)
                                                                                        aacaaTAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAAtagtag                                                <  1:791411/141‑1 (MQ=255)
                                                                                              aaaaaaGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAAtt                                          <  2:314477/141‑1 (MQ=255)
                                                                                               aaaaaGGTTGAAAATGATCACAAAGATCAACCGAGACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAGAGCGTTAGCAGCTTTTAATAGAGATGATGATCGCAAGAATAGTAGTCAATTa                                         <  1:574165/141‑1 (MQ=255)
                                                                                                aaaaGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTAt                                        >  1:787496/1‑141 (MQ=255)
                                                                                                   aGGTTGAAAATAAACACCAAGAGAAACCTATCCATAAAGAGGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTa                                     <  1:19319/141‑1 (MQ=255)
                                                                                                    ggTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTGATAGAGATGATGATAGCAAGAATAGTAGTCAATAATCTAg                                    >  2:480973/1‑141 (MQ=255)
                                                                                                    ggTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAg                                    >  1:732231/1‑141 (MQ=255)
                                                                                                    ggTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAg                                    >  2:706318/1‑141 (MQ=255)
                                                                                                             taatCACAAAGAGCAACCAAAACAATATGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTaa                           <  1:775035/141‑1 (MQ=255)
                                                                                                              aatCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTaaa                          <  1:504872/141‑1 (MQ=255)
                                                                                                              aatCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTaaa                          <  1:914013/141‑1 (MQ=255)
                                                                                                                  acaAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTTATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAAc                      <  1:627720/141‑1 (MQ=255)
                                                                                                                     aaGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTtg                   >  2:333360/1‑141 (MQ=255)
                                                                                                                        agCAACCGAAACATAAAGACAAAAAATCAAAAAAGGAAAAGAAAAATGAAGTCGAAAAAGAAAAAAACTGACAAGCGTTTAATAGAGATGATGAAACCAAGAATAATAGTCAATTATCTAGTGATATTAAAGAACTtgata                >  2:271610/1‑140 (MQ=255)
                                                                                                                        agCAACCGAAACATAAAGAAGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTtgatg                <  2:898333/141‑1 (MQ=255)
                                                                                                                             ccGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCa           <  1:790551/141‑1 (MQ=255)
                                                                                                                               aaaaCATAAAAAAAAAAAAAAAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCaaa         <  2:851831/140‑1 (MQ=255)
                                                                                                                                aaaCATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAAt        <  1:790557/141‑1 (MQ=255)
                                                                                                                                aaaCATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAAt        >  2:854258/1‑141 (MQ=255)
                                                                                                                                aaaCATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGAAGTCAATTATCAAGTGATATTAAAGAACTTGATGAACCAAAt        >  2:854234/1‑141 (MQ=255)
                                                                                                                                 aaCATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATAATAAAGAACTTGATGAACCAAATc       >  1:679655/1‑141 (MQ=255)
                                                                                                                                   cATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCAt     >  1:350802/1‑141 (MQ=255)
                                                                                                                                   cATAAAGATGAAAAATAAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCAt     <  2:148206/141‑1 (MQ=255)
                                                                                                                                   cATAAAGATAAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCAt     <  2:900324/141‑1 (MQ=255)
                                                                                                                                   aaaaaaGATAAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTGTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCAt     <  1:928608/138‑1 (MQ=255)
                                                                                                                                      aaaGATGAAAAATCAGAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCATaaa  >  1:701904/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                      
AAGATAAAGTGCAACAAGACAGTAATAAAAAGATTGAAAATGAGCGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCATAAA  >  USA300TCH1516_ALE/1141088‑1141362

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