BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	806,422	T→G	7.0%	I127I (ATA→ATC)	ribN ←	Riboflavin transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	806,422	0	T	G	7.0%	392.9 / 16.2	114	I127I (ATA→ATC)	ribN	Riboflavin transporter
Reads supporting (aligned to +/- strand): ref base T (56/50); new base G (4/4); total (60/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 117 total aligned reads displayed.

AAAAATACCTGAGAATAAGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGAAAGGATTGAGTTTCATTAA  >  USA300TCH1516_ALE/806332‑806511

aaaaaTACCTGAGAATAAGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAAtt                                                                                           >  1:4108450/1‑91 (MQ=255)

    aTACCTGAGAATAAGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCa                                                                                       >  1:906536/1‑91 (MQ=255)

    aTACCTGAGAATAAGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCa                                                                                       >  1:906403/1‑91 (MQ=255)

      aCCTGAGAATAAGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCAtt                                                                                     >  1:1513668/1‑91 (MQ=255)

      aCCTGAGAATAAGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCAtt                                                                                     <  1:285057/91‑1 (MQ=255)

       ccTGAGAATAAGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTg                                                                                    >  1:2757069/1‑91 (MQ=255)

       ccTGAGAATAAGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTg                                                                                    <  1:2568456/91‑1 (MQ=255)

       ccTGAGAATAAGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTg                                                                                    >  1:586584/1‑91 (MQ=255)

       ccTGAGAATAAGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTg                                                                                    <  1:1268686/91‑1 (MQ=255)

       ccTGAGAATAAGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTg                                                                                    <  1:1268948/91‑1 (MQ=255)

       ccTGAGAATAAGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATGATGATTg                                                                                    <  1:3276868/91‑1 (MQ=255)

       ccTGAGAATAAGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATGATAATTg                                                                                    <  1:53608/91‑1 (MQ=255)

       ccTGAGAATAAGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATGATAATTg                                                                                    <  1:53435/91‑1 (MQ=255)

       ccTGAGAATAAGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATGATAATTg                                                                                    <  1:1551670/91‑1 (MQ=255)

        cTGAGAATAAGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGc                                                                                   >  1:521699/1‑91 (MQ=255)

        cTGAGAATAAGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGc                                                                                   >  1:68921/1‑91 (MQ=255)

        cTGAGAATAAGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGc                                                                                   <  1:2620363/91‑1 (MQ=255)

        cTGAGAATAAGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGc                                                                                   >  1:3939201/1‑91 (MQ=255)

           agaATAAGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGt                                                                                <  1:1548665/91‑1 (MQ=255)

           agaATAAGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGt                                                                                >  1:3614283/1‑91 (MQ=255)

           agaATAAGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGt                                                                                >  1:3615193/1‑91 (MQ=255)

               tAAGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAAtt                                                                            <  1:1836691/91‑1 (MQ=255)

               tAAGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAAtt                                                                            <  1:3373410/91‑1 (MQ=255)

                aaGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAAttt                                                                           >  1:2029965/1‑91 (MQ=255)

                 aGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTg                                                                          >  1:3067443/1‑91 (MQ=255)

                  gCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGa                                                                         >  1:1484282/1‑91 (MQ=255)

                  gCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGa                                                                         >  1:2281715/1‑91 (MQ=255)

                          atgatgGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCt                                                                 <  1:3364579/91‑1 (MQ=255)

                          atgatgGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCt                                                                 <  1:2995212/91‑1 (MQ=255)

                          atgatgGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCt                                                                 <  1:2932600/91‑1 (MQ=255)

                           tgatgGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCtt                                                                >  1:2497059/1‑91 (MQ=255)

                              tgGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTAcc                                                             >  1:8809/1‑91 (MQ=255)

                                gAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCtt                                                           <  1:1954003/91‑1 (MQ=255)

                                gAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCtt                                                           <  1:3266630/91‑1 (MQ=255)

                                 aaTCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCttt                                                          >  1:2539776/1‑91 (MQ=255)

                                   tCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTc                                                        <  1:1710556/91‑1 (MQ=255)

                                   tCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTc                                                        >  1:3963288/1‑91 (MQ=255)

                                   tCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTc                                                        <  1:2287896/91‑1 (MQ=255)

                                     aTAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCAt                                                      <  1:2841498/91‑1 (MQ=255)

                                     aTAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCAt                                                      <  1:4090193/91‑1 (MQ=255)

                                     aTAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCAt                                                      <  1:1074571/91‑1 (MQ=255)

                                      tAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATg                                                     >  1:1587501/1‑91 (MQ=255)

                                        gatgatGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTa                                                   >  1:2129162/1‑91 (MQ=255)

                                         atgatgAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTaa                                                  >  1:2852353/1‑91 (MQ=255)

                                         atgatgAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTaa                                                  >  1:2848744/1‑91 (MQ=255)

                                          tgatgaAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAg                                                 >  1:679679/1‑91 (MQ=255)

                                           gatgaAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGa                                                >  1:330386/1‑91 (MQ=255)

                                             tgaAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGaaa                                              <  1:939171/91‑1 (MQ=255)

                                               aaaaTTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAAc                                            <  1:742395/91‑1 (MQ=255)

                                               aaaaTTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAAc                                            <  1:740095/91‑1 (MQ=255)

                                                aaaTTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACa                                           <  1:2848259/91‑1 (MQ=255)

                                                aaaTTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACa                                           >  1:409847/1‑91 (MQ=255)

                                                 aaTTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACaa                                          <  1:4039169/91‑1 (MQ=255)

                                                 aaTTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACaa                                          <  1:1446328/91‑1 (MQ=255)

                                                   ttCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATg                                        <  1:2181463/91‑1 (MQ=255)

                                                   ttCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATg                                        <  1:2181517/91‑1 (MQ=255)

                                                    tCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGc                                       <  1:1407938/91‑1 (MQ=255)

                                                     cAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCt                                      >  1:1962343/1‑91 (MQ=255)

                                                      aGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCtt                                     <  1:457249/91‑1 (MQ=255)

                                                        gTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTtaa                                   >  1:456207/1‑91 (MQ=255)

                                                         tttCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTtaat                                  <  1:56718/91‑1 (MQ=255)

                                                            cACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATaaa                               <  1:1681778/91‑1 (MQ=255)

                                                            cACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATaaa                               >  1:2026215/1‑91 (MQ=255)

                                                             aCTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATaaaa                              >  1:1194953/1‑91 (MQ=255)

                                                                 ttAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGt                          >  1:2186678/1‑91 (MQ=255)

                                                                 ttAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGt                          >  1:2523678/1‑91 (MQ=255)

                                                                   aaaaGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCATACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTcc                        >  1:401014/1‑91 (MQ=255)

                                                                   aaaaGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCATACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTcc                        >  1:401028/1‑91 (MQ=255)

                                                                    aaaGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCa                       >  1:789636/1‑91 (MQ=255)

                                                                     aaGCATACCTAAGATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAg                      >  1:3190783/1‑91 (MQ=255)

                                                                     aaGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAg                      >  1:4052889/1‑91 (MQ=255)

                                                                      aGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGa                     >  1:1359099/1‑91 (MQ=255)

                                                                      aGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGa                     <  1:4087763/91‑1 (MQ=255)

                                                                        cATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGaaa                   >  1:2806660/1‑91 (MQ=255)

                                                                        cATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGaaa                   >  1:934813/1‑91 (MQ=255)

                                                                         aTACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGAAAg                  >  1:3112029/1‑91 (MQ=255)

                                                                          tACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGAAAgg                 <  1:1477462/91‑1 (MQ=255)

                                                                          tACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGAAAgg                 <  1:1477701/91‑1 (MQ=255)

                                                                            ccTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGAAAGGAt               <  1:1312428/91‑1 (MQ=255)

                                                                             cTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGAAAGGAtt              <  1:2381406/91‑1 (MQ=255)

                                                                             cTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGAAAGGAtt              <  1:3498084/91‑1 (MQ=255)

                                                                              tAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGAAAGGATTg             >  1:1266197/1‑91 (MQ=255)

                                                                                 aaTAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCNTGTAAGAAAACNATGCNTAATAAAATCGTCCAGAAAGGATTGAGt          >  1:1135607/1‑91 (MQ=255)

                                                                                 aaTAACAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGAAAGGATTGAGt          <  1:3074524/91‑1 (MQ=255)

                                                                                  acaGCAATGATAATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGAAAGGATTGAGtt         >  1:182907/3‑91 (MQ=255)

                                                                                  aTAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTGTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGAAAGGATTGAGtt         <  1:2697575/91‑1 (MQ=255)

                                                                                   caGCAATGATAATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGAAAGGATTGAGttt        >  1:1969264/2‑91 (MQ=255)

                                                                                    aGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGAAAGGATTGAGTTTc       >  1:380531/1‑91 (MQ=255)

                                                                                    aGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGAAAGGATTGAGTTTc       >  1:3044358/1‑91 (MQ=255)

                                                                                    aGCAATGATAATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGAAAGGATTGAGTTTc       >  1:2440289/1‑91 (MQ=255)

                                                                                    aGCAATGATAATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGAAAGGATTGAGTTTc       >  1:2440290/1‑91 (MQ=255)

                                                                                     gCAATTATCATTGCTGTAATTTGATACTTCCTTGCCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGAAAGGATTGAGTTTCa      <  1:1248922/91‑1 (MQ=255)

                                                                                     gCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGAAAGGATTGAGTTTCa      >  1:500410/1‑91 (MQ=255)

                                                                                     gCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGAAAGGATTGAGTTTCa      <  1:2866417/91‑1 (MQ=255)

                                                                                     gCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGAAAGGATTGAGTTTCa      <  1:1248854/91‑1 (MQ=255)

                                                                                     gCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGAAAGGATTGAGTTTCa      <  1:2288946/91‑1 (MQ=255)

                                                                                      cAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGAAAGGATTGAGTTTCAt     >  1:3640172/1‑91 (MQ=255)

                                                                                      cAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGAAAGGATTGAGTTTCAt     >  1:3638812/1‑91 (MQ=255)

                                                                                        aTTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGAAAGGATTGAGTTTCATTa   >  1:1675223/1‑91 (MQ=255)

                                                                                         ttATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGAAAGGATTGAGTTTCATTaa  <  1:4020000/91‑1 (MQ=255)

AAAAATACCTGAGAATAAGCCAGCTAATGATGGAATCATAGATGATGAAAATTCAGGTTTCACTATTAAAAGCATACCTAAAATAGCAATTATCATTGCTGTAATTTGATACTTCCTTACCTTTTCATGTAAGAAAACAATGCTTAATAAAATCGTCCAGAAAGGATTGAGTTTCATTAA  >  USA300TCH1516_ALE/806332‑806511

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