Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 110,660 T→A 11.4% I173I (ATT→ATA)  USA300HOU_RS00515 → putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE110,6600TA11.4% 184.1 / 13.2 70I173I (ATT→ATA) USA300HOU_RS00515putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (23/39);  new base A (3/5);  total (26/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

GGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGAA  >  USA300TCH1516_ALE/110522‑110797
                                                                                                                                          |                                                                                                                                         
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                                                           <  2:290376/140‑1 (MQ=255)
                                                                               aaGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaa                                                           <  1:17187/140‑1 (MQ=255)
                                                                               aaGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaa                                                           <  1:662320/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)
                                                                                                                               gatggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAgg           <  2:748415/137‑1 (MQ=255)
                                                                                                                                  gaaaatatatCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTc        <  2:624559/136‑1 (MQ=255)
                                                                                                                                   gagaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAACAGGTTCa       <  2:167730/140‑1 (MQ=255)
                                                                                                                                      atatatCATATAATCCAGAGGTGCAGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAggg    >  1:570142/1‑140 (MQ=255)
                                                                                                                                      atatTTCATATAATCCAGAGTTGGCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACGATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAggg    <  2:196420/140‑1 (MQ=255)
                                                                                                                                      atatTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAggg    <  1:585382/140‑1 (MQ=255)
                                                                                                                                        atTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGaa  >  1:60037/1‑140 (MQ=255)
                                                                                                                                        atTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATAAGCGAATATAGCTTAAGATTTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGCCATTGTTAAAAGGTCCAGGGaa  >  2:342879/1‑140 (MQ=18)
                                                                                                                                          |                                                                                                                                         
GGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGAA  >  USA300TCH1516_ALE/110522‑110797

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