BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,653,040	A→C	5.6%	L86V (TTA→GTA)	gcvT ←	Aminomethyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,653,040	0	A	C	5.6%	622.8 / 20.0	189	L86V (TTA→GTA)	gcvT	Aminomethyltransferase
Reads supporting (aligned to +/- strand): ref base A (79/99); new base C (2/9); total (81/108)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 192 total aligned reads displayed.

TAACAACTAATAAATAATTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTAGCATCTTTACCTGTTACTTCAATTTCACC  >  USA300TCH1516_ALE/1652950‑1653130

tAACAACTAATAAATAATTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATaa                                                                                            <  2:222579/91‑1 (MQ=255)

 aacaacTAATAAATAATTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAAt                                                                                           <  2:1860755/91‑1 (MQ=255)

  acaacTAATAAATAATTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATg                                                                                          <  2:1125005/91‑1 (MQ=255)

   caacTAATAAATAATTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGc                                                                                         >  1:16498/1‑91 (MQ=255)

   caacTAATAAATAATTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGc                                                                                         <  1:1742092/91‑1 (MQ=255)

    aacTAATAAATATTTTTACTCAGCTGATTTATATGTTACTATATCATCGATGATACCGTCTTCTTCATTACATAAAGCAGTGTATACTGCt                                                                                        <  2:2087702/91‑1 (MQ=255)

    aacTAATAAATAATTTTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCATTATATACTGCt                                                                                        <  2:293177/91‑1 (MQ=255)

    aacTAATAAATAATTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATACAGCAGTATATAATGCt                                                                                        >  2:3714388/1‑91 (MQ=255)

    aacTAATAAATAATTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCATTATATACTGCt                                                                                        <  2:561190/91‑1 (MQ=255)

    aacTAATAAATAATTGTCGTCAGCTAATTTATATATCACTAAATCATCAATAATACCGCCTTCTTCAGTACATAAAGCAGTATATAATGcc                                                                                        <  2:3508770/91‑2 (MQ=255)

     acTAATAAATAATTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCtt                                                                                       <  1:357742/91‑1 (MQ=255)

      cTAATAAATAATTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCttt                                                                                      >  2:940572/1‑91 (MQ=255)

      cTAATAAATAATTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCttt                                                                                      <  1:1497311/91‑1 (MQ=255)

      cTAATAAATAATTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCATCTTCATTACATAAAGCAGTATATAATGCttt                                                                                      <  2:3816499/91‑1 (MQ=255)

       taataaATAATTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCATTATATACTGCttta                                                                                     <  2:1922854/91‑2 (MQ=255)

       taataaATAATTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCATTATATACTGCttta                                                                                     <  2:3078026/91‑2 (MQ=255)

       taataaATAATTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCtttt                                                                                     >  2:668541/1‑91 (MQ=255)

       taataaATAATTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCtttt                                                                                     <  2:1291163/91‑1 (MQ=255)

        aataaataATTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCATTATATACTGCtttat                                                                                    <  2:2082427/91‑3 (MQ=255)

         ataaataaTTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTGCATAAAGCAGTATATAATGCTTTTGa                                                                                   <  1:1020465/91‑1 (MQ=255)

         ataaataaTTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGa                                                                                   >  1:1781389/1‑91 (MQ=255)

          taaataaTTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCATTATAGACTGCtttatgt                                                                                  <  1:2987680/91‑5 (MQ=255)

          taaataaTTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGaa                                                                                  <  1:2483953/91‑1 (MQ=255)

           aaataaTTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAag                                                                                 >  1:2168193/1‑91 (MQ=255)

             ataatTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAagta                                                                               >  1:1471248/1‑91 (MQ=255)

             ataatTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAagta                                                                               <  1:1564310/91‑1 (MQ=255)

             ataatTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAagta                                                                               >  2:3398411/1‑91 (MQ=255)

               aatTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAagtagt                                                                             <  1:423415/91‑1 (MQ=255)

                atTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGtt                                                                            <  1:2401879/91‑1 (MQ=255)

                   gtcgtcAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTaaa                                                                         <  1:3035763/91‑1 (MQ=255)

                   gtcgtcAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTaaa                                                                         >  1:2586268/1‑91 (MQ=255)

                       tcAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTAt                                                                     >  2:3048251/1‑91 (MQ=255)

                        cAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATc                                                                    <  1:1972124/91‑1 (MQ=255)

                        cAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATc                                                                    <  2:2220048/91‑1 (MQ=255)

                        cAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATc                                                                    <  1:529926/91‑1 (MQ=255)

                        cAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATc                                                                    <  2:1232782/91‑1 (MQ=255)

                         aGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCa                                                                   <  2:444913/91‑1 (MQ=255)

                          gCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAg                                                                  <  2:154018/91‑1 (MQ=255)

                          gCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAg                                                                  <  1:650844/91‑1 (MQ=255)

                             aaTTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTAt                                                               <  2:586417/91‑1 (MQ=255)

                               tttATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCa                                                             <  2:572397/91‑1 (MQ=255)

                               tttATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCa                                                             <  1:1251162/91‑1 (MQ=255)

                               tttATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCa                                                             <  1:18105/91‑1 (MQ=255)

                               tttATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGCTAAATTATCAGTATCa                                                             >  2:2050475/1‑91 (MQ=255)

                                ttATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCAt                                                            >  1:1388429/1‑91 (MQ=255)

                                ttATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCAt                                                            >  1:912056/1‑91 (MQ=255)

                                    atatTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGa                                                        >  1:691395/1‑91 (MQ=255)

                                     tatTACTAAAACATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGat                                                       <  2:16514/91‑1 (MQ=255)

                                      atTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGata                                                      <  2:1783469/91‑1 (MQ=255)

                                       ttACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGataa                                                     <  2:3250259/91‑1 (MQ=255)

                                       ttACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGataa                                                     >  1:254011/1‑91 (MQ=255)

                                       ttACCAAATCATCAATCATACCGCCCTCTTCCCTCCATCCAGCACTCTATAATACTTCTAAAGTAGTTAAATTATCACTATCCCTTGacca                                                     >  2:2503437/1‑88 (MQ=255)

                                         aCTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGataata                                                   <  1:917965/91‑1 (MQ=255)

                                          cTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGAtaataa                                                  <  2:1158955/91‑1 (MQ=255)

                                          cTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGAtaataa                                                  <  2:1636216/91‑1 (MQ=255)

                                           tAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATaataaa                                                 >  2:2161897/1‑91 (MQ=255)

                                           tAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATaataaa                                                 >  2:2738216/1‑91 (MQ=255)

                                                catcaATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTg                                            >  2:2629138/1‑91 (MQ=255)

                                                 atcaATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGc                                           >  1:3108713/1‑91 (MQ=255)

                                                  tcaATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGca                                          <  1:2417196/91‑1 (MQ=255)

                                                   caATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGcac                                         >  2:1298476/1‑91 (MQ=255)

                                                   caATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGcac                                         >  1:2715855/1‑91 (MQ=255)

                                                   caATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGcac                                         >  2:1298916/1‑91 (MQ=255)

                                                     ataataCCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACaa                                       <  2:1351642/91‑1 (MQ=255)

                                                     ataataCCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACaa                                       <  2:673153/91‑1 (MQ=255)

                                                     ataataCCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACaa                                       <  2:746238/91‑1 (MQ=255)

                                                          aCCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGa                                  >  1:1193185/1‑91 (MQ=255)

                                                           ccGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGAc                                 <  2:291789/91‑1 (MQ=255)

                                                           ccGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGAc                                 <  2:852613/91‑1 (MQ=255)

                                                           ccGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATCGCACAAACTGAc                                 <  1:934123/91‑1 (MQ=255)

                                                            cGCCTTCTTCATTACATAAAGCAGTAGATAATGCTTTTGAAGTAGGTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACt                                <  2:1418632/91‑1 (MQ=255)

                                                             gCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTa                               >  1:2906428/1‑91 (MQ=255)

                                                             gCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTa                               >  1:2982317/1‑91 (MQ=255)

                                                             gCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTa                               >  2:3787877/1‑91 (MQ=255)

                                                              ccTTCTTCCTTACACAAAGCAGTATATAATGCTTTAGAAGTAGTCAAATTAACAGTATCATATGCTAATAAATATCGCACACACTGACTAg                              <  2:2791260/91‑1 (MQ=255)

                                                              ccTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTAg                              <  2:1228374/91‑1 (MQ=255)

                                                              ccTTCTTCAGTACATAAAGCAGTAGATAATGCTGTTGAAGGAGTTAAAGTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTAg                              <  1:1817270/91‑1 (MQ=255)

                                                                ttcttcATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTAGCa                            <  2:1261117/91‑1 (MQ=255)

                                                                ttcttcATTACATAAAGCAGTATATAATGCTCTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTAGCa                            <  2:1731584/91‑1 (MQ=255)

                                                                 tcttcATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTAGCAt                           <  1:3183880/91‑1 (MQ=255)

                                                                  cttcATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTCTCATTTGATAATAAATATTGCACAAACTGACTAGCATc                          >  2:899510/1‑91 (MQ=255)

                                                                  cttcATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTAGCATc                          >  2:2868659/1‑91 (MQ=255)

                                                                    tcATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTAGCATCtt                        <  1:252342/91‑1 (MQ=255)

                                                                      aTTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTAGCATCTTTa                      >  2:2805316/1‑91 (MQ=255)

                                                                        tACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTAGCATCTTTAcc                    <  2:2277550/91‑1 (MQ=255)

                                                                            tAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTAGCATCTTTACCTGtt                >  2:2065064/1‑91 (MQ=255)

                                                                             aaaGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTAGCATCTTTACCTGTTa               <  1:3428842/91‑1 (MQ=255)

                                                                              aaGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTAGCATCTTTACCTGTTAc              <  1:1046264/91‑1 (MQ=255)

                                                                               aGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTAGCATCTTTACCTGTTACt             >  2:615517/1‑91 (MQ=255)

                                                                               aGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTAGCATCTTTACCTGTTACt             <  1:494302/91‑1 (MQ=255)

                                                                               aGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTAGCATCTTTACCTGTTACt             <  2:2131675/91‑1 (MQ=255)

                                                                                gCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTAGCATCTTTACCTGTTACtt            >  1:3751772/1‑91 (MQ=255)

                                                                                gCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTAGCATCTTTACCTGTTACtt            <  1:405390/91‑1 (MQ=255)

                                                                                  aGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTAGCATCTTTACCTGTTACTTCa          <  1:3393322/91‑1 (MQ=255)

                                                                                    tatataATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTAGCATCTTTACCTGTTACTTCAAt        >  1:1717251/1‑91 (MQ=255)

                                                                                      tataATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTAGCATCTTTACCTGTTACTTCAAttt      >  1:1888557/1‑91 (MQ=255)

                                                                                      tataATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTAGCATCTTTACCTGTTACTTCAAttt      >  2:1835401/1‑91 (MQ=255)

                                                                                       ataATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTAGCATCTTTACCTGTTACTTCAATTTc     >  2:2382585/1‑91 (MQ=255)

                                                                                       ataATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTAGCACCTTTACCTGTTACTTCGATTTc     >  2:3134727/1‑91 (MQ=255)

                                                                                          aTGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATACTGCACAAACTGACTAGCATCTTTACCTGTTACTTCAATTTCAcc  >  1:3080216/1‑91 (MQ=255)

TAACAACTAATAAATAATTGTCGTCAGCTAATTTATATATTACTAAATCATCAATAATACCGCCTTCTTCATTACATAAAGCAGTATATAATGCTTTTGAAGTAGTTAAATTATCAGTATCATTTGATAATAAATATTGCACAAACTGACTAGCATCTTTACCTGTTACTTCAATTTCACC  >  USA300TCH1516_ALE/1652950‑1653130

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