Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 110,651 C→T 14.0% 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,6510CT14.0% 153.0 / 11.0 67D170D (GAC→GATUSA300HOU_RS00515putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (28/28);  new base T (2/8);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 9.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTT  >  USA300TCH1516_ALE/110515‑110790
                                                                                                                                        |                                                                                                                                           
gATATTAGGTGAAAATGGACGAGAATAAAATCATTACAAGAGAAGAAATTTAACATTAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATGTGGCGACTTTAAAAATTTGAAAAATTATAAAGACGga                                                                                                                                          <  2:575314/140‑1 (MQ=17)
 atatCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATAGCGACTTTAAAAATTTGAAAAATTATAAAGACGgag                                                                                                                                         >  2:455204/1‑140 (MQ=255)
 atatCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGATCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGgag                                                                                                                                         >  2:125495/1‑140 (MQ=255)
  tatCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGgaga                                                                                                                                        <  2:430129/140‑1 (MQ=255)
    tCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGata                                                                                                                                      >  1:633657/1‑140 (MQ=255)
    tCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGata                                                                                                                                      >  1:331310/1‑140 (MQ=255)
     ccGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCAACTTTAAAAAAATGAAAAATTATAAAGACGGAGagat                                                                                                                                     >  2:171995/1‑140 (MQ=255)
      cGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATAtt                                                                                                                                    <  1:351749/140‑1 (MQ=255)
                   cGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCag                                                                                                                       <  1:382931/140‑1 (MQ=255)
                   cGATAATAAAAACATACCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCag                                                                                                                       <  1:560057/140‑1 (MQ=255)
                   cGAGAATAAAATCATACCAACAGAAGAAATGAGACATGAAAAAATACAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCag                                                                                                                       <  2:200410/140‑1 (MQ=17)
                       aataaAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGAAATTTCATATAATCCagagat                                                                                                                   >  2:557388/1‑138 (MQ=255)
                         taaAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATGAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGACAAATTATCAAGACGGAGATATTTCATATAATCCAGAGGTGc                                                                                                                 >  1:18133/1‑140 (MQ=255)
                         taaAATCATTCCAACAAAAGAAATTAAAGATGAAAAAAGAAAAAAAGAAATCGAAAACTTTAAGGTCTTGGGTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGc                                                                                                                 <  1:574892/140‑1 (MQ=255)
                           aaaTCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCg                                                                                                               <  1:589426/140‑1 (MQ=255)
                                 ttCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGAAGGAGATATCTAATATAATCCAAAGGTACCGaataat                                                                                                         >  2:59002/1‑135 (MQ=17)
                                   ccAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTc                                                                                                       <  2:297916/140‑1 (MQ=255)
                                            gAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCAACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGATGCCGAGTTATTCGGCTAAata                                                                                              >  2:685917/1‑140 (MQ=255)
                                               aTTAGAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCa                                                                                           <  2:20635/140‑1 (MQ=255)
                                                    aGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATCGCGACAAACAGGATGTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTaa                                                                                      <  1:697691/140‑1 (MQ=14)
                                                         aaaaaaTAAAAAAAGAAATCGAAAACTTTGAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAat                                                                                 <  2:450327/140‑1 (MQ=255)
                                                           aaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatga                                                                               <  2:331310/140‑1 (MQ=255)
                                                           aaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatga                                                                               <  2:633657/140‑1 (MQ=255)
                                                           aaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACAGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatga                                                                               >  2:247526/1‑140 (MQ=255)
                                                               aaaaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATgatgat                                                                           <  2:123092/139‑1 (MQ=255)
                                                                 aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGAttatta                                                                         >  1:480160/1‑140 (MQ=255)
                                                                 aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATta                                                                         >  1:603936/1‑140 (MQ=255)
                                                                 aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATta                                                                         >  2:166390/1‑140 (MQ=255)
                                                                  aaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCTGCTAAATATCAATTAACTAATGATGATtat                                                                        <  2:216822/140‑1 (MQ=255)
                                                                   aaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCATAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtata                                                                       >  2:512805/1‑140 (MQ=255)
                                                                          aTCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaa                                                                >  2:494772/1‑140 (MQ=255)
                                                                           tCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaaa                                                               >  2:551879/1‑140 (MQ=255)
                                                                           tCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaaa                                                               >  1:472098/1‑140 (MQ=255)
                                                                               aaaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                           <  1:112266/140‑1 (MQ=255)
                                                                               aaaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATACTTCATATAATCCAGAGGTGCCGAGCTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                           <  2:351482/140‑1 (MQ=255)
                                                                                acatttAAGTTATTTGTTCAATATGGGGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAGCAAtt                                                          <  2:229559/137‑1 (MQ=255)
                                                                                   tttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGCTATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGACGATTATACTGTAAAACAATTACg                                                       >  2:272063/1‑140 (MQ=255)
                                                                                   tttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAGAACAATTACg                                                       >  2:635719/1‑140 (MQ=255)
                                                                                    ttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc                                                      <  2:367149/140‑1 (MQ=255)
                                                                                      aaGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTACGTATCAATTAACTAATGATGATTATAATGTATCACAATTTCGCat                                                    >  1:360326/1‑139 (MQ=18)
                                                                                               gTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatatga                                           >  1:176486/1‑140 (MQ=255)
                                                                                                           gACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAg                               <  1:581843/140‑1 (MQ=255)
                                                                                                              tttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTaa                            <  2:326765/140‑1 (MQ=11)
                                                                                                              tttAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTACATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTaa                            <  2:47169/140‑1 (MQ=255)
                                                                                                              tttAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAACATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTaa                            <  2:673157/140‑1 (MQ=255)
                                                                                                                 aaaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGc                         <  2:65547/140‑1 (MQ=255)
                                                                                                                  aaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAGAGCt                        >  1:121044/1‑140 (MQ=255)
                                                                                                                   aaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTc                       >  2:671701/1‑140 (MQ=255)
                                                                                                                    aaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCACAAGATATGATATAACGACGAGTAAAGcccc                      >  2:139458/1‑140 (MQ=255)
                                                                                                                     aTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCa                     <  1:618312/140‑1 (MQ=255)
                                                                                                                       ttGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAAATATGATATACCGACGAGTAAAGCTCCaaa                   >  2:262967/1‑140 (MQ=255)
                                                                                                                        tGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTACCGGTATCGCCGTCTTAATATAAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCCAAg                  >  1:194600/1‑140 (MQ=18)
                                                                                                                         gatgaagtataAAGATGGAGATATATCATATAAGCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCCCCAAAGt                 <  2:305259/133‑1 (MQ=14)
                                                                                                                         aatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGt                 <  2:106459/136‑1 (MQ=255)
                                                                                                                         aatgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGt                 <  2:556564/136‑1 (MQ=255)
                                                                                                                          aaaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGtt                <  1:531901/140‑1 (MQ=255)
                                                                                                                           aaaaTTATAAAGACGGAGATATTTCATATAATCCAGAAGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATCTGATATACCGACGAGTAAAGCTCCACAGTTa               >  2:702216/1‑140 (MQ=255)
                                                                                                                           aaaaTTATAAAGACGGAGATATTTCATATAATACAGAGGTGCCAAGCTATTCGGCTGCATATCTATTAACTCATGATGATTATAATGTAAAACAATTACGCAACAGCTATGCTATAACGCCGATTACAGATCCAAAgctg               >  2:318181/1‑137 (MQ=17)
                                                                                                                            aaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAAGAAAGCTCCAAACTTAt              >  2:144720/1‑140 (MQ=255)
                                                                                                                              aTTATAGAGATGGAGATATATCATATAATCCAGAGGTGTCGAGCTAGTCGTCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTg            >  1:102559/1‑140 (MQ=17)
                                                                                                                              aTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTg            >  1:539115/1‑140 (MQ=255)
                                                                                                                                tgtcaaGACGGAGATATTTCATAGAATCAAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGtt          <  2:409377/136‑1 (MQ=255)
                                                                                                                                  tagaGATGGAGATATGTCATATAATCAAGAGGTGCCGAGTGATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaa        <  2:480160/137‑1 (MQ=18)
                                                                                                                                  taAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaa        <  1:143340/140‑1 (MQ=255)
                                                                                                                                  taAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaa        <  1:533963/140‑1 (MQ=255)
                                                                                                                                  taAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaa        <  1:166390/140‑1 (MQ=255)
                                                                                                                                   aaaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaaa       <  2:538461/140‑1 (MQ=255)
                                                                                                                                        cGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAGGCTCCAAAGTTATTGGTAAAAGGtg  >  1:182051/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                           
GATATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTT  >  USA300TCH1516_ALE/110515‑110790

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: