BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,804,546	T→C	E50G (GAG→GGG)	rfaG ←	glucosyltransferase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,804,546	0	T	C	99.6%	787.0 / ‑4.8	241	E50G (GAG→GGG)	rfaG	glucosyltransferase I
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (116/124); total (117/124)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.90e-01

Only 100 of 244 total aligned reads displayed.

GGATGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCTCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATA  >  NC_000913/3804456‑3804635

ggATGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCa                                                                                         >  1:1150353/1‑93 (MQ=255)

ggATGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCa                                                                                         <  1:4808461/93‑1 (MQ=255)

    gCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATg                                                                                     >  1:813141/1‑93 (MQ=255)

     ctctTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGc                                                                                    <  2:4823680/93‑1 (MQ=255)

      tctTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCt                                                                                   >  2:4929875/1‑93 (MQ=255)

      tctTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCt                                                                                   >  1:433482/1‑93 (MQ=255)

      tctTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCt                                                                                   >  1:5919185/1‑93 (MQ=255)

      tctTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCt                                                                                   >  1:4956204/1‑93 (MQ=255)

       ctTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCtt                                                                                  <  1:2470397/93‑1 (MQ=255)

        tttGATATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCttt                                                                                 <  2:825464/93‑1 (MQ=255)

        tttGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCttt                                                                                 <  2:3194004/93‑1 (MQ=255)

        tttGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCttt                                                                                 <  1:2790004/93‑1 (MQ=255)

         ttGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTc                                                                                <  2:5671511/93‑1 (MQ=255)

         ttGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTc                                                                                >  2:452840/1‑93 (MQ=255)

          tGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCg                                                                               <  1:6204840/93‑1 (MQ=255)

          tGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCg                                                                               >  2:4112693/1‑93 (MQ=255)

          tGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCg                                                                               <  2:2534786/93‑1 (MQ=255)

          tGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGCTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCg                                                                               <  1:4434732/93‑1 (MQ=255)

           gagaTGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCgg                                                                              >  2:5165581/1‑93 (MQ=255)

           gagaTGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCgg                                                                              <  1:3596791/93‑1 (MQ=255)

             gatgatTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCggga                                                                            >  2:3450440/1‑92 (MQ=255)

             gatgatTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGc                                                                            >  2:3769848/1‑93 (MQ=255)

                gatTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAAt                                                                         <  1:561914/93‑1 (MQ=255)

                gatTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAAt                                                                         <  2:4263345/93‑1 (MQ=255)

                gatTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAAt                                                                         <  2:3120945/93‑1 (MQ=255)

                 atTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATc                                                                        <  2:5377313/93‑1 (MQ=255)

                 atTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATc                                                                        <  2:2679702/93‑1 (MQ=255)

                         cccAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTccc                                                                >  1:5016941/1‑93 (MQ=255)

                          ccAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCa                                                               <  1:1802939/93‑1 (MQ=255)

                           cAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCAc                                                              <  1:404480/93‑1 (MQ=255)

                           cAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCAc                                                              <  1:4155717/93‑1 (MQ=255)

                            aGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACg                                                             >  1:6233261/1‑93 (MQ=255)

                            aGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACg                                                             >  2:3690444/1‑93 (MQ=255)

                             ggCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGa                                                            >  2:5110501/1‑93 (MQ=255)

                               cATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACt                                                          <  1:1060210/93‑1 (MQ=255)

                                  aataTTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACtgtg                                                       >  1:2669997/1‑93 (MQ=255)

                                  aataTTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACtgtg                                                       >  1:699837/1‑93 (MQ=255)

                                   atatTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTgtgt                                                      >  2:1066728/1‑93 (MQ=255)

                                    tatTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGta                                                     <  1:4452373/93‑1 (MQ=255)

                                     atTCTGCAGTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAAGAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGtat                                                    <  1:2081744/93‑1 (MQ=255)

                                      ttCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGtata                                                   <  1:1243222/93‑1 (MQ=255)

                                       tCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGtatat                                                  <  2:4325938/93‑1 (MQ=255)

                                       tCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGtatat                                                  <  1:2017900/93‑1 (MQ=255)

                                       tCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGtatat                                                  <  1:3067613/93‑1 (MQ=255)

                                        cTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGtatata                                                 >  1:1756615/1‑93 (MQ=255)

                                            aTTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCg                                             >  1:5687831/1‑93 (MQ=255)

                                            aTTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCg                                             >  1:932494/1‑93 (MQ=255)

                                             ttGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGa                                            >  2:1759623/1‑93 (MQ=255)

                                             ttGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGa                                            >  1:2004008/1‑93 (MQ=255)

                                              tGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGaa                                           <  1:2618481/93‑1 (MQ=255)

                                               gCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAAc                                          <  2:1586554/93‑1 (MQ=255)

                                               gCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAAc                                          <  1:2635013/93‑1 (MQ=255)

                                                cGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACa                                         <  2:6296530/93‑1 (MQ=255)

                                                  tCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtg                                       >  1:5146132/1‑93 (MQ=255)

                                                      tggttggtATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtggtgg                                   >  2:5791624/1‑93 (MQ=255)

                                                       ggttggtATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGc                                  >  2:323859/1‑93 (MQ=255)

                                                       ggttggtATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGc                                  >  2:5940851/1‑93 (MQ=255)

                                                         ttggtATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGccc                                <  1:3346099/93‑1 (MQ=255)

                                                          tggtATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGcccc                               <  2:1841337/93‑1 (MQ=255)

                                                           ggtATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCg                              <  1:2131956/93‑1 (MQ=255)

                                                             tATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTg                            <  1:223110/93‑1 (MQ=255)

                                                             tATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTg                            >  1:6015391/1‑93 (MQ=255)

                                                               tGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCg                          <  2:3979320/93‑1 (MQ=255)

                                                               tGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCg                          >  1:4320092/1‑93 (MQ=255)

                                                               tGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCg                          >  2:4684944/1‑93 (MQ=255)

                                                                 ggACTTAACTGGCACCTGAATAAGCCCAAGTGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGc                        <  2:6190266/93‑1 (MQ=255)

                                                                 ggACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGc                        >  1:5375849/1‑93 (MQ=255)

                                                                 ggACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGc                        <  1:2407641/93‑1 (MQ=255)

                                                                 ggACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGc                        <  1:3380681/93‑1 (MQ=255)

                                                                 ggACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGc                        >  1:5200598/1‑93 (MQ=255)

                                                                   aCTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCaa                      <  2:927059/93‑1 (MQ=255)

                                                                    cTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAAc                     >  2:2959004/1‑93 (MQ=255)

                                                                    cTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAAc                     <  1:5499769/93‑1 (MQ=255)

                                                                     ttAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACt                    <  1:3994047/93‑1 (MQ=255)

                                                                     ttAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACt                    >  2:5532573/1‑93 (MQ=255)

                                                                     ttAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACt                    <  1:1584426/93‑1 (MQ=255)

                                                                     ttAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCACCGTGCGGCAACt                    >  2:1026831/1‑93 (MQ=255)

                                                                      tAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTg                   >  1:3238429/1‑93 (MQ=255)

                                                                       aaCTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGt                  <  1:5292406/93‑1 (MQ=255)

                                                                        aCTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGtt                 <  1:1109763/93‑1 (MQ=255)

                                                                           ggCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGAt              <  1:355042/93‑1 (MQ=255)

                                                                           ggCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGAt              <  2:3791060/93‑1 (MQ=255)

                                                                           ggCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGAt              <  1:1529051/93‑1 (MQ=255)

                                                                           ggCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGAt              >  2:3371049/1‑93 (MQ=255)

                                                                               ccTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCaa          >  1:3742024/1‑93 (MQ=255)

                                                                                cTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAAt         <  2:1970252/93‑1 (MQ=255)

                                                                                cTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAAt         <  2:4677982/93‑1 (MQ=255)

                                                                                 tGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATGCCCGATCATGGTGGCCCCGTGCGGCAACTGTTGATGCAATg        <  1:3536681/93‑1 (MQ=255)

                                                                                 tGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATg        >  1:1596042/1‑93 (MQ=255)

                                                                                 tGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATg        <  1:4859307/93‑1 (MQ=255)

                                                                                 tGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATg        <  2:1276312/93‑1 (MQ=255)

                                                                                 tGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATg        <  2:2081333/93‑1 (MQ=255)

                                                                                 tGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATg        >  2:224506/1‑93 (MQ=255)

                                                                                   aaTAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATgcg      <  2:699837/93‑1 (MQ=255)

                                                                                   aaTAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGCGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATgcg      >  2:2335182/1‑93 (MQ=255)

                                                                                    aTAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATgcgc     <  2:2111606/93‑1 (MQ=255)

                                                                                     tAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCa    <  2:1300945/93‑1 (MQ=255)

                                                                                      aaGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCTATGCGCAt   <  2:6118498/93‑1 (MQ=255)

                                                                                      aaGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCAt   >  1:1669664/1‑93 (MQ=255)

                                                                                       aGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATa  <  1:4013267/93‑1 (MQ=255)

GGATGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCTCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATA  >  NC_000913/3804456‑3804635

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