BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,160,251	A→G	T145T (ACT→ACC)	sufI ←	suppressor of ftsI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,160,251	0	A	G	99.5%	579.7 / NA	184	T145T (ACT→ACC)	sufI	suppressor of ftsI
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (79/104); minor base T (1/0); total (80/104)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 194 total aligned reads displayed.

GTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCG  >  NC_000913/3160161‑3160339

gTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGt                                                                                        >  1:27254/1‑93 (MQ=255)

  tGGGGATAGGCAGCGTCTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTAt                                                                                      <  1:900192/93‑1 (MQ=255)

    gggATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTg                                                                                    <  1:2489841/93‑1 (MQ=255)

    gggATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTg                                                                                    <  2:2635352/93‑1 (MQ=255)

      gATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGc                                                                                  >  1:1180424/1‑93 (MQ=255)

      gATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGc                                                                                  >  2:1848156/1‑93 (MQ=255)

       aTAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTAAGCCGTGCGGTTGGGGGTATTGGCg                                                                                 >  2:780033/1‑93 (MQ=255)

           gCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGAt                                                                             <  2:3551503/93‑1 (MQ=255)

            cAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATa                                                                            <  1:284038/93‑1 (MQ=255)

            cAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATa                                                                            >  2:2877335/1‑93 (MQ=255)

             aGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATAc                                                                           <  1:3404889/93‑1 (MQ=255)

               cGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACca                                                                         <  2:1374014/93‑1 (MQ=255)

               cGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACca                                                                         <  1:1720141/93‑1 (MQ=255)

                gACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACcac                                                                        >  1:1448799/1‑93 (MQ=255)

                 aCTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACcaca                                                                       >  1:3297168/1‑93 (MQ=255)

                      cTGATTTCATCTGCCACCAGCCACCTTCCGGCAAGGCCGTTATAGACGTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATTCCACAGATTa                                                                  <  2:1923577/93‑1 (MQ=255)

                      cTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTa                                                                  <  1:3555720/93‑1 (MQ=255)

                       tGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAg                                                                 >  2:4631164/1‑93 (MQ=255)

                       tGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAg                                                                 >  2:1735034/1‑93 (MQ=255)

                         aCTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCt                                                               <  1:547252/93‑1 (MQ=255)

                         aCTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCt                                                               >  1:641392/1‑93 (MQ=255)

                           ttCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGc                                                             >  1:1061963/1‑93 (MQ=255)

                           ttCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGc                                                             >  1:194781/1‑93 (MQ=255)

                             cATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGt                                                           <  2:2931658/93‑1 (MQ=255)

                               tCTTCCACCAGCCACATTCCTGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTc                                                         <  2:1412580/93‑1 (MQ=255)

                               tCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTc                                                         <  1:1391015/93‑1 (MQ=255)

                               tCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTc                                                         <  1:3980972/93‑1 (MQ=255)

                                cTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCt                                                        <  2:1435046/93‑1 (MQ=255)

                                cTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCt                                                        <  1:1645943/93‑1 (MQ=255)

                                 ttCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTg                                                       <  1:2219504/93‑1 (MQ=255)

                                 ttCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTg                                                       <  1:4457767/93‑1 (MQ=255)

                                 ttCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTg                                                       >  1:1843281/1‑93 (MQ=255)

                                 ttCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTg                                                       >  2:4654570/1‑93 (MQ=255)

                                    caccaGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCTGCGTGGTATGGGGTATTGGCGTGATATCACAGAGTATCTGCGTTCTGGCg                                                    >  2:4685758/1‑93 (MQ=255)

                                      ccaGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGaa                                                  <  1:410137/93‑1 (MQ=255)

                                       caGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAAt                                                 >  2:3464546/1‑93 (MQ=255)

                                       caGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAAt                                                 <  2:2854442/93‑1 (MQ=255)

                                        aGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATg                                                <  1:2675844/93‑1 (MQ=255)

                                         gCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATgg                                               <  1:152171/93‑1 (MQ=255)

                                          ccACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATggg                                              <  2:1465718/93‑1 (MQ=255)

                                          ccACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATggg                                              <  1:1976435/93‑1 (MQ=255)

                                           cacaTTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGc                                             <  1:1128811/93‑1 (MQ=255)

                                           cacaTTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGc                                             >  2:1457794/1‑93 (MQ=255)

                                           cacaTTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACAACAGAGTAGCTGCGTTCTGGCGAATGGGc                                             >  1:156069/1‑93 (MQ=255)

                                            acaTTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCa                                            <  2:161230/93‑1 (MQ=255)

                                            acaTTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCa                                            >  1:1314262/1‑93 (MQ=255)

                                            acaTTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCa                                            >  1:2055615/1‑93 (MQ=255)

                                            acaTTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGGTGGGGGCATTGGCGTGAGATGACAGAGGAGCTGCGCTCTAGCGAATgggaa                                            >  2:3631841/1‑91 (MQ=255)

                                             caTTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAg                                           >  2:4107281/1‑93 (MQ=255)

                                             caTTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAg                                           <  1:2992353/93‑1 (MQ=255)

                                             caTTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAg                                           <  2:1283859/93‑1 (MQ=255)

                                              aTTCCGGAAAGGCCGTTATAGACCTGCTGAGAAGTGCGGTTGGGGGTATTGGCGGGATAGTACAGTGTAGCTGCGTTCTGGCGACTGGGCAGt                                          >  2:3025783/1‑93 (MQ=255)

                                               ttCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTa                                         <  1:1008573/93‑1 (MQ=255)

                                                   ggCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACggg                                     <  2:1028747/93‑1 (MQ=255)

                                                   ggCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACggg                                     >  1:1173365/1‑93 (MQ=255)

                                                    gCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGt                                    <  1:2585262/93‑1 (MQ=255)

                                                    gCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGt                                    <  2:1199965/93‑1 (MQ=255)

                                                       aGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGcc                                 <  2:3392096/93‑1 (MQ=255)

                                                       aGGCCGTTATAGACCTGCTAAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGcc                                 >  1:3085576/1‑93 (MQ=255)

                                                        ggCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGccc                                <  2:1845922/93‑1 (MQ=255)

                                                          ccGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAg                              >  2:2057027/1‑93 (MQ=255)

                                                           cGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGt                             <  2:3332648/93‑1 (MQ=255)

                                                           cGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGt                             >  2:893334/1‑93 (MQ=255)

                                                           cGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGt                             <  2:671841/93‑1 (MQ=255)

                                                           cGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGt                             >  2:1104037/1‑93 (MQ=255)

                                                            gTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTc                            <  1:1635941/93‑1 (MQ=255)

                                                            gTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTc                            >  1:680862/1‑93 (MQ=255)

                                                             ttATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCa                           <  1:276385/93‑1 (MQ=255)

                                                             ttATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCa                           <  1:4650361/93‑1 (MQ=255)

                                                             ttATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCa                           <  2:3070949/93‑1 (MQ=255)

                                                              tataGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAg                          >  2:1672989/1‑93 (MQ=255)

                                                              tataGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAg                          >  1:1258083/1‑93 (MQ=255)

                                                              tataGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAg                          >  1:3775640/1‑93 (MQ=255)

                                                               ataGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGc                         >  1:1965789/1‑93 (MQ=255)

                                                               ataGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGc                         >  2:3184968/1‑93 (MQ=255)

                                                                taGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCg                        <  2:1489264/93‑1 (MQ=255)

                                                                  gACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCggt                      >  2:123184/1‑91 (MQ=255)

                                                                  gACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGtt                      <  1:2201947/93‑1 (MQ=255)

                                                                  gACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGtt                      >  1:2420818/1‑93 (MQ=255)

                                                                  gACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGtt                      >  1:587135/1‑93 (MQ=255)

                                                                   aCCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGttt                     >  2:2962064/1‑93 (MQ=255)

                                                                    ccTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTg                    <  1:3241942/93‑1 (MQ=255)

                                                                     ctgctgAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTgg                   >  1:40352/1‑93 (MQ=255)

                                                                      tgctgAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGc                  <  2:4324427/93‑1 (MQ=255)

                                                                          gAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGAca              <  1:4365500/93‑1 (MQ=255)

                                                                          gAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGAca              >  1:3914019/1‑93 (MQ=255)

                                                                           aGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGAcat             <  2:3967118/93‑1 (MQ=255)

                                                                            gCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGAcatc            >  1:2569137/1‑93 (MQ=255)

                                                                              cGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGAcatcat          <  1:793444/93‑1 (MQ=255)

                                                                              cGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGAcatcat          <  2:3047797/93‑1 (MQ=255)

                                                                               gTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACatcatc         <  2:2436918/93‑1 (MQ=255)

                                                                               gTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACatcatc         <  2:3838641/93‑1 (MQ=255)

                                                                               gTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACatcatc         >  1:2382535/1‑93 (MQ=255)

                                                                                   ggTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTg     >  2:3511783/1‑93 (MQ=255)

                                                                                   ggTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTCCCCAGTCAGCGTTTGGCGACATCATCCGTg     <  1:718239/93‑1 (MQ=255)

                                                                                   gcttGGTAGTAGTGGCGTGATACCACAGAGTAGCTGCTTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTg     <  2:3122396/91‑1 (MQ=255)

                                                                                    gTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGc    <  2:3780805/93‑1 (MQ=255)

                                                                                     ttGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGcc   >  2:4272983/1‑93 (MQ=255)

                                                                                     ttGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGcc   >  2:2086179/1‑93 (MQ=255)

                                                                                      tGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCg  >  1:4241335/1‑93 (MQ=255)

GTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCG  >  NC_000913/3160161‑3160339

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