Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 110,660 T→A 18.9% 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,6600TA18.9% 155.0 / 28.5 68I173I (ATT→ATA) USA300HOU_RS00515putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (17/37);  new base A (6/7);  total (23/45)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGG  >  USA300TCH1516_ALE/110523‑110794
                                                                                                                                         |                                                                                                                                      
gTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTATAAAAATTTGAAAAATCATAAAGACGGAGATATATCa                                                                                                                                     >  2:318030/1‑141 (MQ=255)
 ttAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCat                                                                                                                                    <  1:631664/141‑1 (MQ=255)
   aaaaTGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAAAAATCGAAAACTTTAAGTTCTTTGTTCAATATGACGACTTTAAAAATTTGAAAAATTATAAAAACGGAGATATTTGatat                                                                                                                                  >  1:157816/1‑141 (MQ=255)
     aaTGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATaa                                                                                                                                <  2:465388/141‑1 (MQ=255)
          tCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCag                                                                                                                           <  2:172217/141‑1 (MQ=255)
            gATAATAAAATCATTCCAACAACCGATATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATACaaag                                                                                                                         >  1:185219/1‑141 (MQ=17)
             ataataAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAgg                                                                                                                        <  1:625174/141‑1 (MQ=255)
             ataataAAATCACTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAgg                                                                                                                        <  2:457780/141‑1 (MQ=255)
               aataaAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAACGAAATCAAAAACCGTAGTTACTTTATTCAATAGGGCAAATTTAAAAATATAAAACATTATCAAAACAGAGAAATTTCAAATATTCCaacagta                                                                                                                      >  2:505875/1‑135 (MQ=12)
                   aaaTCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGa                                                                                                                  <  1:607832/141‑1 (MQ=255)
                   aaaTCATTCCAAAAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGa                                                                                                                  <  2:34218/141‑1 (MQ=255)
                    aaccATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATGTTTCATATAATCCAGAGGTTCCGAg                                                                                                                 <  1:658011/138‑1 (MQ=255)
                    aaacATTACAACAAAAGAAATTAAAGCTAAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAg                                                                                                                 <  2:606406/138‑1 (MQ=255)
                    aaTCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAg                                                                                                                 <  2:716725/141‑1 (MQ=255)
                     aTCATGACAACGAAAGAGAAGAAAGAGGAAAAAATAAAAAAAGAAAACGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAAAAAAAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGt                                                                                                                <  2:294761/141‑1 (MQ=17)
                        aTTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTAt                                                                                                             <  1:708679/141‑1 (MQ=255)
                        aTTCCAACAAAAGAAAGTAAAGATGAAAAAATAAAAAAAGAAATCGAAATCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTAt                                                                                                             <  2:1686/141‑1 (MQ=255)
                         ttCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAATTTATAAAGACGGAGATATTTAATATAATCCAGAGGTGCCGAGTTAtt                                                                                                            >  2:11375/1‑141 (MQ=255)
                         ttCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTAtt                                                                                                            <  1:325976/141‑1 (MQ=255)
                           ccAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCg                                                                                                          <  1:35124/141‑1 (MQ=255)
                           ccAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAAGATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCg                                                                                                          <  2:578489/141‑1 (MQ=255)
                           acaacaaAAAAAATAAACAAAGGAAAAAGAAAAAAACAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCTACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCg                                                                                                          <  2:741539/140‑1 (MQ=17)
                            caacaaAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCgg                                                                                                         <  1:620726/141‑1 (MQ=255)
                             aacaaAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGc                                                                                                        >  1:74493/1‑141 (MQ=255)
                              acaaAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGTTATTTCATATAATCCAGAGGTGCATAGTTATTAGGCt                                                                                                       >  1:49974/1‑141 (MQ=255)
                               aaaaaGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAACAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTa                                                                                                      <  2:586775/140‑1 (MQ=255)
                                     aaaTTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTGGGCTAAATATc                                                                                                >  2:369513/1‑141 (MQ=255)
                                         tAAAGTAGAAAAAATAAAAAACGGAGTTGAAAACTTTATGTTTTATGTCCAGTATGGCGACTTTAGAAATGTGAAAAATTATAAAGACGGATATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAtt                                                                                            <  2:395714/141‑1 (MQ=17)
                                         tAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTAAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGATAAATATCAAtt                                                                                            >  1:742588/1‑141 (MQ=255)
                                            aGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGACAAATAAGAGATACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAAc                                                                                         >  1:368898/1‑141 (MQ=17)
                                            aGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCGGAGTTATTCGGCTAAATATAAATTAAc                                                                                         >  2:10792/1‑141 (MQ=255)
                                            aGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAAc                                                                                         <  1:264701/141‑1 (MQ=255)
                                                gAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAat                                                                                     <  1:67965/141‑1 (MQ=255)
                                                aaaaaaaTAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAat                                                                                     <  2:385420/140‑1 (MQ=255)
                                                       tAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAAtgatgact                                                                              >  2:274334/1‑139 (MQ=255)
                                                        aaaaaaaGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATta                                                                             >  2:777127/1‑141 (MQ=255)
                                                          aaaaaGAAATCGAAAACTTTACGTTCTTTGTTCAATAGGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtata                                                                           <  1:102317/141‑1 (MQ=255)
                                                                aaaTCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaa                                                                     <  1:286613/141‑1 (MQ=255)
                                                                      aaaaCTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAGATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                               <  1:431628/141‑1 (MQ=255)
                                                                         aCTTTAAGGTGTTTGTTCAATATGGGGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTAc                                                            <  2:143266/141‑1 (MQ=255)
                                                                          cTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACg                                                           <  1:479887/141‑1 (MQ=255)
                                                                          cTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACg                                                           <  1:416812/141‑1 (MQ=255)
                                                                           tttAAGTTCTTTGTTCAATATGGTGACTGTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGc                                                          <  2:626870/141‑1 (MQ=255)
                                                                             tAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaa                                                        <  1:250534/141‑1 (MQ=255)
                                                                              aaGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaa                                                       <  2:738181/141‑1 (MQ=255)
                                                                                          cAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGAAATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGAtata                                           >  2:413332/1‑141 (MQ=255)
                                                                                                 gCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACcgacga                                    >  2:116501/1‑141 (MQ=255)
                                                                                                     ttttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTaa                                <  2:620727/140‑1 (MQ=11)
                                                                                                     ttttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTaa                                <  1:236000/140‑1 (MQ=11)
                                                                                                        tAAAAAATTGAAAAATTATAAAGGCGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGc                             <  2:257767/141‑1 (MQ=255)
                                                                                                          aaaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTc                           >  2:668436/1‑141 (MQ=255)
                                                                                                          aaaaTTTGAAAAATTACAAAGACGGAGATATTTCATAGAATCCAGAGGTGCCGAGTTATGCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGCTATGATATACCGACGAGAAAAGCTc                           >  1:383742/1‑141 (MQ=18)
                                                                                                           aaaTTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTcc                          <  2:77328/141‑1 (MQ=255)
                                                                                                               ttGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAg                      >  1:285861/1‑141 (MQ=255)
                                                                                                                   tgaaTTATAAAGATGGGGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTAt                  <  1:704395/139‑1 (MQ=255)
                                                                                                                   tgaaTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTAt                  >  2:458103/3‑141 (MQ=255)
                                                                                                                   agaaTTATAAAGATGGAGATATGTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTAt                  <  2:584391/139‑1 (MQ=255)
                                                                                                                          taAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaaa           <  1:314109/141‑1 (MQ=255)
                                                                                                                          taAAGACGGAGACATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaaa           <  1:620016/141‑1 (MQ=255)
                                                                                                                          aaaaGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaaa           <  1:287323/140‑1 (MQ=255)
                                                                                                                           aaaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTaaaa          >  1:430759/1‑141 (MQ=255)
                                                                                                                             aGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGCCGAGAAAATCTCCAAAGTTATTGTTAAAAgg        <  2:166374/141‑1 (MQ=17)
                                                                                                                             aGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAgg        <  1:63315/141‑1 (MQ=255)
                                                                                                                              gatggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGt       >  1:217472/4‑141 (MQ=255)
                                                                                                                               atggAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCCAAACACTGGAAAGAGCGACGAGTAAAGCTCCAAAGTTATTGGTAAAAGGtt      >  1:401701/3‑141 (MQ=17)
                                                                                                                               atggAAATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGtt      <  1:495174/139‑1 (MQ=255)
                                                                                                                                  gagaTATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTAATTGTTAAAAGGTTCAg   >  1:388679/1‑141 (MQ=255)
                                                                                                                                   agaTATTTCATATAAGCCAGAGGTGCCGAGTTATTCGTCTAAAGATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATCATATACCGACGAGTAAAGCTCCAAAGTTATTGTCAAAAGGTTCAgg  <  1:547376/141‑1 (MQ=21)
                                                                                                                                         |                                                                                                                                      
GTTAAAATGGTCGATAATAAAATCATTCCAACAAAAGAAATTAAAGATGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTCTTTGTTCAATATGGCGACTTTAAAAATTTGAAAAATTATAAAGACGGAGATATTTCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACCGACGAGTAAAGCTCCAAAGTTATTGTTAAAAGGTTCAGG  >  USA300TCH1516_ALE/110523‑110794

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