Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 110,651 C→T 8.5% 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.5% 164.8 / 5.9 64D170D (GAC→GAT) USA300HOU_RS00515putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (26/31);  new base T (3/2);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

AAAAGATATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGG  >  USA300TCH1516_ALE/110511‑110788
                                                                                                                                            |                                                                                                                                          
aaaaGATATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCCATATGGCGCCTTTACACATTTGAAAAATTATAAacac                                                                                                                                            >  1:666632/1‑141 (MQ=18)
  aaGATATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGAcaa                                                                                                                                          >  1:218679/1‑139 (MQ=255)
   aGATATTCGGCAAAAATGGAATATAAGAAAATCATTCGAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGga                                                                                                                                         <  1:183955/141‑1 (MQ=17)
    gATATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGga‑g                                                                                                                                       <  1:360271/141‑1 (MQ=255)
    aatatCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGga‑g                                                                                                                                       <  1:491708/140‑1 (MQ=255)
    aatatCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGAACAAATTCTGAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGga‑g                                                                                                                                       <  1:798699/140‑1 (MQ=18)
     aaatgccgTTAAAATGTACAATCATAAAATAATCCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTCTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGga‑ga                                                                                                                                      <  1:143854/134‑1 (MQ=17)
      tatCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACgaggat                                                                                                                                      >  1:697695/1‑136 (MQ=255)
      tatCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAAACGGA‑Gat                                                                                                                                     >  1:50554/1‑141 (MQ=255)
          cGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGTCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATAttt                                                                                                                                 <  1:804848/141‑1 (MQ=255)
           ggTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAAGAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTc                                                                                                                                <  1:286213/141‑1 (MQ=255)
            gTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCa                                                                                                                               <  1:3740/141‑1 (MQ=255)
               aaaaTGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCatat                                                                                                                            <  1:783006/141‑1 (MQ=255)
                aaaTGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCAtata                                                                                                                           >  1:372474/1‑141 (MQ=255)
                              aaaaTCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAATAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATTCAGAGGTGcca                                                                                                             >  1:78546/1‑140 (MQ=255)
                                      tCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTc                                                                                                     <  1:135152/141‑1 (MQ=255)
                                       ccAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGGGGTGCCGAGTTATTCg                                                                                                    <  1:405277/141‑1 (MQ=255)
                                        caacaaAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTGAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCgg                                                                                                   >  1:550318/1‑141 (MQ=255)
                                               aGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAata                                                                                            <  1:35410/141‑1 (MQ=255)
                                               aGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAata                                                                                            <  1:632180/141‑1 (MQ=255)
                                                 aaaTTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTGCTGTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAatata                                                                                          >  1:697792/1‑140 (MQ=255)
                                                     tAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAtt                                                                                      >  1:500360/1‑141 (MQ=255)
                                                     tAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAtt                                                                                      >  1:278024/1‑141 (MQ=255)
                                                     gaaaGATGAAAAAATAAAAAAAGAAATCGAAAACACTAAGTTCTTGGATCAATATGACGACATTAAAAAATGGAAAAAGTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCAGAGTTATAAGGCTAAATATCAAtt                                                                                      >  1:480818/2‑141 (MQ=17)
                                                       aaGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTaa                                                                                    >  1:577201/1‑141 (MQ=255)
                                                       aaGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTaa                                                                                    >  1:603199/1‑141 (MQ=255)
                                                       aaGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTaa                                                                                    >  1:639762/1‑141 (MQ=255)
                                                        aGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTACCGAGTTATTCGGCTAAATATCAATTAAc                                                                                   >  1:368303/1‑141 (MQ=255)
                                                         gATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATGTAATGAAGAGGTGCCGGGTTATTCGGCTAAATATGAATTaatt                                                                                  >  1:13240/1‑139 (MQ=17)
                                                                     aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTACATATCAATTAACTAATGATGATtaa                                                                      >  1:41320/1‑140 (MQ=255)
                                                                     aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                                                      >  1:228997/1‑141 (MQ=255)
                                                                     aaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                                                      <  1:411945/141‑1 (MQ=255)
                                                                                cGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACa                                                           >  1:6300/1‑141 (MQ=255)
                                                                                cGAAAACTTTAAGTTCTTTCTTCAATATGGCGACTTTAAAAATTTGACAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTCACTAATGATGATTATAATGTacaaca                                                           >  1:487995/1‑141 (MQ=255)
                                                                                cGAAAACTTCAAGTTCTGTTTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACa                                                           <  1:644625/141‑1 (MQ=255)
                                                                                  aacattttAAGTTCTTTGTTGAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                         <  1:95380/136‑1 (MQ=255)
                                                                                  aaaaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATCTAGAGGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAAGGATGATTATAATGTAAAACAAt                                                         <  1:388088/141‑1 (MQ=255)
                                                                                  aaaaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                         <  1:251179/141‑1 (MQ=255)
                                                                                       tttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATTGCATCCAGAGGTGCCTAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc                                                    <  1:595544/141‑1 (MQ=255)
                                                                                        ttAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTATCTCAAGTGAAGTAGAATGTAAAACAATTACGCa                                                   <  1:805693/141‑1 (MQ=17)
                                                                                         tAAGTTCTTTGTTCAATATGGCGCCTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCTAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAACTACGCaa                                                  <  1:513665/141‑1 (MQ=255)
                                                                                         tAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaa                                                  <  1:315856/141‑1 (MQ=255)
                                                                                         tAAGGTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaa                                                  <  1:787066/141‑1 (MQ=255)
                                                                                              tCTTTGTTCAATATTGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgat                                             <  1:580832/141‑1 (MQ=255)
                                                                                                    ttCAATATGGCGCATTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatatgata                                       <  1:56840/141‑1 (MQ=255)
                                                                                                      cAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTCATGATGATTATAATGTAGAACAATTACGCAAAAGATATGAtgta                                     >  1:165707/1‑141 (MQ=255)
                                                                                                      cAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTAAGCAAAAGATATGAtata                                     >  1:485106/1‑141 (MQ=255)
                                                                                                          atGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAGCAATTACGCAAAAGATATGATATACcga                                 >  1:135790/1‑141 (MQ=255)
                                                                                                               gACTTTAAAAATTTGAAAAATTTTAAAGACGGA‑TATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGt                            >  1:396661/1‑141 (MQ=255)
                                                                                                                 cTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTaa                          <  1:645632/141‑1 (MQ=255)
                                                                                                                   ttAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAg                        >  1:533351/1‑141 (MQ=255)
                                                                                                                     aaaaaTTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCt                      >  1:378506/1‑141 (MQ=255)
                                                                                                                     aaaaaTTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCt                      >  1:270059/1‑141 (MQ=255)
                                                                                                                     aaaaaTTTGAAAAATTATAAA‑ATGAATTAGATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCt                      >  1:195917/1‑141 (MQ=17)
                                                                                                                          tttGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCaaa                 <  1:356056/141‑1 (MQ=255)
                                                                                                                             gAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGACTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGAAGAGTAAAGGTCCAAAGtt              >  1:286032/1‑141 (MQ=255)
                                                                                                                              atgaaTGATAAAGAATGA‑GATATAACATAAAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTa             <  1:603771/138‑1 (MQ=17)
                                                                                                                              aaaaaTTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTa             <  1:32150/141‑1 (MQ=255)
                                                                                                                              aaaaaTTATAAAGACGGA‑GATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATGAACTAATGATCATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTa             <  1:584258/141‑1 (MQ=255)
                                                                                                                                  aTTATAAAGATGGA‑GATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGt         <  1:360717/141‑1 (MQ=255)
                                                                                                                                  aTTATAAAGATGGA‑GATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGt         <  1:810814/141‑1 (MQ=255)
                                                                                                                                      taAAGACGGA‑GATATTTCATACAATCCAGAGGTTCCGCGTTTCTGAAATCGCAATGGAGTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaaa     <  1:588702/141‑1 (MQ=18)
                                                                                                                                       aaaGACGTA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaaaa    <  1:397579/141‑1 (MQ=255)
                                                                                                                                        aaGATGGA‑GATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAg   >  1:800053/1‑141 (MQ=255)
                                                                                                                                        aaGATGGA‑GATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATGGTTAAAAg   >  1:472422/1‑141 (MQ=255)
                                                                                                                                         aGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAgg  <  1:286864/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                          
AAAAGATATCCGGTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGA‑GATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGG  >  USA300TCH1516_ALE/110511‑110788

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑