BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,230,972	A→G	G188G (GGA→GGG)	umuC →	UmuC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,230,972	0	A	G	99.2%	388.8 / ‑4.8	123	G188G (GGA→GGG)	umuC	UmuC
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (66/56); total (67/56)
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.58e-01

Only 100 of 131 total aligned reads displayed.

GCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGACGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTT  >  NC_000913/1230882‑1231061

gCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGcg                                                                                         >  2:1542226/1‑93 (MQ=255)

 tccaCGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGcgg                                                                                        >  2:4219614/4‑93 (MQ=255)

  aGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGggcggc                                                                                       <  1:3799485/93‑1 (MQ=255)

  aGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGggcggc                                                                                       <  2:1040346/93‑1 (MQ=255)

     cGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGa                                                                                    <  1:690266/93‑1 (MQ=255)

       ggtggggtggTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATc                                                                                  >  2:3891739/1‑93 (MQ=255)

        gtggggtggTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCa                                                                                 <  1:2702563/93‑1 (MQ=255)

        gtggggtggTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCa                                                                                 <  2:1225810/93‑1 (MQ=255)

        gtggggtggTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCa                                                                                 >  1:3507869/1‑93 (MQ=255)

           gggtggTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCACCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCa                                                                              <  1:89930/93‑1 (MQ=255)

              tggtggATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCaaaa                                                                           >  2:4545112/1‑93 (MQ=255)

                gtggATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCaaaaaa                                                                         <  1:354211/93‑1 (MQ=255)

                gtggATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCaaaaaa                                                                         >  2:3659094/1‑93 (MQ=255)

                  ggATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACt                                                                       >  2:2239694/1‑93 (MQ=255)

                  ggATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACt                                                                       >  2:3860896/1‑93 (MQ=255)

                   gATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTg                                                                      <  2:1250096/93‑1 (MQ=255)

                     tttATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGa                                                                    <  2:4709540/93‑1 (MQ=255)

                     tttATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGa                                                                    >  1:3011841/1‑93 (MQ=255)

                     tgtATCAAATCTGGACCGCCAGCGTAATTTAATGTCTGCTCTCCCCGGGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGa                                                                    <  2:4709612/91‑1 (MQ=255)

                      ttATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGAc                                                                   <  1:3647399/93‑1 (MQ=255)

                      ttATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGAc                                                                   <  1:355494/93‑1 (MQ=255)

                       tATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGAcg                                                                  <  2:3592070/93‑1 (MQ=255)

                         tCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGAcgcg                                                                >  2:64914/1‑93 (MQ=255)

                         tCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGAcgcg                                                                >  1:4691922/1‑93 (MQ=255)

                         tCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGAcgcg                                                                >  1:3114323/1‑93 (MQ=255)

                         tCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGCGACTGGGCGGCGGATCAGCAAAAAACTGGAcgcg                                                                >  2:64915/1‑93 (MQ=255)

                          cAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGa                                                               <  1:1961186/93‑1 (MQ=255)

                          cAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGa                                                               <  1:1961185/93‑1 (MQ=255)

                           aaaTCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGAt                                                              >  1:3348401/1‑93 (MQ=255)

                            aaTCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATg                                                             <  2:1397469/93‑1 (MQ=255)

                            aaTCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATg                                                             >  2:1502060/1‑93 (MQ=255)

                             aTCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATgg                                                            <  1:3803015/93‑1 (MQ=255)

                              tCTGGAACGCCAGCGTAAATTAATGTCTGCTCCCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATggg                                                           <  1:1049561/93‑1 (MQ=255)

                               cTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATgggg                                                          <  2:356380/93‑1 (MQ=255)

                               cTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATgggg                                                          <  2:2567760/93‑1 (MQ=255)

                                tGGACCGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATTACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGa                                                         <  2:180381/93‑1 (MQ=255)

                                tGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGa                                                         <  1:2567420/93‑1 (MQ=255)

                                  gAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATc                                                       >  1:356458/1‑93 (MQ=255)

                                    aCGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCaa                                                     >  2:1202277/1‑93 (MQ=255)

                                    aCGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCaa                                                     >  1:2366088/1‑93 (MQ=255)

                                      gCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCaaaa                                                   >  2:2056633/1‑93 (MQ=255)

                                      gCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCaaaa                                                   >  2:1383911/1‑93 (MQ=255)

                                       ccAGTGTAAATTAATGTATGCCCTACCCGTAGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATAAAAAc                                                  <  2:3844612/93‑1 (MQ=255)

                                       ccAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAAc                                                  >  2:211628/1‑93 (MQ=255)

                                          gCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGt                                               <  1:3980361/93‑1 (MQ=255)

                                          gCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGt                                               >  1:628649/1‑93 (MQ=255)

                                           cGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGtt                                              >  1:1001598/1‑93 (MQ=255)

                                            gTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTtc                                             <  2:57784/93‑1 (MQ=255)

                                               aaTTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCg                                          >  1:200678/1‑93 (MQ=255)

                                                aTTAATGTCTGCTCTCCCCGTGGATGACGTCCGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGa                                         >  1:2193647/1‑93 (MQ=255)

                                                 ttAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGAt                                        >  2:3571841/1‑93 (MQ=255)

                                                 ttAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGAt                                        <  1:2384440/93‑1 (MQ=255)

                                                  tAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGAtt                                       >  2:3092495/1‑93 (MQ=255)

                                                   aaTGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAACACCGTTCTCGAttt                                      >  2:2547755/1‑93 (MQ=255)

                                                   aaTGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGAttt                                      >  2:269978/1‑93 (MQ=255)

                                                   aaTGTCTGCTCTCCACGTGGATGACGTCAGGGGGATTGGGCGGCGGATCAGCAAAAACCTGGACGCGATGGGGATCAAAACCGTTCTCGAttt                                      <  1:2930020/93‑1 (MQ=255)

                                                    aTGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTg                                     <  1:1054665/93‑1 (MQ=255)

                                                    aTGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTg                                     <  2:163367/93‑1 (MQ=255)

                                                     tGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGTTTTgg                                    >  1:1210495/1‑93 (MQ=255)

                                                      gTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTggg                                   >  2:1438712/1‑92 (MQ=255)

                                                      gTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGACGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGc                                   >  2:4016237/1‑93 (MQ=255)

                                                      gTCTGCTCGCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGTATCAAAACCGTTCTCGATTTGGc                                   <  2:3904209/93‑1 (MQ=255)

                                                       tCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCg                                  >  1:1883337/1‑93 (MQ=255)

                                                       tCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCg                                  >  1:923412/1‑93 (MQ=255)

                                                         tGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGa                                >  1:1946304/1‑93 (MQ=255)

                                                         tGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGa                                <  1:4222351/93‑1 (MQ=255)

                                                         tGCTCTCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGa                                <  1:536827/93‑1 (MQ=255)

                                                           ctctCCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATa                              <  1:1934422/93‑1 (MQ=255)

                                                            tctcCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATAc                             >  2:3561073/1‑93 (MQ=255)

                                                             ctcCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACa                            >  2:2520875/1‑93 (MQ=255)

                                                              tcCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGCCGCGGTGGGGATACAAACCGTTCTCGCTTTGGCGGATACAg                           >  1:1400659/1‑93 (MQ=255)

                                                              tcCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCTATTTGGCGGATACAg                           >  1:1400660/1‑93 (MQ=255)

                                                              tcCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAg                           >  1:2949778/1‑93 (MQ=255)

                                                              tcCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAg                           <  2:3625740/93‑1 (MQ=255)

                                                              tcCCCGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAg                           >  2:4153473/1‑93 (MQ=255)

                                                               ccccGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGa                          <  1:3841693/93‑1 (MQ=255)

                                                                 ccGTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGata                        >  1:4563388/1‑93 (MQ=255)

                                                                   gTGGATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATc                      <  2:304158/93‑1 (MQ=255)

                                                                     ggATGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCg                    <  1:2848085/93‑1 (MQ=255)

                                                                        tGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAATCCGTTCTCGATTTGGCGGATACAGATATCCGGtt                 <  2:203348/93‑1 (MQ=255)

                                                                        tGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGtt                 <  2:3031832/93‑1 (MQ=255)

                                                                        tGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGtt                 >  2:4622604/1‑93 (MQ=255)

                                                                        tGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGtt                 >  2:4622605/1‑93 (MQ=255)

                                                                        tGACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGtt                 <  1:1522606/93‑1 (MQ=255)

                                                                         gACGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGttt                <  2:3121932/93‑1 (MQ=255)

                                                                          aCGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTa               >  1:3057680/1‑93 (MQ=255)

                                                                           cGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTac              >  1:3659372/1‑92 (MQ=255)

                                                                           cGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTAt              >  1:3659370/1‑93 (MQ=255)

                                                                           cGTCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTAt              <  2:1068546/93‑1 (MQ=255)

                                                                             tCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCTGTTTATcc            >  2:2023871/1‑93 (MQ=255)

                                                                             tCTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATcc            <  1:2659843/93‑1 (MQ=255)

                                                                              cTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCg           >  2:2859923/1‑93 (MQ=255)

                                                                              cTGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCg           <  2:1021609/93‑1 (MQ=255)

                                                                               tGGGGGATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGt          >  2:933780/1‑93 (MQ=255)

                                                                               tGGGGGATTGGGCGACGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGt          >  2:1533079/1‑93 (MQ=255)

                                                                                gggggATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTa         >  1:210896/1‑93 (MQ=255)

                                                                                gggggATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTa         <  1:696176/93‑1 (MQ=255)

                                                                                  gggATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTaaa       <  2:3109906/93‑1 (MQ=255)

                                                                                    gATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACa     <  1:3412975/93‑1 (MQ=255)

                                                                                    gATTGGGCGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACa     <  1:3413091/93‑1 (MQ=255)

GCAGACGGGTGGGGTGGTGGATTTATCAAATCTGGAACGCCAGCGTAAATTAATGTCTGCTCTCCCCGTGGATGACGTCTGGGGGATTGGACGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTT  >  NC_000913/1230882‑1231061

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