BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,440,007	Δ1 bp	84.7%	coding (34/201 nt)	cspA_2 ←	Cold shock protein CspA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,440,007	0	C	.	84.7%	757.9 / 143.2	227	coding (34/201 nt)	cspA_2	Cold shock protein CspA
Reads supporting (aligned to +/- strand): ref base C (15/14); new base . (101/91); total (118/109)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

Only 100 of 227 total aligned reads displayed.

ATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAA  >  USA300TCH1516_ALE/1439920‑1440096

aTGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGc                                                                                        >  1:1289776/1‑91 (MQ=255)

   aTTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAACCCGAATCCTTTTT‑AGCGTTa                                                                                    >  2:520180/1‑91 (MQ=255)

    tttGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCTGAATCCTTTTTCAGCGTTa                                                                                    >  2:1459925/1‑91 (MQ=255)

    tttGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTaa                                                                                   >  1:1780065/1‑91 (MQ=255)

     ttGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTaa                                                                                   <  2:3471510/91‑1 (MQ=255)

     ttGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTaaa                                                                                  >  2:2419739/1‑91 (MQ=255)

      tGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAAc                                                                                 <  1:3004945/91‑1 (MQ=255)

      tGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAAc                                                                                 >  2:89647/1‑91 (MQ=255)

       gTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAAcc                                                                                >  2:2996666/1‑91 (MQ=255)

        tAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCa                                                                               >  1:1572938/1‑91 (MQ=255)

        tAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCa                                                                               >  2:2941424/1‑91 (MQ=255)

         aaCCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCAt                                                                              >  2:1732245/1‑91 (MQ=255)

           ccATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCAttt                                                                            >  2:1538080/1‑91 (MQ=255)

             aagatGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTa                                                                           >  1:917673/5‑91 (MQ=255)

             aTCTTGGTTAATTGCTGAAAAATGTACGAATCCGTCATTTTCTCCCTCCACTTCGATAATGCCGAGTCCTTTTT‑AGCGTCAAACCAtttca                                                                          >  1:3588267/1‑89 (MQ=255)

              tCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAAc                                                                         <  1:1033178/91‑1 (MQ=255)

                ttGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACt                                                                        <  2:70481/91‑1 (MQ=255)

                 tGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGt                                                                      >  2:167262/1‑91 (MQ=255)

                 tGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGt                                                                      <  1:2655693/91‑1 (MQ=255)

                 tGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGt                                                                      <  2:1442182/91‑1 (MQ=255)

                  ggTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGt                                                                      >  2:2191847/1‑91 (MQ=255)

                  ggTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑TTCTTTTAACCATTTAACTGtc                                                                     >  2:1272790/1‑90 (MQ=255)

                  ggTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTa                                                                     <  2:3586440/91‑1 (MQ=255)

                   gTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTAc                                                                    <  2:1595284/91‑1 (MQ=255)

                   gTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTAc                                                                    <  2:908698/91‑1 (MQ=255)

                      aaTTGCTGAAACATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAACGCCGACTCCTTCTT‑AGCGTTAAACCATTTAACTGTACCtt                                                                 >  2:986686/1‑91 (MQ=255)

                      aaTTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCt                                                                  >  1:3879313/1‑91 (MQ=255)

                       aTTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTg                                                                <  2:622064/91‑1 (MQ=255)

                       aTTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTg                                                                <  1:1168282/91‑1 (MQ=255)

                        ttGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGt                                                               >  2:3930775/1‑91 (MQ=255)

                        ttGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGt                                                               >  1:2580613/1‑91 (MQ=255)

                         tGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGtt                                                              >  2:2795850/1‑91 (MQ=255)

                         tGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGtt                                                              <  2:941573/91‑1 (MQ=255)

                          tcTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGttt                                                             <  2:2257991/90‑1 (MQ=255)

                          gCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGttt                                                             >  1:2897106/1‑91 (MQ=255)

                          gCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGttt                                                             >  1:3937846/1‑91 (MQ=255)

                           cTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTc                                                            >  2:3932853/1‑91 (MQ=255)

                           cTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTc                                                            <  2:606182/91‑1 (MQ=255)

                           cTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTc                                                            <  2:744753/91‑1 (MQ=255)

                           cTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTc                                                            <  1:2237779/91‑1 (MQ=255)

                            tGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCa                                                           <  1:1975792/91‑1 (MQ=255)

                            tGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCa                                                           >  1:1425287/1‑91 (MQ=255)

                            tGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCa                                                           <  2:369447/91‑1 (MQ=255)

                               aaaaTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATaa                                                        >  1:2447331/1‑91 (MQ=255)

                               aaaaTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATaa                                                        <  1:1903492/91‑1 (MQ=255)

                                aaaTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAAt                                                       >  1:3771460/1‑91 (MQ=255)

                                 aaTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATc                                                      <  1:152423/91‑1 (MQ=255)

                                 aaTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATc                                                      <  1:3223421/91‑1 (MQ=255)

                                 aaTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATc                                                      <  2:1034938/91‑1 (MQ=255)

                                 aaTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATc                                                      <  2:956796/91‑1 (MQ=255)

                                  aTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCt                                                     <  1:2217154/91‑1 (MQ=255)

                                   tGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATACTCTg                                                    >  1:3293379/1‑91 (MQ=255)

                                   tGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTg                                                    >  2:2556334/1‑91 (MQ=255)

                                    gTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGa                                                   >  2:2176306/1‑91 (MQ=255)

                                     tACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGaa                                                  >  2:3181315/1‑91 (MQ=255)

                                      aCGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGaaa                                                 >  1:2973328/1‑91 (MQ=255)

                                          aTACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTc                                             >  2:1332564/1‑91 (MQ=255)

                                            aCGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCa                                           >  1:2316240/1‑91 (MQ=255)

                                              gTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAg                                         <  1:2842821/91‑1 (MQ=255)

                                                cATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGAc                                       <  1:911727/91‑1 (MQ=255)

                                                cATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGAc                                       <  1:55864/91‑1 (MQ=255)

                                                  ttttCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTa                                     >  1:3119168/1‑91 (MQ=255)

                                                  ttttCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTa                                     >  2:76393/1‑91 (MQ=255)

                                                   tttCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTa                                     <  1:75550/91‑1 (MQ=255)

                                                   tttCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTaa                                    <  2:3074388/91‑1 (MQ=255)

                                                   tttCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTaa                                    <  1:3761007/91‑1 (MQ=255)

                                                      ctcCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAat                                 <  1:3830422/91‑1 (MQ=255)

                                                       tcCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAatt                                >  2:150281/1‑91 (MQ=255)

                                                       tcCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAatt                                >  2:2041597/1‑91 (MQ=255)

                                                       tcCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAatt                                >  1:521447/1‑91 (MQ=255)

                                                        ccTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAattc                               >  1:331275/1‑91 (MQ=255)

                                                          ttCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAATACTAAAattcat                             <  1:3048041/91‑1 (MQ=255)

                                                          ttCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAattcat                             <  1:1227428/91‑1 (MQ=255)

                                                           tCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAattcat                             <  1:694598/91‑1 (MQ=255)

                                                                ttCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAata                       <  1:1137916/91‑1 (MQ=255)

                                                                ttCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAata                       <  1:2813992/91‑1 (MQ=255)

                                                                  cGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATg                     <  2:2778382/91‑1 (MQ=255)

                                                                  cGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATg                     >  1:626677/1‑91 (MQ=255)

                                                                  cGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATg                     >  1:1494087/1‑91 (MQ=255)

                                                                   gATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAAGTGTACCTTGTTTCTTAATCTGAAACCTACAAGACTAAAATTCATTCAGTATGc                    >  1:3684150/1‑91 (MQ=255)

                                                                   gATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATATGAAACCCCCAAGTCTAAAATTCATTCAATATGc                    >  2:3420931/1‑91 (MQ=255)

                                                                    aTAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCt                   <  2:1843090/91‑1 (MQ=255)

                                                                     tAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCtt                  <  1:2052700/91‑1 (MQ=255)

                                                                     tAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCtt                  >  1:3518338/1‑91 (MQ=255)

                                                                      aaaGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTa                 <  2:3538944/91‑1 (MQ=255)

                                                                       aaGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTAt                <  2:1865221/91‑1 (MQ=255)

                                                                       aaGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTAt                <  1:2046878/91‑1 (MQ=255)

                                                                       aaGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTCACATAATCTGAAACCTCCAAGACTACAATTCATTCAACATGCTTAt                >  2:3229146/1‑91 (MQ=255)

                                                                         gCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTc              <  2:1104048/91‑1 (MQ=255)

                                                                         gCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTc              >  1:3359710/1‑91 (MQ=255)

                                                                           cGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGc            >  1:1856587/1‑91 (MQ=255)

                                                                           cGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGc            >  2:2562935/1‑91 (MQ=255)

                                                                                ccTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATAtt       >  2:1189991/1‑91 (MQ=255)

                                                                                 ccttttTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATAtt       <  1:2540233/89‑1 (MQ=255)

                                                                                 cTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATAtt       <  1:2532941/91‑1 (MQ=255)

                                                                                  cttttTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAACATTCATTCAATATGCTTATTCGCATATTa      >  1:3177008/2‑91 (MQ=255)

                                                                                   ttttCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTAc     <  2:2639450/91‑1 (MQ=255)

                                                                                     ttCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTTCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACaa   >  2:1796015/1‑91 (MQ=255)

                                                                                     ttCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACaa   >  2:2876197/1‑91 (MQ=255)

                                                                                      tCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACaaa  >  1:3104442/1‑91 (MQ=255)

ATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAA  >  USA300TCH1516_ALE/1439920‑1440096

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