BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,915,434	G→A	A385T (GCT→ACT)	b1833 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,915,434	0	G	A	99.6%	844.0 / ‑6.0	244	A385T (GCT→ACT)	b1833	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (124/119); total (124/120)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

Only 100 of 249 total aligned reads displayed.

CGGATGGTGACCTGGATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTT  >  NC_000913/1915344‑1915523

cGGATGGTGACCTGGATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACt                                                                                         <  2:7244454/93‑1 (MQ=255)

  gATGGTGACCTGGATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACttt                                                                                       >  1:2182692/1‑93 (MQ=255)

    tGGTGACCTGGATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTa                                                                                     >  2:2191302/1‑93 (MQ=255)

       tGACCTGGATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTa                                                                                  <  2:3537592/93‑1 (MQ=255)

         aCCTGGATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATg                                                                                <  2:4658142/93‑1 (MQ=255)

          ccTGGATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATgg                                                                               <  2:3243971/93‑1 (MQ=255)

          ccTGGATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATgg                                                                               >  1:1815924/1‑93 (MQ=255)

           cTGGATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATggg                                                                              >  1:343481/1‑93 (MQ=255)

           cTGGATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATggg                                                                              >  2:4901022/1‑93 (MQ=255)

           cTGGATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATggg                                                                              >  2:6098290/1‑93 (MQ=255)

             ggATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGAc                                                                            <  1:1041590/93‑1 (MQ=255)

             ggATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGAc                                                                            <  1:6142051/93‑1 (MQ=255)

             ggATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGAc                                                                            >  2:2965924/1‑93 (MQ=255)

              gATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGAcc                                                                           <  1:5914107/93‑1 (MQ=255)

              gATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGAcc                                                                           >  2:4374198/1‑93 (MQ=255)

               aTTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCg                                                                          >  2:5853111/1‑93 (MQ=255)

               aTTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCg                                                                          >  2:6892025/1‑93 (MQ=255)

               aTTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCg                                                                          >  1:3544558/1‑93 (MQ=255)

                   gTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATACTCACTTTTACTATGGGACCGGCTg                                                                      >  1:6389971/1‑93 (MQ=255)

                    tCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGc                                                                     <  1:697603/93‑1 (MQ=255)

                       cTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGtt                                                                  <  1:7520724/93‑1 (MQ=255)

                       cTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGtt                                                                  >  2:7135801/1‑93 (MQ=255)

                        tGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGttt                                                                 >  2:458654/1‑93 (MQ=255)

                          gTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTa                                                               <  1:2678836/93‑1 (MQ=255)

                          gTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTa                                                               <  1:1487422/93‑1 (MQ=255)

                            cGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTAtt                                                             <  1:4551308/93‑1 (MQ=255)

                             gATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGatttattt                                                            >  1:2617893/1‑93 (MQ=255)

                             gATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTAttt                                                            <  2:5819815/93‑1 (MQ=255)

                              aTGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTc                                                           <  1:1977608/93‑1 (MQ=255)

                               tGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCg                                                          <  1:5668820/93‑1 (MQ=255)

                                gCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCgg                                                         >  2:6011357/1‑93 (MQ=255)

                                gCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCgg                                                         <  1:1339215/93‑1 (MQ=255)

                                 cTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGgc                                                        <  2:6372283/93‑1 (MQ=255)

                                    ggccTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCa                                                     <  2:7575932/91‑1 (MQ=255)

                                    gACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCa                                                     <  2:5021738/93‑1 (MQ=255)

                                    gACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCa                                                     <  1:5511470/93‑1 (MQ=255)

                                    gACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCa                                                     <  1:3993828/93‑1 (MQ=255)

                                      ccTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGc                                                   <  1:4336542/93‑1 (MQ=255)

                                      ccTGTTTGTCATATCTTTAACCATGTCGCTGATGAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGc                                                   <  1:7639401/93‑1 (MQ=255)

                                       cTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCg                                                  >  1:3337833/1‑93 (MQ=255)

                                       cTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCg                                                  >  2:1543516/1‑93 (MQ=255)

                                         gTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGt                                                >  2:4104076/1‑93 (MQ=255)

                                         gTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGt                                                >  2:3636692/1‑93 (MQ=255)

                                         gTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGt                                                <  2:3165308/93‑1 (MQ=255)

                                         gTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGt                                                <  1:1932376/93‑1 (MQ=255)

                                           ttGTCATATCTTTAACTATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTGTATGTCTGCGCATCGGTaa                                              >  2:7033052/1‑93 (MQ=255)

                                           ttGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTaa                                              >  1:5232958/1‑93 (MQ=255)

                                           ttGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTaa                                              <  2:891073/93‑1 (MQ=255)

                                            tGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAAt                                             >  2:7898362/1‑93 (MQ=255)

                                              tCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAAttt                                           <  1:3939080/93‑1 (MQ=255)

                                              tCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAAttt                                           <  1:4823332/93‑1 (MQ=255)

                                               cATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAAtttt                                          >  1:1544431/1‑93 (MQ=255)

                                               cATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAAtttt                                          >  2:4806919/1‑93 (MQ=255)

                                                atatCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGAGCAGCGGTAATTTTg                                         >  2:6816022/1‑93 (MQ=255)

                                                 tatCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGa                                        <  2:4439190/93‑1 (MQ=255)

                                                 tatCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGa                                        <  1:7174797/93‑1 (MQ=255)

                                                  atCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGAc                                       >  1:5757835/1‑93 (MQ=255)

                                                     tttAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTAt                                    >  1:2452021/1‑93 (MQ=255)

                                                     tttAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTAt                                    >  2:1623683/1‑93 (MQ=255)

                                                     tttAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTAt                                    <  1:7742448/93‑1 (MQ=255)

                                                      ttAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTAtt                                   >  2:7715920/1‑93 (MQ=255)

                                                       tAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTc                                  <  2:1204708/93‑1 (MQ=255)

                                                       tAACCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTc                                  <  1:6961627/93‑1 (MQ=255)

                                                        aaCCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCt                                 <  1:6502858/93‑1 (MQ=255)

                                                         aCCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCtt                                <  2:5674723/93‑1 (MQ=255)

                                                         aCCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCtt                                >  2:1171216/1‑93 (MQ=255)

                                                         aCCATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCtt                                >  2:3921510/1‑93 (MQ=255)

                                                           cATGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGc                              >  1:5984474/1‑93 (MQ=255)

                                                            aTGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCt                             >  1:1203946/1‑93 (MQ=255)

                                                            aTGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCt                             >  2:1443775/1‑93 (MQ=255)

                                                            aTGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCt                             >  2:2564062/1‑93 (MQ=255)

                                                            aTGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGGAGTTGTGACTTTTCCTGCt                             >  2:5062957/1‑93 (MQ=255)

                                                             tGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTCTGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCtg                            >  1:5034968/1‑93 (MQ=255)

                                                             tGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCtg                            >  2:6274417/1‑93 (MQ=255)

                                                             tGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCtg                            >  1:5088137/1‑93 (MQ=255)

                                                             tGTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCtg                            >  1:6856707/1‑93 (MQ=255)

                                                              gTCGCTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCtgt                           <  2:149847/93‑1 (MQ=255)

                                                                  cTGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGaa                       <  1:7608379/93‑1 (MQ=255)

                                                                   tGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAAt                      <  2:7304461/93‑1 (MQ=255)

                                                                   tGATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAAt                      >  1:2896871/1‑93 (MQ=255)

                                                                    gATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATg                     >  2:7196769/1‑93 (MQ=255)

                                                                    gATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATg                     <  2:5367846/93‑1 (MQ=255)

                                                                    gATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATg                     >  2:3113559/1‑93 (MQ=255)

                                                                    gATTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATg                     <  1:2962093/93‑1 (MQ=255)

                                                                     aTTAATCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATgg                    <  2:4084900/93‑1 (MQ=255)

                                                                        aaTCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTg                 <  2:551852/93‑1 (MQ=255)

                                                                          tCGCGATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGa               <  2:1815924/93‑1 (MQ=255)

                                                                           cgcgATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGAc              <  2:2122516/93‑1 (MQ=255)

                                                                           cgcgATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGAc              <  2:3848463/93‑1 (MQ=255)

                                                                              gATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCCGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTAGACAga           >  2:4842922/1‑92 (MQ=255)

                                                                              gATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGc           >  1:874675/1‑93 (MQ=255)

                                                                              gATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGc           >  1:3246416/1‑93 (MQ=255)

                                                                              gATCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGc           <  1:7246816/93‑1 (MQ=255)

                                                                               aTCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGcc          <  2:34782/93‑1 (MQ=255)

                                                                                tCAGATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCg         >  2:2659969/1‑93 (MQ=255)

                                                                                   gATCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCtt      >  1:4679724/1‑93 (MQ=255)

                                                                                    aTCCTCGCTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTa     <  1:7429646/93‑1 (MQ=255)

                                                                                    aTCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTa     <  2:874288/93‑1 (MQ=255)

                                                                                    aTCCTCACTTTTACTATGGGACCGGCTGCGGTTTATTTTGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTa     >  2:3580715/1‑93 (MQ=255)

                                                                                     tCCTCACTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTAc    >  1:2083954/1‑93 (MQ=255)

CGGATGGTGACCTGGATTGGTCGCTGGTCGATGCTCGACCTGTTTGTCATATCTTTAACCATGTCGCTGATTAATCGCGATCAGATCCTCGCTTTTACTATGGGACCGGCTGCGTTTTATTTCGGCGCAGCGGTAATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTT  >  NC_000913/1915344‑1915523

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