BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,493,567	T→C	F86L (TTT→CTT)	b1423 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,493,567	0	T	C	99.5%	610.3 / ‑5.7	196	F86L (TTT→CTT)	b1423	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (85/110); total (86/110)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 200 total aligned reads displayed.

AATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCATTGAT  >  NC_000913/1493477‑1493656

aaTGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGctt                                                                                         >  1:5625789/1‑93 (MQ=255)

  tGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGc                                                                                       >  1:5691630/1‑93 (MQ=255)

      gACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCAtt                                                                                   >  1:1755852/1‑93 (MQ=255)

      gACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCAtt                                                                                   >  2:1786396/1‑93 (MQ=255)

       aCGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTa                                                                                  <  1:5037245/93‑1 (MQ=255)

       aCGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTa                                                                                  <  1:3721617/93‑1 (MQ=255)

        cGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTAt                                                                                 >  2:4921072/1‑93 (MQ=255)

          aGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATgg                                                                               <  1:7695748/93‑1 (MQ=255)

          aGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATgg                                                                               >  1:4887501/1‑93 (MQ=255)

           gCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATggg                                                                              <  2:1544045/93‑1 (MQ=255)

             tGGCGCGACTGAATGTTGAGCGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTTCTCTTCGGCGGATGCTTGCCATTATGGGGa                                                                            <  2:4782624/93‑1 (MQ=255)

             tGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGa                                                                            <  1:7419408/93‑1 (MQ=255)

             tGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGa                                                                            <  2:4782626/93‑1 (MQ=255)

             tGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGa                                                                            >  2:1411502/1‑93 (MQ=255)

                 gcgACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTa                                                                        <  2:4489179/93‑1 (MQ=255)

                  cgACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTAc                                                                       <  1:6995967/93‑1 (MQ=255)

                  cgACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTAc                                                                       <  2:7886584/93‑1 (MQ=255)

                   gACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTAcc                                                                      >  2:6698830/1‑93 (MQ=255)

                   gACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTAcc                                                                      >  1:2895479/1‑93 (MQ=255)

                    aCTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGGGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTAccc                                                                     <  1:7638252/93‑1 (MQ=255)

                    aCTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTAccc                                                                     >  2:7768589/1‑93 (MQ=255)

                    aCTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTAccc                                                                     >  2:7419069/1‑93 (MQ=255)

                     cTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCg                                                                    >  1:4931215/1‑93 (MQ=255)

                      tGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCCTTATGGGGATGTACCCgg                                                                   >  2:2191697/1‑93 (MQ=255)

                      tGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCgg                                                                   >  1:1777036/1‑93 (MQ=255)

                        aaTGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTg                                                                 >  2:6069289/1‑93 (MQ=255)

                        aaTGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTg                                                                 <  2:1161013/93‑1 (MQ=255)

                         aTGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTgg                                                                >  1:3670580/1‑92 (MQ=255)

                         aTGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGa                                                                <  2:2562590/93‑1 (MQ=255)

                          tGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAg                                                               >  2:1743570/1‑93 (MQ=255)

                          tGTTGAACGTCAAGGGGCGATGCGCGGTGGGACTGCACAAGAGCTTGCTACTCTTCTGCGGATGCTGGCCATTATGGGGATGTACCCGGTGAg                                                               >  2:7019888/1‑93 (MQ=255)

                           gTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGc                                                              >  2:3954641/1‑93 (MQ=255)

                             tGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTa                                                            <  1:5441182/93‑1 (MQ=255)

                             tGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTa                                                            <  1:2319312/93‑1 (MQ=255)

                              gAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTAt                                                           >  1:5546706/1‑93 (MQ=255)

                              gAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTAt                                                           <  1:5109327/93‑1 (MQ=255)

                              gAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTAt                                                           <  2:2115761/93‑1 (MQ=255)

                               aaCGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTAtt                                                          <  1:3362391/93‑1 (MQ=255)

                                 cGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGTCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTAc                                                        <  1:206737/93‑1 (MQ=255)

                                 cGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTAc                                                        <  1:6819281/93‑1 (MQ=255)

                                 cGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTAc                                                        <  1:1081385/93‑1 (MQ=255)

                                  gTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACg                                                       <  2:7540067/93‑1 (MQ=255)

                                  gTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACg                                                       >  1:535557/1‑93 (MQ=255)

                                   tCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGa                                                      <  2:310941/93‑1 (MQ=255)

                                     ctGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGAtc                                                    <  2:3697795/92‑1 (MQ=255)

                                     aTGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGAtc                                                    <  1:1202508/93‑1 (MQ=255)

                                      tGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGAtct                                                   >  1:900202/1‑93 (MQ=255)

                                       ggggCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGAtctc                                                  <  1:7680875/93‑1 (MQ=255)

                                       ggggCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGAtctc                                                  >  1:6062901/1‑93 (MQ=255)

                                       ggggCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGAtctc                                                  <  2:6830355/93‑1 (MQ=255)

                                        gggCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGGCATTATGGGGATGTAACCGGTGAGCTATTACGAtctct                                                 >  2:6537207/1‑93 (MQ=255)

                                        gggCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGAtctct                                                 <  2:3611733/93‑1 (MQ=255)

                                         ggCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGAtctctc                                                <  2:1923635/93‑1 (MQ=255)

                                         ggCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGAtctctc                                                <  1:2542374/93‑1 (MQ=255)

                                          gCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCg                                               <  2:4697311/93‑1 (MQ=255)

                                            gATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCgca                                             <  2:2807739/93‑1 (MQ=255)

                                            gATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCgca                                             <  2:2682185/93‑1 (MQ=255)

                                             aTTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCgcag                                            >  2:4351739/1‑93 (MQ=255)

                                                cgcgTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCgcaggca                                         <  1:5378736/93‑1 (MQ=255)

                                                cgcgTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCgcaggca                                         <  2:3376187/93‑1 (MQ=255)

                                                 gcgTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCgcaggcag                                        >  1:6388839/1‑93 (MQ=255)

                                                   gTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAggg                                      <  1:4780418/93‑1 (MQ=255)

                                                    ttGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAgggg                                     >  2:4275913/1‑93 (MQ=255)

                                                     tGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGt                                    <  1:281652/93‑1 (MQ=255)

                                                     tGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGt                                    <  1:148923/93‑1 (MQ=255)

                                                         aCTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCg                                <  2:134656/93‑1 (MQ=255)

                                                         aCTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCg                                >  1:3247800/1‑93 (MQ=255)

                                                         aCTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCg                                <  2:5688688/93‑1 (MQ=255)

                                                           tGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGt                              <  1:2127681/93‑1 (MQ=255)

                                                            gCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTa                             <  1:4733409/93‑1 (MQ=255)

                                                              acaAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACa                           >  2:2635754/1‑93 (MQ=255)

                                                              acaAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACa                           <  1:7755593/93‑1 (MQ=255)

                                                                aaGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACAtt                         <  1:7079911/93‑1 (MQ=255)

                                                                 agagCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTc                        >  1:6191430/1‑93 (MQ=255)

                                                                  gagCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCg                       >  2:270725/1‑93 (MQ=255)

                                                                    gCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGAc                     >  1:6347125/1‑93 (MQ=255)

                                                                      ttGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAg                   >  2:2490861/1‑93 (MQ=255)

                                                                       tGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGc                  <  2:6064559/93‑1 (MQ=255)

                                                                       tGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGc                  >  2:938144/1‑93 (MQ=255)

                                                                        gCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCa                 <  1:3472113/93‑1 (MQ=255)

                                                                         cTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCAt                <  2:1726065/93‑1 (MQ=255)

                                                                         cTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCAt                <  1:3356139/93‑1 (MQ=255)

                                                                          tACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCAtt               <  2:1126373/93‑1 (MQ=255)

                                                                           tctctTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCAttt              <  2:5743432/92‑1 (MQ=255)

                                                                            ctctTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTc             <  1:6849848/93‑1 (MQ=255)

                                                                            ctctTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTc             >  2:838561/1‑93 (MQ=255)

                                                                            ctctTCGGCGGATGCTTGCCATTACGGGGATGTACCCGGTGAGCTATTACGACCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTc             <  1:3963859/93‑1 (MQ=255)

                                                                             tctTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCg            >  2:6354444/1‑93 (MQ=255)

                                                                               ttCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGc          >  2:3993611/1‑93 (MQ=255)

                                                                               ttCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGc          >  1:6438155/1‑93 (MQ=255)

                                                                               ttCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGc          >  1:970017/1‑93 (MQ=255)

                                                                                tCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGGTATTACGATCTCTCGCAGGCAGGGGTGCCGGTAGATTCGACAGCATTTGGGcc         >  2:5657924/1‑93 (MQ=255)

                                                                                tCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGcc         >  1:3173736/1‑93 (MQ=255)

                                                                                  ggcggATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCa       >  1:4524473/1‑93 (MQ=255)

                                                                                   gcggATGCTTGCCATTATGGAGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCAt      <  2:5154347/93‑1 (MQ=255)

                                                                                     ggATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCATtg    >  2:62679/1‑93 (MQ=255)

                                                                                     ggATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCATtg    <  1:4896129/93‑1 (MQ=255)

                                                                                      gATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCATtga   >  1:7850066/1‑93 (MQ=255)

                                                                                      gATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCATtga   >  2:7640676/1‑93 (MQ=255)

                                                                                      gATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCATtga   >  1:3149145/1‑93 (MQ=255)

AATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCATTGAT  >  NC_000913/1493477‑1493656

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