Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE111,5120AG26.4% 120.9 / 37.0 73S194S (TCA→TCG) USA300HOU_RS00520putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (28/25);  new base G (9/10);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 7.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATA  >  USA300TCH1516_ALE/111380‑111645
                                                                                                                                    |                                                                                                                                     
ggAAGTTAAAGAAGAAAACATATAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCTGc                                                                                                                                     <  1:364947/135‑1 (MQ=255)
 gAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCt                                                                                                                                    <  1:32696/135‑1 (MQ=255)
  aaGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTc                                                                                                                                   >  2:554820/1‑135 (MQ=255)
  aaGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGCTGCCGAGTTATTCGGCta                                                                                                                                   >  2:20297/1‑134 (MQ=255)
         agaagaAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATAACa                                                                                                                            <  1:363862/135‑1 (MQ=255)
          gaagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATGTATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCaa                                                                                                                           >  2:585204/1‑135 (MQ=255)
           aagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAAt                                                                                                                          <  2:321854/135‑1 (MQ=255)
             gcaaaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTa                                                                                                                        >  1:120219/3‑135 (MQ=255)
                      aaaaaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatgat                                                                                                               >  1:638709/1‑135 (MQ=255)
                       aaaaaGAAATCGAAAACTTTAAGTTTTTTGCGCATTATAGCTATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgatg                                                                                                              <  2:397085/135‑1 (MQ=255)
                       aaaaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgatg                                                                                                              <  1:192779/135‑1 (MQ=255)
                        aaaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTTAATATCAATTAACTAAtgatga                                                                                                             <  1:142422/135‑1 (MQ=255)
                         aaaGCAATCGAAAACATTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATgatgat                                                                                                            >  1:176303/1‑135 (MQ=255)
                         aaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATgatgat                                                                                                            >  2:261316/1‑135 (MQ=255)
                            gAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATtat                                                                                                         <  2:666775/135‑1 (MQ=255)
                               aTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACAAATGATGATTATAAt                                                                                                      >  1:428916/1‑135 (MQ=255)
                               aTCGAAAACTTTAAGTTTTTTGCGAAATATAGCAATTTTAAAGATTTAATCACTATTAAAGATGGAGTTATATCATATAATCAAGAAGTGCAGAGTTATTCGGTTAAATATCAATTAACTAATGATTATTATAAt                                                                                                      >  2:381685/1‑135 (MQ=18)
                                   aaaaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTATCTAATGATGATTATAATGTaa                                                                                                  <  1:4754/135‑1 (MQ=255)
                                   aaaaCCTTAAGTTTTTTGCGCAATATAGAAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTaa                                                                                                  <  2:490620/135‑1 (MQ=255)
                                    aaaCTTTAAGTTTTTTTCGAAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTaaa                                                                                                 <  2:411367/135‑1 (MQ=255)
                                    aaaCTTTAAGTTTTTTGCGAAATATAGCAATTTTAAAGATTTAATGAACTATAACGATGGAGATCTATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTACTGATGATTATAATGTaaa                                                                                                 <  2:186194/135‑1 (MQ=18)
                                       cTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACa                                                                                              <  1:305248/135‑1 (MQ=255)
                                       cTTTAAGTTTTTTGCGCAATAGAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACa                                                                                              <  2:128695/135‑1 (MQ=255)
                                          tAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTCTACAGAAGAAGATATATCATATAATAAAGGGTTTACTATGCAGTCAGCACAATACCAATTAACTAATGATGATTATAATGTAAAACAAtt                                                                                           >  2:636710/1‑135 (MQ=21)
                                          tAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAAtt                                                                                           >  1:82918/1‑135 (MQ=255)
                                          tAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACACtt                                                                                           >  2:367592/1‑135 (MQ=255)
                                                ttttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTCTAATGTAAAACAATTACgcaa                                                                                     <  2:359092/135‑4 (MQ=17)
                                                    gcgcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGa                                                                                 <  2:583881/135‑1 (MQ=255)
                                                      gcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGGGTTATTCGGCTCAATATCAATTAACTAATGCTGCTTATAATGTACACCCATTACGTTAAAGatc                                                                               >  1:599848/1‑134 (MQ=18)
                                                      gcAATATAGCAATTTTAAAGATTTAATGAATGATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGata                                                                               <  2:74240/135‑1 (MQ=14)
                                                         atatAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGa                                                                            <  2:215517/135‑1 (MQ=255)
                                                            taGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAGTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATAGGatat                                                                         >  1:21636/1‑135 (MQ=14)
                                                             aGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATa                                                                        <  1:63678/135‑1 (MQ=255)
                                                             aGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGAAATATGATATATTGCGGAGTAGCCGAGTTATTCAGCTCAATATCAATTAAATAATGATAATTATAATGTAAAACAATTACGTAAAAGATATGACATa                                                                        <  1:486721/135‑1 (MQ=18)
                                                                aaTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATAccg                                                                     <  2:499032/135‑2 (MQ=11)
                                                                aaTTTTAAAGATTTAAGGGATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTGATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGTTATTACATACca                                                                     <  2:33337/135‑1 (MQ=18)
                                                                  ttttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaac                                                                   <  2:638709/135‑1 (MQ=255)
                                                                  ttttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAACTCAATATCAATTAACTAATGATGATTATAATGTAAAACCATAACGTAAAAGATATGACATACcaac                                                                   >  2:4404/1‑135 (MQ=255)
                                                                  ttttAAAGATTTAATGAATTATAAAGATGGAGATATATCAGATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaac                                                                   >  2:590935/1‑135 (MQ=255)
                                                                   tttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGATATACcgacg                                                                  <  2:466169/135‑2 (MQ=1)
                                                                   tttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAATATATGACATACcaaca                                                                  <  2:310371/135‑1 (MQ=255)
                                                                   tttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaaaa                                                                  >  2:606335/1‑133 (MQ=255)
                                                                   tttAAAGATTTAATGAATTATAAAGATGAAAATATATCATATAATCCAGAGGTACCGAGTTATTCGGCTAAATATCAATTAACTAATTATGAATATAATTAAAAACAATTACGCAAAAGATATGATATACcaacg                                                                  >  1:249710/1‑134 (MQ=11)
                                                                         gATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATaaa                                                            <  1:222322/135‑1 (MQ=255)
                                                                           tttAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGCTTATAATGTAAAACAATTACGTAAAAGATATGACATACCACAAAATAAAgc                                                          >  2:318397/1‑135 (MQ=255)
                                                                           tttAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAGACAATTACGTAAAAGATATGACATACCAGCAAATAAAgc                                                          >  2:445241/1‑135 (MQ=255)
                                                                            ttAATGAATTATAAAGATGGAGATACATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgcg                                                         >  2:350320/1‑135 (MQ=255)
                                                                                tGAATTATAAAGATTGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAATCCATTACTTAAAAGATATGACATACCAACAAATAAAGCGCCGa                                                     <  1:620202/135‑1 (MQ=255)
                                                                                tGAATTATAAAAATGGAGATATATCATAGAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGa                                                     <  2:655770/135‑1 (MQ=255)
                                                                                  aactataAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACTAACAAATAAAGCGCCGAAg                                                   <  1:161825/132‑1 (MQ=255)
                                                                                   aTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGc                                                  <  2:399482/135‑1 (MQ=255)
                                                                                   aTTATAAAGATGGAGATACAGCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGc                                                  <  2:681097/135‑1 (MQ=255)
                                                                                    ttATAAATATGGAGATATATCATATAATCCAGAAGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATTACATACCAACTAATAAAGCGCCGAAGCt                                                 <  1:583866/135‑1 (MQ=255)
                                                                                        aaaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATATCAACAAATAAAGCGCCGAAGCtgttg                                             <  1:136827/135‑1 (MQ=255)
                                                                                        aaaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgttg                                             <  1:566331/135‑1 (MQ=255)
                                                                                        aaaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgttg                                             <  1:499526/135‑1 (MQ=255)
                                                                                        aaaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgttg                                             >  2:531681/1‑135 (MQ=255)
                                                                                        aaaGATGGAGATATATCATATAATCCAGAGGTACCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgttg                                             >  1:658260/1‑135 (MQ=255)
                                                                                                agaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAgg                                     >  2:425363/1‑135 (MQ=255)
                                                                                                       tCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAggg                              >  2:71148/1‑135 (MQ=255)
                                                                                                       tCATATAATCCAGAGGTGCCGAGATATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGATACAggg                              >  1:540768/1‑135 (MQ=255)
                                                                                                           ataATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAAtt                          >  2:476090/1‑135 (MQ=255)
                                                                                                           ataATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAAtt                          >  2:569625/1‑135 (MQ=255)
                                                                                                            taATCCAGAGGTGCCGAGTTATTCAGCTAAATATAAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATCACATATCCACAAATAAAGCGCCGAAGCGGTTGTTGAAAGGTAGAGGGAAttt                         >  1:31314/1‑135 (MQ=14)
                                                                                                                  agagGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATAAAGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAgg                   >  1:334505/1‑135 (MQ=255)
                                                                                                                     ggTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTAAAGGGAATTTAAAAGGTtc                >  2:279234/1‑135 (MQ=255)
                                                                                                                       tGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtcat              <  2:176303/135‑1 (MQ=255)
                                                                                                                       tGCCGAGTTATTCAGCTCAATATCAATTAACTAATACTGATTATAATGTAAAAAAATTAAGTAAAAGATATAAAATAAAGACAAAGCAAGAGCCGAACCTGTTGTTGAAAAGTAACGCGCATATAAAAGGTtcat              >  2:277599/1‑135 (MQ=255)
                                                                                                                          cGAGTTATTCAGCTCAATATCAATAAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAg           >  1:12392/1‑135 (MQ=255)
                                                                                                                             gTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGATTCATCAGTTg        >  2:652855/1‑135 (MQ=255)
                                                                                                                              ttATTCAGCTCAATATCAAATAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTgg       >  1:489850/1‑135 (MQ=255)
                                                                                                                                   cAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGAtata  >  1:220917/1‑135 (MQ=255)
                                                                                                                                   cAGCTCAATATCAATTAACTAATGATAAAAAAAAAGTAAAACAATAACGAAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGATGGAtata  >  1:222156/1‑135 (MQ=14)
                                                                                                                                    |                                                                                                                                     
GGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATA  >  USA300TCH1516_ALE/111380‑111645

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