BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	946,689	A→G	T144A (ACG→GCG)	b0899 →	putative transport

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	946,689	0	A	G	99.6%	851.6 / ‑4.8	252	T144A (ACG→GCG)	b0899	putative transport
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (123/128); total (123/129)
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.09e-01

Only 100 of 268 total aligned reads displayed.

CTCTATTTCACACCTTATGCGCTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACACGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATCTGGCACAAAAACCCCAGGCAATTCTG  >  NC_000913/946599‑946778

ctctATTTCACACCTTATGCGCTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACAcgcg                                                                                         <  1:2194282/93‑1 (MQ=255)

 tctATTTCACACCTTATGCGCTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAGGGGCGGGGTCAGTACCTGGATTAAACACGCGa                                                                                        <  1:7558981/93‑1 (MQ=255)

 tctATTTCACACCTTATGCGCTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGa                                                                                        <  1:6027814/93‑1 (MQ=255)

 tctATTTCACACCTTATGCGCTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGa                                                                                        >  1:5210283/1‑93 (MQ=255)

  ctATTTCACACCTTATGCGCTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGAt                                                                                       <  2:2161631/93‑1 (MQ=255)

   tATTTCACACCTTATGCGCTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATg                                                                                      >  1:4244926/1‑93 (MQ=255)

   tATTTCACACCTTATGCGCTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATg                                                                                      >  2:7273368/1‑93 (MQ=255)

   tATTTCACACCTTATGCGCTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACAAGCGATg                                                                                      >  1:4618621/1‑93 (MQ=255)

    aTTTCACACCTTATGCGCTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATgg                                                                                     <  2:1569792/93‑1 (MQ=255)

    aTTTCACACCTTATGCGCTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATgg                                                                                     >  1:3200978/1‑93 (MQ=255)

       tCACACCTTATGCGCTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTGAGTACCTGGATTAAACACGCGAGGGGAc                                                                                  >  1:7002787/1‑93 (MQ=255)

       tCACACCTTATGCGCTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGAc                                                                                  <  2:1396734/93‑1 (MQ=255)

       tCACAACTTATGCGCTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGAc                                                                                  <  1:4466979/93‑1 (MQ=255)

         acacCTTATGCGCTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGAccc                                                                                >  1:6337391/1‑93 (MQ=255)

         acacCTTATGCGCTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGAccc                                                                                <  2:6515202/93‑1 (MQ=255)

         acacCTTATGCGCTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGAccc                                                                                >  1:3936102/1‑93 (MQ=255)

         acacCTTATGCGCTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGAccc                                                                                >  2:186357/1‑93 (MQ=255)

          cacCTTATGCGCTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCg                                                                               >  1:4351522/1‑93 (MQ=255)

          cacCTTATGCGCTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCAGGGTCATTACCTGGATTAAACACGCCATGGGACCCg                                                                               <  2:4148698/93‑1 (MQ=255)

           acCTTATGCGCTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCgg                                                                              <  1:7435316/93‑1 (MQ=255)

            ccTTATGCGCTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGa                                                                             <  1:7745698/93‑1 (MQ=255)

             cTTATGCGCTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGAc                                                                            <  2:8302919/93‑1 (MQ=255)

               tATGCGCTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTg                                                                          >  1:6163308/1‑93 (MQ=255)

                 tGCGCTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGc                                                                        <  2:2452942/93‑1 (MQ=255)

                   cgcTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGGCAGTAACTGGATTAAACACGGGATGGGACCCGGACTGTCtt                                                                      >  1:2334476/1‑93 (MQ=255)

                    gcTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCtga                                                                     >  1:1591621/1‑91 (MQ=255)

                    gcTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTa                                                                     >  2:8448804/1‑93 (MQ=255)

                     cTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTAt                                                                    >  1:6997526/1‑93 (MQ=255)

                      tAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGTGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTAtc                                                                   >  1:4399113/1‑93 (MQ=255)

                       aaTTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTAtct                                                                  <  1:3999108/93‑1 (MQ=255)

                       aaTTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTAtct                                                                  >  2:7898583/1‑93 (MQ=255)

                       aaTTGTTCGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTAtct                                                                  <  1:2893833/93‑1 (MQ=255)

                        aTTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTAtctc                                                                 <  1:8755652/93‑1 (MQ=255)

                          tgttgGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTcgc                                                               >  2:2083027/1‑93 (MQ=255)

                          tgttgGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTcgc                                                               <  1:2721120/93‑1 (MQ=255)

                           gttgGTGAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTcgcc                                                              <  1:1729708/93‑1 (MQ=255)

                             tgGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTcgccgc                                                            >  1:6686867/1‑93 (MQ=255)

                              ggTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCcgcg                                                           <  1:1153076/93‑1 (MQ=255)

                              ggTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCcgcg                                                           >  2:6717706/1‑93 (MQ=255)

                                 cAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTgg                                                        <  1:7240260/93‑1 (MQ=255)

                                  aGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTggc                                                       >  2:2228962/1‑92 (MQ=255)

                                  aGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGa                                                       <  2:5920566/93‑1 (MQ=255)

                                   gTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGAc                                                      >  1:4723348/1‑93 (MQ=255)

                                    ttAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGAcc                                                     >  2:4049712/1‑93 (MQ=255)

                                    ttAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGAcc                                                     <  2:3188900/93‑1 (MQ=255)

                                     tAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCt                                                    >  1:8543199/1‑93 (MQ=255)

                                      aGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTa                                                   >  2:1888774/1‑93 (MQ=255)

                                       ggCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTAc                                                  <  1:5691686/93‑1 (MQ=255)

                                         cTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTg                                                <  2:455449/93‑1 (MQ=255)

                                         cTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTg                                                <  1:260967/93‑1 (MQ=255)

                                          tCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTgg                                               <  2:1032237/93‑1 (MQ=255)

                                            gACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGggt                                             <  1:4924845/93‑1 (MQ=255)

                                            gACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGggt                                             <  2:5416931/93‑1 (MQ=255)

                                            gACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGggt                                             <  2:2599169/93‑1 (MQ=255)

                                              ccTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGggtgg                                           <  1:5292927/93‑1 (MQ=255)

                                               cTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGggtggt                                          <  2:2988895/93‑1 (MQ=255)

                                                ttCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGgtggtg                                         >  2:4957363/1‑93 (MQ=255)

                                                ttCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGgtggtg                                         >  2:3811891/1‑93 (MQ=255)

                                                 tCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGc                                        >  2:5677061/1‑93 (MQ=255)

                                                 tCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGc                                        >  2:5174914/1‑93 (MQ=255)

                                                  cAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCCACTGGGTGGTGCa                                       >  2:3394117/1‑93 (MQ=255)

                                                   aaaGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCat                                      <  2:8211691/93‑1 (MQ=255)

                                                   aaaGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCat                                      >  2:936619/1‑93 (MQ=255)

                                                   aaaGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCat                                      >  2:1515160/1‑93 (MQ=255)

                                                    aaGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCata                                     >  1:5077184/1‑93 (MQ=255)

                                                    aaGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCata                                     >  2:3619481/1‑93 (MQ=255)

                                                     aGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCatat                                    <  1:5780368/93‑1 (MQ=255)

                                                     aGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCatat                                    >  2:4744582/1‑93 (MQ=255)

                                                      gATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATAtt                                   <  1:6599399/93‑1 (MQ=255)

                                                           gAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTa                              <  2:1901787/93‑1 (MQ=255)

                                                           gAAGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTa                              <  1:3855266/93‑1 (MQ=255)

                                                            aaGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTAt                             >  1:3584586/1‑93 (MQ=255)

                                                             aGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATc                            <  2:6548141/93‑1 (MQ=255)

                                                             aGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATc                            <  2:3508972/93‑1 (MQ=255)

                                                             aGGGCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATc                            >  2:4932937/1‑93 (MQ=255)

                                                              gggCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATCt                           >  2:5387178/1‑93 (MQ=255)

                                                              gggCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATCt                           <  2:7967347/93‑1 (MQ=255)

                                                               ggCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATCTg                          >  1:6462602/1‑93 (MQ=255)

                                                               ggCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATCTg                          <  2:8531539/93‑1 (MQ=255)

                                                               ggCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATCTg                          <  1:7392078/93‑1 (MQ=255)

                                                                gCGGGGTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATCTgg                         <  2:2002920/93‑1 (MQ=255)

                                                                   gggTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATCTGGcac                      >  1:4968985/1‑93 (MQ=255)

                                                                   gggTCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATCTGGcac                      >  2:6859653/1‑93 (MQ=255)

                                                                      tCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATCTGGCACaaa                   >  2:6347980/1‑93 (MQ=255)

                                                                      tCAGTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATCTGGCACaaa                   <  1:2999407/93‑1 (MQ=255)

                                                                         gTACCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATCTGGCACAAAAAc                <  2:1560370/93‑1 (MQ=255)

                                                                           aCCTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATCTGGCACAAAAAccc              >  2:1354403/1‑93 (MQ=255)

                                                                            ccTGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATCTGGCACAAAAAcccc             <  1:7513374/93‑1 (MQ=255)

                                                                             cTGGATTAAACACGCGTTGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGTTGCACATTCCCTATCTGGCACAAAAACCCCa            >  2:2784931/1‑93 (MQ=255)

                                                                              tGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATCTGGCACAAAAACCCCAg           <  1:1078019/93‑1 (MQ=255)

                                                                              tGGATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATCTGGCACAAAAACCCCAg           <  1:3237469/93‑1 (MQ=255)

                                                                               ggATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATCTGGCACAAAAACCCCAgg          >  1:5359/1‑93 (MQ=255)

                                                                               ggATTAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATCTGGCACAAAAACCCCAgg          >  1:7302845/1‑93 (MQ=255)

                                                                                  ttAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATCTGGCACAAAAACCCCAGGCaa       >  1:1033786/1‑93 (MQ=255)

                                                                                  ttAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATCTGGCACAAAAACCCCAGGCaa       >  1:588605/1‑93 (MQ=255)

                                                                                   tAAACACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATCTGGCACAAAAACCCCAGGCAAt      <  1:1882787/93‑1 (MQ=255)

                                                                                   tAAACACGCGATGGGACCCGGACGGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATCTGGCACAAAAACCCCAGGCAAt      <  2:5977893/93‑1 (MQ=255)

                                                                                    aaaCACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATCTGGCACAAAAACCCCAGGCAAtt     <  1:8299037/93‑1 (MQ=255)

                                                                                     aaCACGCGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATCTGGCACAAAAAACCCAGGCAATTc    >  2:7573415/1‑93 (MQ=255)

                                                                                      acacgcGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTCGGCGGTGCATCTTCCCTATCTGGCACAAAAACCCCACGCAATTCt   >  1:6288951/1‑93 (MQ=255)

CTCTATTTCACACCTTATGCGCTAATTGTTGGTCAGTTAGGCTCGACCTTCAAAGATGGGAAGGGCGGGGTCAGTACCTGGATTAAACACACGATGGGACCCGGACTGGCTTATCTCGCCGCGTGGACCTACTGGGTGGTGCATATTCCCTATCTGGCACAAAAACCCCAGGCAATTCTG  >  NC_000913/946599‑946778

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