BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,722,918	T→A	N247Y (AAC→TAC)	b1647 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,722,918	0	T	A	100.0%	698.6 / NA	192	N247Y (AAC→TAC)	b1647	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); new base A (110/82); total (110/82)

Only 100 of 203 total aligned reads displayed.

TTTCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTTAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACA  >  NC_000913/1722828‑1723007

tttCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTaaa                                                                                         >  1:6003617/1‑93 (MQ=255)

tttCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTaaa                                                                                         >  1:6023959/1‑93 (MQ=255)

tttCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTaaa                                                                                         >  1:2070996/1‑93 (MQ=255)

tttCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTaaa                                                                                         >  2:6257330/1‑93 (MQ=255)

 ttCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAac                                                                                        >  1:5881285/1‑93 (MQ=255)

     gCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAaccacc                                                                                    <  1:1746615/93‑1 (MQ=255)

     gCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAaccacc                                                                                    >  2:3603940/1‑93 (MQ=255)

     gCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAaccacc                                                                                    >  2:3373998/1‑93 (MQ=255)

      cTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGCAATACCCAGGCGTAAACCAGCt                                                                                   >  2:894294/1‑93 (MQ=255)

       ttCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTg                                                                                  >  2:3314128/1‑93 (MQ=255)

        tCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGt                                                                                 <  2:5014339/93‑1 (MQ=255)

          gACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTc                                                                               >  2:3864490/1‑93 (MQ=255)

             aGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAAt                                                                            >  2:1582116/1‑93 (MQ=255)

              gCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATc                                                                           <  2:422728/93‑1 (MQ=255)

                tGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAACCCACCTGTTCAATCGa                                                                         <  2:4428083/93‑1 (MQ=255)

                tGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGa                                                                         >  2:3545788/1‑93 (MQ=255)

                tGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGa                                                                         <  1:2490438/93‑1 (MQ=255)

                 gCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAg                                                                        >  2:704973/1‑93 (MQ=255)

                 gCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAg                                                                        >  1:828246/1‑93 (MQ=255)

                  cccGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGc                                                                       >  1:2373254/1‑93 (MQ=255)

                  cccGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGc                                                                       <  1:2038092/93‑1 (MQ=255)

                   ccGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGcc                                                                      >  2:5413991/1‑93 (MQ=255)

                    cGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGccc                                                                     >  2:5788207/1‑93 (MQ=255)

                     gTACGCGCTCCATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGGTCAATCGAGcccc                                                                    >  1:4306659/1‑92 (MQ=255)

                     gTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCcg                                                                    >  1:4463705/1‑93 (MQ=255)

                      tACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCcgc                                                                   >  2:5820552/1‑93 (MQ=255)

                      tACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCcgc                                                                   <  1:5187917/93‑1 (MQ=255)

                        cgcgCTCAATTTTCCCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGt                                                                 <  1:3981714/93‑1 (MQ=255)

                        cgcgCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGt                                                                 >  2:6245879/1‑93 (MQ=255)

                        cgcgCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGt                                                                 >  2:5124008/1‑93 (MQ=255)

                         gcgcTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGtt                                                                <  1:2671875/93‑1 (MQ=255)

                           gcTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGCGCCCGCGATTa                                                              >  1:475224/1‑93 (MQ=255)

                           gcTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTa                                                              >  1:3990700/1‑93 (MQ=255)

                           gcTAAATTTCACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTa                                                              <  2:1185801/93‑1 (MQ=255)

                            cTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTaa                                                             >  2:2439056/1‑93 (MQ=255)

                            cTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTaa                                                             >  1:2674597/1‑93 (MQ=255)

                             tCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAAc                                                            >  1:3429392/1‑93 (MQ=255)

                             tCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAAc                                                            >  1:5379937/1‑93 (MQ=255)

                             tCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAAc                                                            <  1:2394190/93‑1 (MQ=255)

                                  tttACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActcctc                                                       <  1:1465719/93‑1 (MQ=255)

                                     aCCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGc                                                    >  1:5909653/1‑93 (MQ=255)

                                       cTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGcca                                                  <  2:4046481/93‑1 (MQ=255)

                                       cTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGcca                                                  >  1:5682569/1‑93 (MQ=255)

                                        tGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGccac                                                 >  2:1023583/1‑93 (MQ=255)

                                           aCCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCcaccac                                              <  1:5884916/93‑1 (MQ=255)

                                           aCCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCcaccac                                              <  2:5447310/93‑1 (MQ=255)

                                            ccGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACAACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACcaca                                             <  1:3383219/93‑1 (MQ=255)

                                              gATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGc                                           >  2:223797/1‑93 (MQ=255)

                                              gATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGc                                           <  2:1532441/93‑1 (MQ=255)

                                                tAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCa                                         >  2:2285842/1‑93 (MQ=255)

                                                 aaTTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGTCCGCGTTTAACTCCTCTTCCACCATATCCAg                                        >  2:6038900/1‑93 (MQ=255)

                                                       cAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtcatc                                  <  1:4660223/93‑1 (MQ=255)

                                                        aGCGGCTGCGATGGTAAACGTAATCCCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtcatcc                                 >  1:5719646/1‑92 (MQ=255)

                                                        aGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtcatca                                 >  1:401581/1‑93 (MQ=255)

                                                        aGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtcatca                                 >  1:5592227/1‑93 (MQ=255)

                                                        aGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtcatca                                 <  1:3319893/93‑1 (MQ=255)

                                                          cGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACg                               >  1:2737488/1‑93 (MQ=255)

                                                           ggCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGc                              >  1:2785698/1‑93 (MQ=255)

                                                           ggCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGc                              <  2:419761/93‑1 (MQ=255)

                                                            gCTGCGATGGTAAACGTAATACCCAGGCGTAAACTACCTGTTCAATCGAGCACGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCa                             >  1:3859578/1‑93 (MQ=255)

                                                            gCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCa                             >  2:4438204/1‑93 (MQ=255)

                                                            gCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCa                             <  2:1344454/93‑1 (MQ=255)

                                                             cTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAg                            <  2:2042219/93‑1 (MQ=255)

                                                             cTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAg                            <  1:2830081/93‑1 (MQ=255)

                                                             cTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAg                            >  2:4863756/1‑93 (MQ=255)

                                                               gCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCg                          >  2:3460066/1‑93 (MQ=255)

                                                               gCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCg                          <  1:1048983/93‑1 (MQ=255)

                                                                 gATGGTAAACGTAATACCCAGGCTTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGc                        <  2:5246295/93‑1 (MQ=255)

                                                                 gATGGTAAACGTAATACCCAGGCGTAATCCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGc                        <  1:2530511/93‑1 (MQ=255)

                                                                 gATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGc                        >  1:1333246/1‑93 (MQ=255)

                                                                  aTGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCt                       >  2:2856799/1‑93 (MQ=255)

                                                                   tGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTg                      <  1:98233/93‑1 (MQ=255)

                                                                   tGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTg                      >  1:1859363/1‑93 (MQ=255)

                                                                    ggTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTgg                     >  2:2792951/1‑93 (MQ=255)

                                                                    ggTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTgg                     >  1:5918465/1‑93 (MQ=255)

                                                                     gTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGa                    <  1:1263892/93‑1 (MQ=255)

                                                                     gTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGa                    >  1:3356841/1‑93 (MQ=255)

                                                                      tAAACGTAATACCCAGGCGTAACCCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTCCCACAAAAGCCAGTTCATCACGCAGCGGCTGGaa                   >  2:5142396/1‑93 (MQ=255)

                                                                      tAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGaa                   <  1:1705716/93‑1 (MQ=255)

                                                                      tAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGaa                   <  1:2176282/93‑1 (MQ=255)

                                                                        aaCGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGATTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGgaaa                 >  2:5559417/1‑92 (MQ=255)

                                                                         aCGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAaata                <  1:6306056/93‑1 (MQ=255)

                                                                          cGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAaataa               >  2:5648738/1‑93 (MQ=255)

                                                                           gTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAaataat              <  1:2220888/93‑1 (MQ=255)

                                                                           gTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAaataat              <  1:4134917/93‑1 (MQ=255)

                                                                            tAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAatc             >  1:6070789/1‑93 (MQ=255)

                                                                             aaTACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAatca            <  2:636314/93‑1 (MQ=255)

                                                                               tACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAatcatc          >  1:15455/1‑93 (MQ=255)

                                                                               tACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAatcatc          >  2:3028214/1‑93 (MQ=255)

                                                                                aCCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAAtcatca         >  2:401654/1‑93 (MQ=255)

                                                                                 cccAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATcatcat        <  1:5769761/93‑1 (MQ=255)

                                                                                  ccAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCAtt       >  2:3774419/1‑93 (MQ=255)

                                                                                   cAGGCGTAAACCACC‑GTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTaa     <  2:4761724/93‑1 (MQ=255)

                                                                                    aGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTaa     <  1:3109531/93‑1 (MQ=255)

                                                                                    aGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTaa     >  2:1741295/1‑93 (MQ=255)

                                                                                    aGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTaa     >  2:3768892/1‑93 (MQ=255)

                                                                                     tgCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTaaa    <  1:6208540/92‑1 (MQ=255)

                                                                                     ggCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTaaa    <  1:2091284/93‑1 (MQ=255)

                                                                                     ggCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTaaa    >  2:788095/1‑93 (MQ=255)

                                                                                      gCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAAc   <  1:4227345/93‑1 (MQ=255)

TTTCTGCTTCGACAGCTGCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTTAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACA  >  NC_000913/1722828‑1723007

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