BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,486,394	T→C	A278A (GCA→GCG)	yfdE ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,486,394	0	T	C	90.3%	753.2 / 59.7	279	A278A (GCA→GCG)	yfdE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (13/14); new base C (136/115); total (149/130)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 280 total aligned reads displayed.

CGCTCAATATATTGTTTAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGTGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATG  >  NC_000913/2486302‑2486484

cGCTCAATATATTGTTTAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGt                                                                                             <  1:6753205/93‑1 (MQ=255)

   tCAATATATTGTTTAAGAATAGCCTGGTTTTTTACGCGGAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGcc                                                                                          >  1:6269481/1‑93 (MQ=255)

   tCAATATATTGTTTAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGcc                                                                                          >  1:3913820/1‑93 (MQ=255)

   tCAATATATTGTTTAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGcc                                                                                          <  2:1314933/93‑1 (MQ=255)

    cAATATATTGTTTAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGTGCCt                                                                                         <  1:1000234/93‑1 (MQ=255)

     aaTATATTGTTTAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTg                                                                                        >  1:7393045/1‑93 (MQ=255)

     aaTATATTGTTTAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTg                                                                                        <  2:6217288/93‑1 (MQ=255)

      atatatTGTTTAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTgg                                                                                       >  2:5354463/1‑93 (MQ=255)

        atatTGTTTAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCa                                                                                     >  1:3722685/1‑93 (MQ=255)

         tatTGTTTAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCAt                                                                                    >  1:3637011/1‑93 (MQ=255)

          atTGTTTAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATa                                                                                   >  1:2862528/1‑93 (MQ=255)

            tGTTTAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAAc                                                                                 <  2:6471047/93‑1 (MQ=255)

              tttAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGc                                                                               <  2:2076734/93‑1 (MQ=255)

               ttAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCa                                                                              <  1:5166067/93‑1 (MQ=255)

                tAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAg                                                                             >  1:656699/1‑93 (MQ=255)

                 aaGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑a                                                                           <  2:435520/93‑1 (MQ=255)

                 aaGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑a                                                                           <  1:3309960/93‑1 (MQ=255)

                   gAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑aaa                                                                         <  2:3038915/93‑1 (MQ=255)

                     aTAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑aaaaa                                                                       >  2:1006677/1‑93 (MQ=255)

                      tAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑aaaaaa                                                                      >  1:2681201/1‑93 (MQ=255)

                      tAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑aaaaaa                                                                      <  1:4236067/93‑1 (MQ=255)

                       aGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAg                                                                     <  1:5479415/93‑1 (MQ=255)

                        gCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGc                                                                    >  2:7030392/1‑93 (MQ=255)

                        gCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGc                                                                    >  1:4705420/1‑93 (MQ=255)

                         ccTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCt                                                                   >  2:2283521/1‑93 (MQ=255)

                          cTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCtt                                                                  >  2:7333586/1‑93 (MQ=255)

                          cTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCtt                                                                  <  2:1707789/93‑1 (MQ=255)

                          cTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCtt                                                                  >  1:3077775/1‑93 (MQ=255)

                           tGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTg                                                                 >  1:6500780/1‑93 (MQ=255)

                            ggTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGt                                                                >  1:1193008/1‑93 (MQ=255)

                            ggTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGt                                                                >  1:5322830/1‑93 (MQ=255)

                            ggTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGt                                                                <  2:2828186/93‑1 (MQ=255)

                             gTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTc                                                               >  1:4532330/1‑93 (MQ=255)

                             gTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTc                                                               >  1:6027898/1‑93 (MQ=255)

                              ttttGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCa                                                              <  2:6702957/93‑1 (MQ=255)

                               tttGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGTGCCTGGCATAACGCAG‑AAAAAAGCTTGTCAt                                                             >  2:4235147/1‑93 (MQ=255)

                               tttGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCAt                                                             <  1:5583168/93‑1 (MQ=255)

                                ttGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCAtt                                                            <  1:3140024/93‑1 (MQ=255)

                                 tGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTa                                                           <  1:7481097/93‑1 (MQ=255)

                                   tACGCGTAAAATATTGCTGCTAAATCAGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTAcc                                                         >  2:555814/1‑93 (MQ=255)

                                     cgcgTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACcac                                                       <  1:4994243/93‑1 (MQ=255)

                                     cgcgTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACcac                                                       >  2:3069899/1‑93 (MQ=255)

                                      gcgTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACcaca                                                      <  1:1864748/93‑1 (MQ=255)

                                      gcgTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACcaca                                                      <  2:3753662/93‑1 (MQ=255)

                                       cgTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCacaa                                                     >  2:7221478/1‑93 (MQ=255)

                                       cgTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCacaa                                                     <  1:3455775/93‑1 (MQ=255)

                                        gTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCacaac                                                    >  1:4225594/1‑93 (MQ=255)

                                           aaaTATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACaaa                                                 <  1:790369/93‑1 (MQ=255)

                                            aaTATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAAt                                                <  1:5375315/93‑1 (MQ=255)

                                            aaTATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAAt                                                <  1:2953781/93‑1 (MQ=255)

                                            aaTATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAAt                                                <  1:1878072/93‑1 (MQ=255)

                                            aaTATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAAt                                                <  2:4426141/93‑1 (MQ=255)

                                             atatTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATc                                               <  1:591138/93‑1 (MQ=255)

                                             atatTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATc                                               >  2:6073386/1‑93 (MQ=255)

                                              tatTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCg                                              >  1:491479/1‑93 (MQ=255)

                                              tatTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCg                                              >  1:1832213/1‑93 (MQ=255)

                                                 tgctgcTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTaa                                           <  2:1816769/93‑1 (MQ=255)

                                                 tgctgcTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTaa                                           <  2:5009070/93‑1 (MQ=255)

                                                   ctgctAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATc                                         <  1:1500203/93‑1 (MQ=255)

                                                    tgctAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCg                                        >  1:6671078/1‑93 (MQ=255)

                                                    tgctAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCg                                        <  1:5040938/93‑1 (MQ=255)

                                                    tgctAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCg                                        >  2:3415506/1‑93 (MQ=255)

                                                      ctAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGc                                      >  2:1447179/1‑93 (MQ=255)

                                                       tAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGTGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCt                                     >  2:7424432/1‑93 (MQ=255)

                                                       tAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCt                                     <  2:4017844/93‑1 (MQ=255)

                                                       tAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCt                                     <  2:379987/93‑1 (MQ=255)

                                                       tAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCt                                     >  1:4048055/1‑93 (MQ=255)

                                                        aaaTCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCtt                                    <  2:2571473/93‑1 (MQ=255)

                                                         aaTCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTa                                   <  1:1702248/93‑1 (MQ=255)

                                                         aaTCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTa                                   <  1:6834199/93‑1 (MQ=255)

                                                         aaTCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTa                                   >  2:5629470/1‑93 (MQ=255)

                                                             ggggATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATcca                               >  2:2923590/1‑92 (MQ=255)

                                                             ggggATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCt                               <  2:809265/93‑1 (MQ=255)

                                                              gggATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTg                              >  1:4860636/1‑93 (MQ=255)

                                                               ggATCATTAACCAGTTCCGTAAGCTCCAGTGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGa                             <  1:7097277/93‑1 (MQ=255)

                                                               ggATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGa                             <  2:5866966/93‑1 (MQ=255)

                                                                gATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAg                            >  1:3374602/1‑93 (MQ=255)

                                                                  tCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTa                          >  1:7524950/1‑93 (MQ=255)

                                                                    aTTAACCAGTTCCGTAAGCTCCAGTGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTAtt                        >  2:2171779/1‑93 (MQ=255)

                                                                    aTTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTAtt                        >  2:6817383/1‑93 (MQ=255)

                                                                        aCCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGaaa                    <  2:869407/93‑1 (MQ=255)

                                                                        aCCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGaaa                    >  1:6924139/1‑93 (MQ=255)

                                                                         ccAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGaaaa                   <  2:4096094/93‑1 (MQ=255)

                                                                         ccAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGaaaa                   <  1:5791984/93‑1 (MQ=255)

                                                                          cAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATTGTAATCGGCTTATCCTGAGTATTGAAAAc                  >  2:1147890/1‑93 (MQ=255)

                                                                           aGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACa                 <  1:2749015/93‑1 (MQ=255)

                                                                            gTTCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACAt                <  1:5923353/93‑1 (MQ=255)

                                                                             ttCCGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATc               >  2:4679802/1‑93 (MQ=255)

                                                                               ccGTAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCaa             <  1:3168432/93‑1 (MQ=255)

                                                                               ccGTAAGCTCCAGCGCCGGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCaa             >  2:4900079/1‑93 (MQ=255)

                                                                                  tAAGCTCCAGTGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAg          >  1:2343861/1‑93 (MQ=255)

                                                                                  tAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAg          <  1:22011/93‑1 (MQ=255)

                                                                                  tAAGCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAg          >  1:6151444/1‑93 (MQ=255)

                                                                                     gCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTg       >  1:103995/1‑93 (MQ=255)

                                                                                     gCTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTg       >  1:4338831/1‑93 (MQ=255)

                                                                                      cTCCAGTGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGc      <  2:23476/93‑1 (MQ=255)

                                                                                      cTCCAGCGCCTGGCATAACGCAGAAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATAAAAAGGTg       >  1:7313727/1‑93 (MQ=255)

                                                                                      cTCCAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGc      <  2:5708470/93‑1 (MQ=255)

                                                                                        ccAGCGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCa    <  1:3771493/93‑1 (MQ=255)

                                                                                          aGTGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATg  >  2:5232931/1‑93 (MQ=255)

CGCTCAATATATTGTTTAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGTGCCTGGCATAACGCAG‑AAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATG  >  NC_000913/2486302‑2486484

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