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.8%	1410.2 / ‑4.8	437	E50G (GAG→GGG)	rfaG	glucosyltransferase I
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (216/220); minor base G (0/1); total (216/221)
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.53e-01

Only 100 of 442 total aligned reads displayed.

GGATGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCTCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATA  >  NC_000913/3804456‑3804635

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

  aTGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCaaa                                                                                       <  1:6015942/93‑1 (MQ=255)

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

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

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

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

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

     ctctTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGc                                                                                    >  2:4997022/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

         ttGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCttta                                                                                <  2:7044830/93‑2 (MQ=255)

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

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

         ttGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTAAACTGGCACCTGAATAAGCCCAAATGCTTTc                                                                                >  2:3323622/1‑93 (MQ=255)

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

           gagaTGATTTTGTACCCAGGCATAATATCCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCgg                                                                              <  1:5257471/93‑1 (MQ=255)

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

            agaTGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCagg                                                                             >  1:817332/1‑93 (MQ=255)

               tgatTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCaa                                                                          >  2:2459357/1‑93 (MQ=255)

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

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

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

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

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

                    ttGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGcc                                                                     >  2:6306791/1‑93 (MQ=255)

                     tGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCt                                                                    <  2:5466398/93‑1 (MQ=255)

                     tGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCt                                                                    <  1:2054171/93‑1 (MQ=255)

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

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

                        aCCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTcc                                                                 <  2:618827/93‑1 (MQ=255)

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

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

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

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

                               cATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACt                                                          >  1:7089427/1‑93 (MQ=255)

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

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

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

                                 taataTTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACtgt                                                        <  1:4129645/93‑1 (MQ=255)

                                  aataTTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACtgtg                                                       <  2:1877446/93‑1 (MQ=255)

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

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

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

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

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

                                     atTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGtat                                                    >  1:3151986/1‑93 (MQ=255)

                                     atTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGtat                                                    >  1:1937601/1‑93 (MQ=255)

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

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

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

                                          gCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAAt                                                                         <  2:7839175/67‑1 (MQ=255)

                                          gCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAAt                                                                         >  1:7839176/1‑67 (MQ=255)

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

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

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

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

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

                                                   ccATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtgg                                      >  1:4038438/1‑93 (MQ=255)

                                                   ccATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtgg                                      >  2:7302208/1‑93 (MQ=255)

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

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

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

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

                                                     aTGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtggtg                                    <  2:4564998/93‑1 (MQ=255)

                                                     aTGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtggtg                                    <  2:1004454/93‑1 (MQ=255)

                                                     aTGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtggtg                                    <  2:1746991/93‑1 (MQ=255)

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

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

                                                         ttggtATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGccc                                >  1:1174354/1‑93 (MQ=255)

                                                              aTGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGc                           >  1:3334715/1‑93 (MQ=255)

                                                                 ggACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCGACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCgga                        >  1:5732890/1‑92 (MQ=255)

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

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

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

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

                                                                   aCTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCaa                      >  2:2261868/1‑93 (MQ=255)

                                                                   aCTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCaa                      >  2:2156966/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                                                                             cACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGc            <  1:5635709/93‑1 (MQ=255)

                                                                                cTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAAt         >  2:5582966/1‑93 (MQ=255)

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

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

                                                                                    aTAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGGAATgcgc     >  1:4773847/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` `< 30 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`