BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,562,886	A→G	V299A (GTT→GCT)	b1489 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,562,886	0	A	G	98.2%	520.4 / ‑0.5	167	V299A (GTT→GCT)	b1489	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/3); new base G (83/81); total (83/84)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-01

Only 100 of 172 total aligned reads displayed.

GCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAACATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGA  >  NC_000913/1562796‑1562975

gCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCa                                                                                         <  2:3705040/93‑1 (MQ=255)

gCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCa                                                                                         <  2:798937/93‑1 (MQ=255)

gCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCa                                                                                         <  1:3484232/93‑1 (MQ=255)

  tGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATg                                                                                       >  2:728064/1‑93 (MQ=255)

   gACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGc                                                                                      >  1:1223980/1‑93 (MQ=255)

    aCCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCg                                                                                     >  2:2546910/1‑93 (MQ=255)

     ccAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGa                                                                                    <  1:1181904/93‑1 (MQ=255)

     ccAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGa                                                                                    >  1:188891/1‑93 (MQ=255)

      cAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGAt                                                                                   >  1:1996381/1‑93 (MQ=255)

      cAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGAt                                                                                   >  1:1697074/1‑93 (MQ=255)

        gCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTc                                                                                 <  1:229923/93‑1 (MQ=255)

          tttGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGt                                                                               >  2:3582372/1‑93 (MQ=255)

          tttGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGt                                                                               <  1:4544435/93‑1 (MQ=255)

          tttGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGt                                                                               <  1:1636608/93‑1 (MQ=255)

           ttGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGtt                                                                              <  1:2204885/93‑1 (MQ=255)

           ttGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGtt                                                                              >  1:4082619/1‑93 (MQ=255)

            tGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGAGCGAACAGCGAAGCATGCGATTCGTTg                                                                             >  2:223585/1‑93 (MQ=255)

             gcgcATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGa                                                                            <  2:4036491/93‑1 (MQ=255)

             gcgcATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGa                                                                            >  2:3991135/1‑93 (MQ=255)

              cgcATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAg                                                                           <  1:2073425/93‑1 (MQ=255)

               gcATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGt                                                                          >  1:3260614/1‑93 (MQ=255)

                 aTTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTAc                                                                        >  1:1173610/1‑93 (MQ=255)

                  ttttGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACa                                                                       >  2:2997154/1‑93 (MQ=255)

                  ttttGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACa                                                                       >  2:1383782/1‑93 (MQ=255)

                  ttttGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACa                                                                       <  1:1474962/93‑1 (MQ=255)

                    ttGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGa                                                                     <  2:2421704/93‑1 (MQ=255)

                     tGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAACATGCGATTCGTTGAGTACAGAt                                                                    <  1:4108963/93‑1 (MQ=255)

                      gAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGAtt                                                                   >  2:894581/1‑93 (MQ=255)

                       aaTTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTc                                                                  >  2:577599/1‑93 (MQ=255)

                       aaTTTCTGCACAGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTc                                                                  <  1:2968486/93‑1 (MQ=255)

                         tttCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGa                                                                <  1:2921876/93‑1 (MQ=255)

                           tCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATg                                                              <  1:3332705/93‑1 (MQ=255)

                            cTGCACCGTTGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGt                                                             <  1:1101748/93‑1 (MQ=255)

                             tGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTGCAGATTCGATGtt                                                            >  2:3757116/1‑93 (MQ=255)

                             tGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGtt                                                            >  1:3978112/1‑93 (MQ=255)

                             tGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGtt                                                            >  1:4386936/1‑93 (MQ=255)

                             tGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGtt                                                            >  2:4376560/1‑93 (MQ=255)

                              gCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTa                                                           >  2:3160992/1‑93 (MQ=255)

                               cACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTAc                                                          <  2:2838032/93‑1 (MQ=255)

                                aCCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACg                                                         >  2:1367557/1‑93 (MQ=255)

                                 ccGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGa                                                        >  1:535252/1‑93 (MQ=255)

                                    tGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACaa                                                     >  2:1185328/1‑93 (MQ=255)

                                    tGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACaa                                                     >  1:128115/1‑93 (MQ=255)

                                     gggAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACaaa                                                    <  2:1007570/93‑1 (MQ=255)

                                     gggAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACaaa                                                    >  2:3622063/1‑93 (MQ=255)

                                       gAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATg                                                  <  2:4005497/93‑1 (MQ=255)

                                        aaGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGa                                                 >  1:439080/1‑93 (MQ=255)

                                         aGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGAt                                                >  2:1566019/1‑93 (MQ=255)

                                         aGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGAt                                                <  1:2736878/93‑1 (MQ=255)

                                         aGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGAt                                                <  2:290941/93‑1 (MQ=255)

                                          gATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGAtt                                               <  2:3119388/93‑1 (MQ=255)

                                          gATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGAtt                                               >  1:2507149/1‑93 (MQ=255)

                                          gATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGAtt                                               >  1:2679109/1‑93 (MQ=255)

                                               cGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTcccc                                          <  2:1864248/93‑1 (MQ=255)

                                                gCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTccccc                                         <  1:4614882/93‑1 (MQ=255)

                                                 cccAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTGCGACAAATGATTTCCCCCa                                        >  1:2744911/1‑93 (MQ=255)

                                                  ccAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCAt                                       <  2:817116/93‑1 (MQ=255)

                                                   cAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCAtt                                      >  2:1726260/1‑93 (MQ=255)

                                                     gtgtATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTc                                    >  2:1031878/1‑93 (MQ=255)

                                                     gtgtATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTc                                    >  1:3633118/1‑93 (MQ=255)

                                                     gtgtATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAACATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTc                                    <  2:2762539/93‑1 (MQ=255)

                                                      tgtATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCa                                   <  1:2842922/93‑1 (MQ=255)

                                                          tCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGaa                               >  2:397984/1‑93 (MQ=255)

                                                           cGGCGTCCCGTTGCGAAGTGCGAACAGCGAAGCATGCGATTCGTTGATTACAGATTCGTTGTTACGACAATTGATTTCCCCCATTTCATGaaa                              <  2:4582198/93‑1 (MQ=255)

                                                           cGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGaaa                              <  1:270257/93‑1 (MQ=255)

                                                            ggCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAAt                             <  1:1620742/93‑1 (MQ=255)

                                                               aTCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGc                          >  1:2523836/1‑93 (MQ=255)

                                                               aTCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGc                          >  1:4413665/1‑93 (MQ=255)

                                                                tCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCg                         >  2:4139896/1‑93 (MQ=255)

                                                                tCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCg                         <  1:1863628/93‑1 (MQ=255)

                                                                tCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCg                         <  1:4022289/93‑1 (MQ=255)

                                                                tCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCg                         >  2:1215586/1‑93 (MQ=255)

                                                                 cccGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCgg                        >  1:863456/1‑93 (MQ=255)

                                                                 cccGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCgg                        >  2:2062173/1‑93 (MQ=255)

                                                                  ccGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCggg                       >  2:4340502/1‑93 (MQ=255)

                                                                   cGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGgc                      <  2:4551315/93‑1 (MQ=255)

                                                                   cGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGgc                      <  2:1976474/93‑1 (MQ=255)

                                                                   cGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGgc                      <  1:2105720/93‑1 (MQ=255)

                                                                    gTTGCGCAGTGCGACCAGAGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGgct                     <  1:705236/93‑1 (MQ=255)

                                                                    gTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGgct                     <  1:3703696/93‑1 (MQ=255)

                                                                     ttGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTACCCCATTTCATGAAATGGCGGgctg                    >  2:2783825/1‑93 (MQ=255)

                                                                      tGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGgctgc                   <  2:4664425/93‑1 (MQ=255)

                                                                         gcAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATCATTTCCCCCATTTCATGAAATGGCGGGCtgctgc                <  2:2976105/93‑1 (MQ=255)

                                                                          cAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGca               >  1:1012131/1‑93 (MQ=255)

                                                                            gTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGcaca             >  2:4195480/1‑93 (MQ=255)

                                                                             tGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACAt            <  1:3456386/93‑1 (MQ=255)

                                                                             ggCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACAt            <  1:1388067/92‑1 (MQ=255)

                                                                              gCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATg           <  2:1157414/93‑1 (MQ=255)

                                                                                gAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATggc         <  1:728064/93‑1 (MQ=255)

                                                                                gAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATggc         <  1:3200033/93‑1 (MQ=255)

                                                                                gAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATggc         <  1:1777428/93‑1 (MQ=255)

                                                                                gAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATggc         <  1:1330074/93‑1 (MQ=255)

                                                                                 aaCAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATggcg        >  2:2245280/1‑93 (MQ=255)

                                                                                 aaCAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATggcg        <  2:4299939/93‑1 (MQ=255)

                                                                                  aCAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATggggg       >  2:3048212/1‑93 (MQ=255)

                                                                                   cAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATggcggg      >  2:974697/1‑92 (MQ=255)

                                                                                     gCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCga    >  2:3391183/1‑93 (MQ=255)

                                                                                     gCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCga    >  2:1675406/1‑93 (MQ=255)

                                                                                     gCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCga    <  2:150351/93‑1 (MQ=255)

                                                                                     gCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCga    <  1:1035459/93‑1 (MQ=255)

GCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAACATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGA  >  NC_000913/1562796‑1562975

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