BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,273,643	G→A	P198P (CCG→CCA)	yjcC →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,273,643	0	G	A	97.2%	794.8 / 4.2	217	P198P (CCG→CCA)	yjcC	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (2/4); new base A (106/105); total (108/109)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

Only 100 of 222 total aligned reads displayed.

TGACCCGACATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGT  >  NC_000913/4273553‑4273732

tgaCCCGACATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACt                                                                                         >  1:1808722/1‑93 (MQ=255)

tgaCCCGACATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACt                                                                                         <  1:2352649/93‑1 (MQ=255)

 gaCCCGACATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTg                                                                                        >  1:401529/1‑93 (MQ=255)

 gaCCCGACATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTg                                                                                        >  1:2434956/1‑93 (MQ=255)

    ccGACATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGAtt                                                                                     >  2:5214671/1‑93 (MQ=255)

    ccGACATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCCCCACTGAtt                                                                                     <  1:2141831/93‑1 (MQ=255)

      gACATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCa                                                                                   >  2:4397479/1‑93 (MQ=255)

       aCATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGAtttct                                                                                  >  1:74265/1‑90 (MQ=255)

       aCATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCAt                                                                                  <  1:3864620/93‑1 (MQ=255)

        cATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCAtt                                                                                 <  2:263446/93‑1 (MQ=255)

        cATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCAtt                                                                                 <  2:534180/93‑1 (MQ=255)

           tGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGa                                                                              >  2:5226453/1‑93 (MQ=255)

           tGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGa                                                                              <  2:614771/93‑1 (MQ=255)

             cAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAg                                                                            <  2:724179/93‑1 (MQ=255)

             cAATGGGGTGGGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAg                                                                            <  1:663442/93‑1 (MQ=255)

              aaTGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAgg                                                                           <  1:1701275/93‑1 (MQ=255)

               aTGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGa                                                                          >  1:1400836/1‑93 (MQ=255)

                  gggTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGAttt                                                                       <  2:3521354/93‑1 (MQ=255)

                   ggTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTa                                                                      <  2:3392622/93‑1 (MQ=255)

                    gtgtgtATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTaa                                                                     <  1:5692882/93‑1 (MQ=255)

                    gtgtgtATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTaa                                                                     <  2:3465344/93‑1 (MQ=255)

                     tgtgtATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAAc                                                                    <  1:3483189/93‑1 (MQ=255)

                      gtgtATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAAcc                                                                   >  1:2389259/1‑93 (MQ=255)

                      gtgtATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAAcc                                                                   <  2:3190328/93‑1 (MQ=255)

                      gtgtATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAAcc                                                                   <  1:5633617/93‑1 (MQ=255)

                       tgtATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCg                                                                  >  2:2573994/1‑93 (MQ=255)

                       tgtATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCg                                                                  >  2:620516/1‑93 (MQ=255)

                        gtATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGt                                                                 <  1:4798531/93‑1 (MQ=255)

                          aTGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCTGGAGCAACGTTTTCTCCACTGATTCATTTGAATGATTTAACCGTat                                                               >  1:592736/1‑92 (MQ=255)

                          aTGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTAc                                                               >  1:4093969/1‑93 (MQ=255)

                           tGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACa                                                              <  2:5159396/93‑1 (MQ=255)

                           tGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACa                                                              <  1:2834680/93‑1 (MQ=255)

                            gATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACaa                                                             >  1:2515891/1‑93 (MQ=255)

                              tACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACaaaa                                                           >  1:3557903/1‑92 (MQ=255)

                              tACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAg                                                           >  1:535996/1‑93 (MQ=255)

                              tACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAg                                                           >  2:1449882/1‑93 (MQ=255)

                               aCGGTGACGAATACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGa                                                          <  2:670171/93‑1 (MQ=255)

                               aCGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGa                                                          >  2:5731680/1‑93 (MQ=255)

                                 ggTGACGAAAACCTTTTTCTCGTTGAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGaaa                                                        <  2:5488307/93‑1 (MQ=255)

                                 ggTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGaaa                                                        >  1:5863971/1‑93 (MQ=255)

                                   tGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATg                                                      <  2:4156129/93‑1 (MQ=255)

                                   tGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATg                                                      >  2:4491100/1‑93 (MQ=255)

                                    gACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATgg                                                     <  2:1197298/93‑1 (MQ=255)

                                     aCGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGc                                                    >  1:3169119/1‑93 (MQ=255)

                                        aaaaCCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTAt                                                 >  1:1930861/1‑93 (MQ=255)

                                         aaaCCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTTAAGGATTTAACCGTACAAAGAAATGGCTAtt                                                >  2:1041418/1‑93 (MQ=255)

                                         aaaCCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTAtt                                                <  2:2750582/93‑1 (MQ=255)

                                          aaCCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTAttt                                               <  2:2043089/93‑1 (MQ=255)

                                           aCCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTta                                              <  1:4694910/93‑1 (MQ=255)

                                            ccTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTtat                                             <  2:6023631/93‑1 (MQ=255)

                                            ccTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTtat                                             >  1:685421/1‑93 (MQ=255)

                                            ccTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTtat                                             >  1:6092486/1‑93 (MQ=255)

                                            ccTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTACACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTtat                                             >  2:2263163/1‑93 (MQ=255)

                                             cTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTtata                                            >  1:108088/1‑93 (MQ=255)

                                             cTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTtata                                            <  1:4297890/93‑1 (MQ=255)

                                             cTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTtata                                            <  1:901820/93‑1 (MQ=255)

                                              tttttCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTatat                                           <  2:1871349/93‑1 (MQ=255)

                                              ttgttCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTatat                                           <  1:2662365/90‑1 (MQ=255)

                                                tttCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGc                                         <  2:1663425/93‑1 (MQ=255)

                                                tttCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGc                                         <  2:1225035/93‑1 (MQ=255)

                                                tttCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGc                                         >  2:5521259/1‑93 (MQ=255)

                                                 ttCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCg                                        >  1:1834267/1‑93 (MQ=255)

                                                 ttCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCg                                        <  2:3665618/93‑1 (MQ=255)

                                                   ctcGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGAc                                      <  2:3023536/93‑1 (MQ=255)

                                                    tcGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCCTTTGAAGGCTTTAACCCTACAAAGAAATGGCTATTTAAATGCGACa                                     >  2:4712367/1‑93 (MQ=255)

                                                    tcGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACa                                     >  1:4148527/1‑93 (MQ=255)

                                                    tcGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACa                                     >  1:6047621/1‑93 (MQ=255)

                                                     cGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAg                                    <  1:1989463/93‑1 (MQ=255)

                                                      gTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGt                                   <  2:371345/93‑1 (MQ=255)

                                                          aGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTAt                               >  1:2776535/1‑93 (MQ=255)

                                                           acAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTAtt                              <  1:1450530/92‑1 (MQ=255)

                                                            cAAAGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTc                             <  2:453506/93‑1 (MQ=255)

                                                             aaaGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCg                            >  1:1045220/1‑93 (MQ=255)

                                                             aaaGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCg                            >  2:4184624/1‑93 (MQ=255)

                                                              aaGAGGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGa                           >  2:1350540/1‑93 (MQ=255)

                                                               agagGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGAc                          <  1:497165/93‑1 (MQ=255)

                                                               agagGCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGAc                          >  2:2375380/1‑93 (MQ=255)

                                                                   gCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACaaaa                      <  2:3766994/93‑1 (MQ=255)

                                                                   gCCTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACaaaa                      <  1:1612230/93‑1 (MQ=255)

                                                                    ccTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAAc                     <  2:4426116/93‑1 (MQ=255)

                                                                    ccTCGGCAGCAACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAAc                     <  2:3629038/93‑1 (MQ=255)

                                                                         gcagcaACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCa                >  2:512565/1‑93 (MQ=255)

                                                                         gcagcaACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCa                >  2:5801915/1‑93 (MQ=255)

                                                                         gcagcaACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCa                >  2:6126620/1‑93 (MQ=255)

                                                                         gcagcaACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCa                >  1:3436251/1‑93 (MQ=255)

                                                                         gcagcaACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCa                >  2:4812890/1‑93 (MQ=255)

                                                                          cgtcaACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCaa               <  2:5895406/90‑1 (MQ=255)

                                                                           agcaACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAAt              >  2:4341885/1‑93 (MQ=255)

                                                                            gcaACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAAtt             >  2:307657/1‑93 (MQ=255)

                                                                             caACGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTg            >  2:5005330/1‑93 (MQ=255)

                                                                              aaCGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGc           >  1:5278529/1‑93 (MQ=255)

                                                                                cGTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAg         >  2:1594778/1‑93 (MQ=255)

                                                                                 gTTTTCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGc        >  1:97980/1‑93 (MQ=255)

                                                                                  ttttCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGATTATTCGACAAAACGCCCAATTGCAGcc       >  2:1109731/1‑93 (MQ=255)

                                                                                   tttCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATACGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCa      <  1:6116815/93‑1 (MQ=255)

                                                                                    ttCTCCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCAt     >  2:1503641/1‑93 (MQ=255)

                                                                                      ctcCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCAttg   >  1:2005511/1‑93 (MQ=255)

                                                                                       tcCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCAttgt  >  2:2697566/1‑93 (MQ=255)

                                                                                       tcCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCAttgt  <  2:5331951/93‑1 (MQ=255)

                                                                                       tcCACTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCAttgt  >  1:340307/1‑93 (MQ=255)

TGACCCGACATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGT  >  NC_000913/4273553‑4273732

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