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.0%	687.2 / ‑3.1	210	E50G (GAG→GGG)	rfaG	glucosyltransferase I
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (97/110); total (99/111)
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.60e-01

Only 100 of 213 total aligned reads displayed.

GGATGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCTCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATA  >  NC_000913/3804456‑3804635

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

  aTGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCaaa                                                                                       >  1:3956916/1‑93 (MQ=255)

  aTGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCaaa                                                                                       >  2:453748/1‑93 (MQ=255)

   tGCTCTTTGAGATGATTTTGTACCCAGGCCTAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCACAAAt                                                                                      >  1:2887141/1‑93 (MQ=255)

   tGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAAt                                                                                      >  1:3809923/1‑93 (MQ=255)

   tGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAAt                                                                                      <  1:3391132/93‑1 (MQ=255)

   tGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAAt                                                                                      >  2:4675683/1‑93 (MQ=255)

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

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

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

        tttGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCttt                                                                                 >  1:1930110/1‑93 (MQ=255)

        tttGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCttt                                                                                 >  2:3728330/1‑93 (MQ=255)

        tttGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCttt                                                                                 >  1:2765190/1‑93 (MQ=255)

         ttGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCTCAAATGCTTTc                                                                                >  2:4740664/1‑93 (MQ=255)

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

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

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

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

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

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

            agaTGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCggg                                                                             <  2:3839562/93‑1 (MQ=255)

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

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

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

              atgatTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCa                                                                           >  1:3469678/1‑93 (MQ=255)

               tgatTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCaa                                                                          <  2:1087237/93‑1 (MQ=255)

                gatTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGTACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAAt                                                                         <  1:2958266/93‑1 (MQ=255)

                gatTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAAt                                                                         >  1:238360/1‑93 (MQ=255)

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

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

                  ttttGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCg                                                                       >  1:1068981/1‑93 (MQ=255)

                    ttGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGTTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGcc                                                                     <  1:4553639/93‑1 (MQ=255)

                    ttGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGcc                                                                     <  2:3326347/93‑1 (MQ=255)

                      gTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCtt                                                                   <  2:3262368/93‑1 (MQ=255)

                      gTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCtt                                                                   >  2:4299336/1‑93 (MQ=255)

                       tACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTc                                                                  >  2:770404/1‑93 (MQ=255)

                         cccAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTccc                                                                <  1:412962/93‑1 (MQ=255)

                          ccAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCa                                                               >  1:1405264/1‑93 (MQ=255)

                          ccAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCa                                                               >  1:1538989/1‑93 (MQ=255)

                           cAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCAc                                                              >  2:3355806/1‑93 (MQ=255)

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

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

                              gCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGAc                                                           <  2:4449876/93‑1 (MQ=255)

                              gCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGAc                                                           <  1:943639/93‑1 (MQ=255)

                              gCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGAc                                                           <  1:210935/93‑1 (MQ=255)

                                ataataTTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACtg                                                         >  1:1815408/1‑93 (MQ=255)

                                ataataTTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACtg                                                         <  2:2966049/93‑1 (MQ=255)

                                ataataTTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACtg                                                         <  1:2246762/93‑1 (MQ=255)

                                ataataTTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACtg                                                         >  2:1160864/1‑93 (MQ=255)

                                ataataTTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACtg                                                         >  2:831115/1‑93 (MQ=255)

                                 taataTTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACtgt                                                        >  2:2534610/1‑93 (MQ=255)

                                 taataTTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCGTCCCACGACtgt                                                        >  2:1927501/1‑93 (MQ=255)

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

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

                                    tatTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGta                                                     >  2:4377870/1‑93 (MQ=255)

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

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

                                      tgcTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGtata                                                   <  2:1299238/91‑1 (MQ=255)

                                      atCTTCATTGCGTCCATTGTTGGTATGTGATTTAACTCGCACCTGAATAAGCTCAAATGCTTTCTTTCAATCGCCTTCCCACGACTGTGtata                                                   <  1:4021968/92‑1 (MQ=255)

                                       tCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGtatat                                                  >  1:888968/1‑93 (MQ=255)

                                        cTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGtatata                                                 <  2:2859432/93‑1 (MQ=255)

                                         tGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATAc                                                >  2:1494747/1‑93 (MQ=255)

                                          gCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATAcc                                               <  2:3440831/93‑1 (MQ=255)

                                          gCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATAcc                                               <  2:2723230/93‑1 (MQ=255)

                                           cATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATAccc                                              <  2:2850209/93‑1 (MQ=255)

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

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

                                                 gTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAt                                        >  1:1886529/1‑93 (MQ=255)

                                                 gTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAt                                        <  1:1010809/93‑1 (MQ=255)

                                                  tCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtg                                       <  2:1783671/93‑1 (MQ=255)

                                                   ccATGGTTGGTATGGGACTTAACTGGCACCTTAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtgg                                      <  2:2080684/93‑1 (MQ=255)

                                                   ccATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtgg                                      <  2:873856/93‑1 (MQ=255)

                                                   ccATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtgg                                      <  1:3210761/93‑1 (MQ=255)

                                                    cATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtggt                                     >  2:1999706/1‑93 (MQ=255)

                                                    cATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtggt                                     <  2:1526195/93‑1 (MQ=255)

                                                    cATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtggt                                     <  1:2185606/93‑1 (MQ=255)

                                                    cATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtggt                                     >  2:2818505/1‑93 (MQ=255)

                                                         ttggtATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACTTGGTGGccc                                >  2:3611351/1‑93 (MQ=255)

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

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

                                                                 ggACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCACGTGCGGc                        >  2:4097740/1‑93 (MQ=255)

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

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

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

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

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

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

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

                                                                         cTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTg                <  2:660396/93‑1 (MQ=255)

                                                                          tGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGa               <  2:3770830/93‑1 (MQ=255)

                                                                          tGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGa               <  1:3569757/93‑1 (MQ=255)

                                                                          tGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGa               >  1:4288683/1‑93 (MQ=255)

                                                                          tGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATCGTGGCCCCGTGCGGCAACTGTTGa               >  1:1970649/1‑93 (MQ=255)

                                                                             cACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGACCCCGTGCGGCAACTGTTGATGc            >  2:3244232/1‑93 (MQ=255)

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

                                                                               ccTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCaa          <  1:3076361/93‑1 (MQ=255)

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

                                                                                   aaTAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATgcg      >  1:1087146/1‑93 (MQ=255)

                                                                                     tAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCa    >  2:2460740/1‑93 (MQ=255)

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

GGATGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCTCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATA  >  NC_000913/3804456‑3804635

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