Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 110,651 C→T 8.1% D170D (GAC→GAT USA300HOU_RS00515 → putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE110,6510CT 8.1% 181.6 / 2.2 69D170D (GAC→GATUSA300HOU_RS00515putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (23/39);  new base T (3/4);  total (26/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AAGATATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAG  >  USA300TCH1516_ALE/110513‑110787
                                                                                                                                          |                                                                                                                                        
aaGATATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACg                                                                                                                                         >  1:10353/1‑140 (MQ=255)
 aGATAGACGGATAAAATGGTCGATAATAAAATAATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGATCAATATGGCGACTTTAAGAATTTGAAAAATTATAAAGACgg                                                                                                                                        <  2:404241/140‑1 (MQ=18)
    tatCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGgaga                                                                                                                                     <  2:376354/140‑1 (MQ=255)
    tatCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGgaga                                                                                                                                     >  1:294358/1‑140 (MQ=255)
         ggTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATAttt                                                                                                                                <  1:628334/140‑1 (MQ=255)
                agagacaataataAAATCATTCCAACAAAAGAAAATAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATaa                                                                                                                         <  2:286836/133‑1 (MQ=255)
                  ggTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAACAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAAta                                                                                                                       >  2:733756/1‑139 (MQ=255)
                  ggTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATc                                                                                                                       >  1:262295/1‑140 (MQ=255)
                   gTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAACATTTGAAAAATTATAAAGACGGAGATCTTTCATATAATcc                                                                                                                      >  1:122501/1‑140 (MQ=255)
                   gTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATcc                                                                                                                      <  1:344640/140‑1 (MQ=255)
                     cGAAAATAAAAACATGCCAACAAAAGAACTTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATTATCCag                                                                                                                    <  2:642302/140‑1 (MQ=18)
                      gATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCaga                                                                                                                   >  1:247833/1‑140 (MQ=255)
                             aaaTCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCg                                                                                                            <  1:700373/140‑1 (MQ=255)
                              aaTCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTGTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCgc                                                                                                           >  2:36248/1‑139 (MQ=255)
                                tcaatcaaACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGt                                                                                                         <  2:58519/140‑1 (MQ=255)
                                tCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTAATATAATCAAGAGATGCAGAGt                                                                                                         >  2:64769/1‑140 (MQ=255)
                                   ttCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTAt                                                                                                      >  2:199041/1‑140 (MQ=255)
                                   ttCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCAAGAGGTGCCGAGTTAt                                                                                                      >  1:76999/1‑140 (MQ=255)
                                    tCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTAGGTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGAGCCGAGTTAtt                                                                                                     >  2:373844/1‑140 (MQ=255)
                                       aacaaAAGAAAGTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCgg                                                                                                  <  1:676796/140‑1 (MQ=255)
                                          aaaaGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTa                                                                                               <  1:383598/140‑1 (MQ=255)
                                           aaaGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTaa                                                                                              <  1:94817/140‑1 (MQ=255)
                                            aaGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTaaa                                                                                             <  1:229652/140‑1 (MQ=255)
                                              gAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAata                                                                                           <  2:66028/140‑1 (MQ=255)
                                              aaaaaaaaaGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAata                                                                                           <  2:606842/134‑1 (MQ=255)
                                               aaaTTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTCTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAatat                                                                                          <  1:449452/140‑1 (MQ=255)
                                               aaaCTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAAGTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAatat                                                                                          <  1:303802/140‑1 (MQ=255)
                                                aaTTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCAAGAGGTGCCGAGTTATTCGGCTAAACATc                                                                                         >  2:188873/1‑140 (MQ=255)
                                                  ttAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCaa                                                                                       <  2:481468/140‑1 (MQ=255)
                                                   tAAAAATGAAAAAAAAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAt                                                                                      <  2:262838/140‑1 (MQ=255)
                                                         tGAAAAAATAAAAAAAGAAACTGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTa                                                                                <  1:736148/140‑1 (MQ=255)
                                                             aaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatga                                                                            <  1:272522/140‑1 (MQ=255)
                                                                aTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAAtgatga                                                                         <  1:644421/140‑1 (MQ=255)
                                                                  aaaaaacGAAATGGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGAtt                                                                       <  2:263101/140‑1 (MQ=255)
                                                                  aaaaaaaGAAATTGAAAACTTTAAGTTCTTTGTTCAAGATGGCGACTTTAAAAATTTGAAAAATTATAAAGATGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGAtt                                                                       <  2:102515/140‑1 (MQ=255)
                                                                   aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATta                                                                      >  1:644356/1‑140 (MQ=255)
                                                                   aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATAGCGACTTTAAAAATTAGAAAAATTATAAAGACGGAGATATTTCATATAATTCTGCGCTGAGGAAATCTACCCCTAGATATCAATTAACTAATGAAGAGta                                                                      >  2:536992/1‑140 (MQ=14)
                                                                    gaaaaGAAATCGAAAACTCTAAGTTATTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                                                     <  1:413848/139‑1 (MQ=255)
                                                                    aaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATAGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                                                     >  2:120006/1‑140 (MQ=255)
                                                                       aaGAAAACGAAAACTTTAAGTTCTTTGTTCAATACGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAAt                                                                  <  2:32326/140‑1 (MQ=255)
                                                                           aaTCGAACACGGTAAGTTTTTTGTTTAATATTGCGACTTTAAAAATTTGAAACACTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaa                                                              <  2:457811/140‑1 (MQ=17)
                                                                                  aaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATCCGGCTAAATATCAATTACCTAATGATGACTATAACGTAAAACactc                                                       >  1:663369/1‑137 (MQ=18)
                                                                                        aaGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaa                                                 <  1:662320/140‑1 (MQ=255)
                                                                                        aaGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaa                                                 <  2:290376/140‑1 (MQ=255)
                                                                                        aaGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaa                                                 <  1:17187/140‑1 (MQ=255)
                                                                                             cgttGTTCAATACGCCGATTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgat                                            <  2:733775/138‑1 (MQ=255)
                                                                                                   acAATATAGCGACTTTAAAAATTTAAAAAAATATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGACTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatatgata                                      <  2:470303/139‑1 (MQ=255)
                                                                                                    cAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatatgatat                                     >  2:585382/1‑140 (MQ=255)
                                                                                                    cAATATGGCGACTTTAAAAATTGGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatatgatat                                     >  2:275777/1‑140 (MQ=255)
                                                                                                      atatGGCGACTTTAAAAATTTGAAAAATTATACAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGACTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATAc                                   <  1:322779/140‑1 (MQ=255)
                                                                                                           gCGACTTTACAAAGTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGACGTGCCGAGTTATTCGGCTAAATATCAATTCACTCATGATGATTATAATGTAAAACAAGTACGCAAAAGATATGATATACcgacg                              <  1:205375/140‑1 (MQ=18)
                                                                                                           gCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACcgacg                              >  1:478688/1‑140 (MQ=255)
                                                                                                           gCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACagaag                              >  2:223383/1‑138 (MQ=255)
                                                                                                                tttAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTaa                         >  1:317875/1‑140 (MQ=255)
                                                                                                                     aaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTc                    <  1:636319/140‑1 (MQ=255)
                                                                                                                     aaaTTTAATAAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTc                    <  2:391627/140‑1 (MQ=255)
                                                                                                                      aattgtaaaaaGTATAAAGACGGAGATAGTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTcc                   <  1:54007/134‑1 (MQ=255)
                                                                                                                      aaattGAGAAATGATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTcc                   <  2:262295/137‑1 (MQ=255)
                                                                                                                      aaTCTGAAAACTTATAAAGACGGAGCTATTTCATATAACCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATAGACCGACGAGTAAAGCTcc                   <  1:36248/140‑1 (MQ=18)
                                                                                                                        tttGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCaa                 <  1:261157/140‑1 (MQ=255)
                                                                                                                            aaaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGtt             >  1:539252/1‑140 (MQ=255)
                                                                                                                             tgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTa            >  1:433290/3‑140 (MQ=255)
                                                                                                                               aataaacaaGATGGAAATAAATCATATAAACCACAGGGGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTAtt          <  2:457732/133‑1 (MQ=17)
                                                                                                                                aTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAACTATCAATTAACTAATGATGATTATACTGTAAAACAATTACGCAACAGATAGGATATACCGCCGAACAAAACCCCACAATTATTg         >  2:3335/1‑140 (MQ=11)
                                                                                                                                    taAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaa     <  1:241711/140‑1 (MQ=255)
                                                                                                                                    taAAGACGGAGAGATTTCATATAATCCAGAGGCGCCGAGTTATTCGGCCAAATAGCAAGTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaa     <  1:210008/140‑1 (MQ=18)
                                                                                                                                     aaaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTAGTGTTaaa    >  1:464434/1‑140 (MQ=255)
                                                                                                                                      aaGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaaaa   <  1:491139/140‑1 (MQ=255)
                                                                                                                                       aGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGAGTATAAGGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAGGCTCCAAAGGTATTGGTGAAAg  >  1:215817/1‑140 (MQ=17)
                                                                                                                                          |                                                                                                                                        
AAGATATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAG  >  USA300TCH1516_ALE/110513‑110787

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: