Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 110,660 T→A 22.3% 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,6600TA22.3% 135.3 / 36.5 71I173I (ATT→ATA) USA300HOU_RS00515putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (20/35);  new base A (9/7);  total (29/42)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGAA  >  USA300TCH1516_ALE/110521‑110797
                                                                                                                                           |                                                                                                                                         
cGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATAttt                                                                                                                                          <  1:318086/141‑1 (MQ=255)
   ttAAAATGGTCGAAAATAAAATCATACCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCat                                                                                                                                       <  2:177719/141‑1 (MQ=255)
    tAAAACGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGAAGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCata                                                                                                                                      <  2:602119/141‑1 (MQ=255)
    tAAAACAGTCAATAATAAAACTATTCCAAAGAATAAAATTAAAGATTAACAAATAATAAAAGAAATCTAAAACTTTAAGTTCGTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCata                                                                                                                                      <  2:704983/141‑1 (MQ=17)
       aaTGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGACGCAAAAATAAACACAGACATTGAAAAATAGAACCTCTTTGTTCAAGATAGCGACTTGAAAGATTGGAAAAGTTATAAAGACGGAGATATTTCATATaa                                                                                                                                   >  2:765191/1‑141 (MQ=14)
         tGGTCGATAATAAAATCATTCCAATAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATc                                                                                                                                 >  1:494431/1‑141 (MQ=255)
            tcaataataAAAACATTCCAACAAAAGACATTAAAGATAAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAAGATTGCGACTTTAAAAATTTGAAAAATTATAAAGACGAAGATATATTATATAATCCag                                                                                                                              <  2:51342/138‑1 (MQ=17)
                taataaAATCATTCCAACAAAAGAAGTAAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCATTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAAACAGAGATATTTAATATAATCCaaaaaa                                                                                                                          >  1:394735/1‑136 (MQ=17)
                        ccATTCCAACAAAAAAAATCAAAAAACAAAAAATAAAAAAAGAAATCGAAAACCTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGtt                                                                                                                  <  2:599684/140‑1 (MQ=17)
                               aacaaAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGc                                                                                                           <  2:147373/141‑1 (MQ=255)
                                acgaaaGAAATGAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCt                                                                                                          <  1:247110/138‑1 (MQ=255)
                                acaaAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCt                                                                                                          <  2:232181/141‑1 (MQ=255)
                                acaaAAGAAAACAAATATGTAAAAAGAAAAAAAGAGATCGAAAACTTTAAGGTCTGTGTTCAATATGGCGACTTTAAAAATATGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGACATATACGGCt                                                                                                          >  2:184260/1‑141 (MQ=17)
                                 caaAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTa                                                                                                         >  2:24268/1‑141 (MQ=255)
                                  aaaaGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTaa                                                                                                        <  1:465006/141‑1 (MQ=255)
                                    aaGAAATTAAAGAGGAAAAAATAAAAAAACAAATCCAAAACTTTAAATTCTTTGTTCAATATGGCGAATAAAAAAAAAAAAAAAATTATAAAGACGGAGATATTTCATATAATCCAGTGGTGCCGAGTTATTCGGCTAAat                                                                                                      <  2:748688/141‑1 (MQ=17)
                                              aGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAAc                                                                                            <  2:78005/141‑1 (MQ=255)
                                              aGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAAc                                                                                            <  1:26001/141‑1 (MQ=255)
                                              aGATAAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAAc                                                                                            <  2:220608/141‑1 (MQ=255)
                                               gATGAAAAAATAAAAAAAGAAATCGAAAAATTAGAGTTCTTTGATAAATAGGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAAACAAAAGGTGCAGAGTTATTCGGCTAAATATAAATTAAca                                                                                           >  2:336073/1‑140 (MQ=17)
                                                aagcaacaaTAAAAAATGCAACCGAAAATATGAAGTTCTTTGTGCAATATGGCGACTGTAAAAATTTGAAAAATTATAAATACAGAGATATTTCAGATAATCCATAGGGGCCGAGTTATTCGGCTAAATATCAATTAACTa                                                                                          <  2:640528/134‑1 (MQ=17)
                                                  gAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAat                                                                                        <  1:311789/141‑1 (MQ=255)
                                                   aaaaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGCGGTGCCGAGTTGTTCGGCTAAATATCAATTCACTAatg                                                                                       >  2:650898/1‑141 (MQ=255)
                                                      aaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgatg                                                                                    <  1:444125/141‑1 (MQ=255)
                                                      aaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgatg                                                                                    <  1:759245/141‑1 (MQ=255)
                                                        aTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGACAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATTGCAATCAACTAATgaggat                                                                                  >  2:172766/1‑141 (MQ=14)
                                                             aaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTATCTAATGATGATTATaa                                                                             <  1:383059/141‑1 (MQ=255)
                                                                 gAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTa                                                                         >  2:758350/1‑141 (MQ=255)
                                                                 aaaaTCTAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTa                                                                         <  1:390983/140‑1 (MQ=255)
                                                                         aaaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAtt                                                                 >  1:140765/1‑141 (MQ=255)
                                                                         aaaCTTTAAGTTCTTAGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAtt                                                                 >  1:575525/1‑141 (MQ=255)
                                                                            cTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTAAAAAGACGAAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATGATAATGTAAAACAATTACg                                                              >  1:644359/1‑141 (MQ=255)
                                                                             tttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc                                                             <  1:701788/141‑1 (MQ=255)
                                                                             tttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc                                                             <  2:587219/141‑1 (MQ=255)
                                                                                    tCTGTGTTCAATATGGCGACTTTAAAAATTCACTAGATTTGAAAGACTTGTATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATCATTATAATCTAACCCCATTACGCACAAcat                                                      <  2:542759/141‑1 (MQ=9)
                                                                                             aaTATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATAc                                             >  1:314700/1‑141 (MQ=255)
                                                                                                   gCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACcgacga                                       >  1:68929/1‑141 (MQ=255)
                                                                                                    cGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAg                                      <  1:766676/141‑1 (MQ=255)
                                                                                                    cGACTTTAAAAATTTGAAAAAGTATAAAGACGTAGATATTTCATATAATCCAGAGGGGCCGAGTTATTCTGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAg                                      <  1:738530/141‑1 (MQ=255)
                                                                                                       cTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTaa                                   <  1:732286/141‑1 (MQ=255)
                                                                                                        tttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTaaa                                  <  2:633037/141‑1 (MQ=11)
                                                                                                          tAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGc                                <  2:394735/141‑1 (MQ=255)
                                                                                                            aaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTc                              >  2:361640/1‑141 (MQ=255)
                                                                                                             aaaaTTGAAAAATAATAAAAACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTcc                             <  2:328481/141‑1 (MQ=255)
                                                                                                             aaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTcc                             <  2:769101/141‑1 (MQ=255)
                                                                                                                attaatgaaTTAAAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAAAGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCaaa                          <  2:48219/134‑1 (MQ=14)
                                                                                                                 ttGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATCTACCGACGAGTAAAGCTCCAAAg                         >  2:191653/1‑141 (MQ=255)
                                                                                                                    aaaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTa                      <  2:42796/141‑1 (MQ=255)
                                                                                                                    aaaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTa                      <  2:212178/141‑1 (MQ=255)
                                                                                                                    aaaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTa                      <  2:597467/141‑1 (MQ=255)
                                                                                                                     tgaaTTATAAAGATGGAGCTCTATCATATAATCAAGAGGTGCCGAGGTATTCGGCGAAATTTCAATTAACTAATGAGGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTAt                     >  2:534705/3‑141 (MQ=17)
                                                                                                                     tgaaTTATAAAGATGGAGATATATCATATAATCAAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGAAAAAGATATGATATACCGACGAGTAAAGCTACAAAGTTAt                     >  2:708309/3‑141 (MQ=14)
                                                                                                                     aaaaTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTTTTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTAt                     <  1:59804/141‑1 (MQ=255)
                                                                                                                        aTTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCAACGAGTAAAGCTCCAAAGGTaaggt                  >  2:53968/1‑137 (MQ=255)
                                                                                                                         ttATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGtt                 <  2:289029/141‑1 (MQ=255)
                                                                                                                               aGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAgg           <  2:428757/141‑1 (MQ=255)
                                                                                                                               aGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAgg           <  1:260646/141‑1 (MQ=255)
                                                                                                                               aGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAgg           <  2:33456/141‑1 (MQ=255)
                                                                                                                                gatggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGt          >  1:119327/4‑141 (MQ=255)
                                                                                                                                gACGGAGATATTTCATATAATCCAGAGGTGCCGAGTAATTCGGCTAAATATCAATTAATTAATGCTGATAAAAAGGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGATCCAAAGTTATTGTTAAAAGGt          >  1:534965/1‑141 (MQ=18)
                                                                                                                                 atggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGtt         <  2:333820/139‑1 (MQ=255)
                                                                                                                                 aCGGAGATATTTCATATAATCCAGCGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGtt         <  2:179290/141‑1 (MQ=255)
                                                                                                                                   ggAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATAATATACCGACGAGGAAAGCTCCAAAGTTATTGTTAAAAGGTACa       >  2:53045/1‑141 (MQ=255)
                                                                                                                                     agaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCGATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAgg     >  2:500541/1‑141 (MQ=255)
                                                                                                                                     agaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAgg     >  2:437673/1‑141 (MQ=255)
                                                                                                                                     agaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGAAAAAGATATGATATACCAACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCggg     >  2:398789/1‑141 (MQ=255)
                                                                                                                                      gaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCTAAAGTTATTGTTAAAAGGTTCAggg    >  1:245451/1‑141 (MQ=255)
                                                                                                                                      gaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCGCTAAAGTGATTGTTAAAAGGTTCAggg    >  1:245459/1‑141 (MQ=255)
                                                                                                                                      gaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCCCCAAAGTTATTGTTAAAAGGTTCAggg    >  2:208852/1‑141 (MQ=255)
                                                                                                                                       atatTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCCACGTTATTGTTAAAAGGTTCAGGGa   >  1:662513/1‑141 (MQ=255)
                                                                                                                                        tatTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGGTAAAAGGTTCAGGGaa  <  1:382626/141‑1 (MQ=255)
                                                                                                                                        tatTCCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGCTAAAAGGTTCAGGGaa  <  1:382646/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                         
CGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGGGAA  >  USA300TCH1516_ALE/110521‑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: ‑