Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 111,512 A→G 30.2% S194S (TCA→TCG)  USA300HOU_RS00520 → putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE111,5120AG30.2% 97.5 / 48.5 66S194S (TCA→TCG) USA300HOU_RS00520putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (27/19);  new base G (7/13);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

CTACTGAGGAAGTTAAAGAAGAAAACATAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAA  >  USA300TCH1516_ALE/111373‑111650
                                                                                                                                            |                                                                                                                                          
cTACTGAGGAAGTTAAAGAAGAAAACATAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGCTATTCa                                                                                                                                            >  2:534864/1‑140 (MQ=255)
  aCTGAGGAAGTTAAAGAAGAAAACATAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGc                                                                                                                                          <  2:60026/140‑1 (MQ=255)
  aCTGAGGAAGTTAAAGAAGAAAACATAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGc                                                                                                                                          <  2:274767/140‑1 (MQ=255)
  aCTGAGGAAGTTAAAGAAGAAAACATAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAAAATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGc                                                                                                                                          <  1:267274/140‑1 (MQ=255)
    ggAGGAAGTTAAAGAAGAAAACATAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCta                                                                                                                                        <  2:22419/139‑2 (MQ=255)
         aaGTTACAGAAGAAAACATAAA‑AAACGCGCTCGAAAACCTCTTGTCACACACACGAGAGAGCAATTTTAAAGATTTAATCAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAatat                                                                                                                                   <  1:264173/140‑1 (MQ=9)
           gTTAAAGAAGAAAACATAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCa                                                                                                                                 <  1:186073/140‑1 (MQ=255)
            tcaaaGAAGAAAACGTAAA‑AAAAGAAATCGAACACTTTAAGTTTTTTGCGCAAGATAGCAATTTTAAAGAGTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCaa                                                                                                                                <  2:462002/138‑1 (MQ=18)
             tAAAGAAGAAAACATAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATCTAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAAt                                                                                                                               >  2:214123/1‑140 (MQ=255)
             tAAAGAAAAAAAAAAAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAataagatc                                                                                                                               <  2:550026/140‑6 (MQ=21)
              atagcagaaTACATAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAAtt                                                                                                                              <  2:381008/135‑1 (MQ=255)
               aagaagAAAACATAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATATCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTa                                                                                                                             <  2:92467/140‑1 (MQ=255)
                agaagaAAACATAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTaa                                                                                                                            >  2:205004/1‑140 (MQ=255)
                 gaagaaAACATAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAAc                                                                                                                           <  2:175927/140‑1 (MQ=255)
                    gaaAACATAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTAATCAGCGTAATATCAATTAACTaa                                                                                                                        >  2:131075/1‑140 (MQ=255)
                    cgaaaaATAAACAAAAGAAAGCGAAAACTTTAAGTTTTTTGCGCAATAAAGCAATTTTATAGATGTAATGAATTATAAAGGTGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTa                                                                                                                         <  2:545843/138‑1 (MQ=17)
                     aaaaactAAA‑AAAAGAAAGCGAAAACTTAAAGTTTTTTGCGCAATATAGCNATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAat                                                                                                                       <  2:140503/134‑1 (MQ=255)
                     aaaaCATAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAat                                                                                                                       >  1:401428/1‑140 (MQ=255)
                            aaa‑aaaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGAtt                                                                                                                <  2:242226/140‑1 (MQ=255)
                             aa‑aaaaGAAATCGAAAACTCTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATta                                                                                                               >  1:224248/1‑140 (MQ=255)
                              a‑aaaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                                                                                              >  1:408018/1‑140 (MQ=255)
                                   aGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATg                                                                                                          >  2:433165/1‑140 (MQ=255)
                                    gAAATTGAAAACTTTAAGTTTTTGGCGCAATATAGCAATTTTAAATATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGt                                                                                                         <  2:164254/140‑1 (MQ=255)
                                      agtCGAGAACTGCAAGGTTTTTGCGCAATATAGTAATTGTAAAGATTTAATGAATTATAAAGATGGGGATATATGATATGATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATTATTATAATGTaa                                                                                                       <  2:173698/138‑1 (MQ=17)
                                      aaTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGAAATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaa                                                                                                       <  1:147424/140‑1 (MQ=255)
                                       aTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaa                                                                                                      <  2:207953/140‑1 (MQ=255)
                                        gcGAAAACTTTAAGTGTTTTGCGCAATATAGCAATTTTAAAGAATTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaaa                                                                                                     <  2:69240/139‑1 (MQ=17)
                                          gAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACa                                                                                                   <  2:368429/140‑1 (MQ=255)
                                          gAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACa                                                                                                   <  1:80283/140‑1 (MQ=255)
                                            aaaCTTTAAGTTTTTTGCGCAATATAGTAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAt                                                                                                 <  2:101244/140‑1 (MQ=255)
                                             aaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAtt                                                                                                >  1:484481/1‑140 (MQ=255)
                                             aaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAtt                                                                                                >  2:419426/1‑140 (MQ=255)
                                             aaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAAtc                                                                                                >  1:85838/1‑139 (MQ=255)
                                              aCTTGAAGTTTTTAGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTa                                                                                               <  2:196022/140‑1 (MQ=255)
                                               cTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGAGTATAATGTAAAACAATTAc                                                                                              >  2:564909/1‑140 (MQ=255)
                                                 ttAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAACACAATTACgc                                                                                            >  1:63437/1‑139 (MQ=255)
                                                 ttAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACgc                                                                                            >  2:168039/1‑139 (MQ=255)
                                                   aaGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTCAACCAATTACGTac                                                                                          >  1:123702/1‑139 (MQ=255)
                                                        tgttGCGGGATAGAGCAATTTTAAAGATGTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGat                                                                                     <  2:547981/138‑1 (MQ=14)
                                                              gcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATAACa                                                                               >  2:47937/1‑140 (MQ=255)
                                                               cAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATACAGCGGTGCCGAGTTATTCGGCTAAATATCAATTAACTCAGGATGATTATAATGTAAAACAATTACGCAAAAgatatgatat                                                                              >  2:347031/1‑140 (MQ=14)
                                                                  tataGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATAcc                                                                           >  1:101335/1‑140 (MQ=14)
                                                                        aaTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACATAAAAAATATGACATAACAACaaa                                                                     >  2:273817/1‑140 (MQ=255)
                                                                         aTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAACGTAAATCAATTACGTAAAAGATATGACATACCAACAAAt                                                                    <  1:297807/140‑1 (MQ=255)
                                                                             tAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGACAATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAg                                                                >  1:565305/1‑140 (MQ=255)
                                                                                aGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgcgc                                                             >  1:17825/1‑140 (MQ=255)
                                                                                aGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgcgc                                                             >  1:318301/1‑140 (MQ=255)
                                                                                 gATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCCACAAATAAAGCGcc                                                            >  1:152633/1‑140 (MQ=255)
                                                                                       attaaTTACAAAAATGTATAAATATCATATAACCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCt                                                      <  2:524057/137‑1 (MQ=17)
                                                                                         gAATTATAAAGATGGAGATATATCATATAATACAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATAAGACATAAAAACAAACAAAGCACCGAAGCagt                                                    >  2:326779/1‑140 (MQ=17)
                                                                                                       gagaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCAGTTGTTGAAAGGTACa                                      >  2:266129/1‑140 (MQ=255)
                                                                                                        agaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAg                                     >  1:19536/1‑140 (MQ=255)
                                                                                                         gaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAgg                                    <  1:294527/140‑1 (MQ=255)
                                                                                                          atatatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCAGTTGTTGAAAGGTACAggg                                   >  2:21915/1‑140 (MQ=255)
                                                                                                          atatatCATATAAGCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAggg                                   <  1:87269/140‑1 (MQ=255)
                                                                                                                       tCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGtt                      >  1:404279/1‑140 (MQ=255)
                                                                                                                         gagagGTGCCGAGTTATTCAGATCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtca                    <  1:276543/139‑1 (MQ=255)
                                                                                                                           gagGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtcatc                  >  1:179620/1‑140 (MQ=255)
                                                                                                                            agGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtcatca                 >  1:42363/1‑140 (MQ=255)
                                                                                                                            aaacaccGAGTAATTCAGCTCAATATCAATTAACTAATGACGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtcatca                 <  1:80059/135‑1 (MQ=255)
                                                                                                                                gCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTg             >  2:572155/1‑140 (MQ=255)
                                                                                                                                accGAGTTAGTCAGCCCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTg             <  1:67787/139‑1 (MQ=255)
                                                                                                                                     gTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACACTTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGatat        >  1:198058/1‑140 (MQ=255)
                                                                                                                                     gTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGatat        <  2:501657/140‑1 (MQ=255)
                                                                                                                                     gTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGATACAGGGAATGTAAAAGGTTCATCAGTTGGatat        >  2:308058/1‑140 (MQ=255)
                                                                                                                                      ttATTCAGCTCAATATCAATTAACTAATGATGACTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGAtata       >  2:449546/1‑140 (MQ=255)
                                                                                                                                           cAGCTCAATATCAATTAACTAATGATGATTATAATGTAACACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCAACAGTGGGATGTaaaaaa  >  2:300219/1‑140 (MQ=255)
                                                                                                                                            |                                                                                                                                          
CTACTGAGGAAGTTAAAGAAGAAAACATAAA‑AAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAA  >  USA300TCH1516_ALE/111373‑111650

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