BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,488,636	T→C	E445E (GAA→GAG)	b2373 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,488,636	0	T	C	87.9%	529.9 / 56.9	191	E445E (GAA→GAG)	b2373	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (15/8); new base C (125/42); total (141/50)
Fisher's exact test for biased strand distribution p-value = 3.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.30e-01

Only 100 of 196 total aligned reads displayed.

CGCCATTATTAAAAATAACGATCGTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCTTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGCCCATGACACCCCAGGTACCACAATCCA  >  NC_000913/2488545‑2488727

 gccATTATTAAAAATAACGATCGTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTc                                                                                           >  1:5609964/1‑93 (MQ=255)

  ccATTATTAAAAATAACGATCGTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCa                                                                                          >  1:5021845/1‑93 (MQ=255)

  ccATTATTAAAAATAACGATCGTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCa                                                                                          >  2:3618028/1‑93 (MQ=255)

    attattAAAAATAACGATCGTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAAt                                                                                        >  1:4086190/1‑93 (MQ=255)

    attattAAAAATAACGATCGTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAAt                                                                                        >  2:108575/1‑93 (MQ=255)

     ttattaAAAATAACGATCGTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATg                                                                                       >  1:6338491/1‑93 (MQ=255)

      tattaAAAATAACGATCGTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATgg                                                                                      >  1:4971273/1‑93 (MQ=255)

      tattaAAAATAACGATCGTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATgg                                                                                      >  1:577442/1‑93 (MQ=255)

       attaAAAATAACGATCGTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGc                                                                                     >  1:232621/1‑93 (MQ=255)

       attaAAAATAACGATCGTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGc                                                                                     <  1:2815561/93‑1 (MQ=255)

       attaAAAATAACGATCGTCACCGGCAGGTTATATCGACAAATCGGTTCAATTTCAATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGc                                                                                     >  2:3357123/1‑93 (MQ=255)

        ttaAAAATAACGATCGTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCg                                                                                    >  1:584493/1‑93 (MQ=255)

            aaaTAACGATCGTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACaa                                                                                >  1:5140649/1‑93 (MQ=255)

            aaaTAACGATCGTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACaa                                                                                <  2:2856714/93‑1 (MQ=255)

             aaTAACGATCGTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAAc                                                                               <  1:2453561/93‑1 (MQ=255)

             aaTAACGATCGTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAAc                                                                               <  2:1735985/93‑1 (MQ=255)

              aTAACGATCGTCGCCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCGCCCTCAATGGCGACAAcc                                                                              <  2:1318014/93‑1 (MQ=255)

                  cGATCGTCAGCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCTTCAATGGCGACAACCGgag                                                                          <  2:2689657/93‑1 (MQ=255)

                  cGATCGTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGgag                                                                          >  1:1915305/1‑93 (MQ=255)

                     tCGTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGTACTATCACCCTCAATGGAGACAACCGGAGAAc                                                                       >  2:4554336/1‑93 (MQ=255)

                     tCGTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCTTCAATGGCGACAACCGGAGAAc                                                                       >  2:3958464/1‑93 (MQ=255)

                     tCGTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAAc                                                                       >  1:3053/1‑93 (MQ=255)

                     tCGTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAAc                                                                       >  1:1008030/1‑93 (MQ=255)

                      cGTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAAcc                                                                      >  2:5149360/1‑93 (MQ=255)

                      cGTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAAcc                                                                      >  2:3451164/1‑93 (MQ=255)

                      cGTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAAcc                                                                      >  2:5324555/1‑93 (MQ=255)

                      cGTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATCGCGACAACCGGAGAAcc                                                                      >  2:1600944/1‑93 (MQ=255)

                       gTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCa                                                                     >  2:4208845/1‑93 (MQ=255)

                       gTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCa                                                                     <  1:3145153/93‑1 (MQ=255)

                       gTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCa                                                                     <  2:6011896/93‑1 (MQ=255)

                       gTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCa                                                                     >  1:1373597/1‑93 (MQ=255)

                        tCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAg                                                                    >  2:2806268/1‑93 (MQ=255)

                        tCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAg                                                                    >  2:375426/1‑93 (MQ=255)

                        tCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAg                                                                    >  1:5072369/1‑93 (MQ=255)

                          aCCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGaa                                                                  >  2:1145658/1‑93 (MQ=255)

                          aCCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGaa                                                                  <  1:1760803/93‑1 (MQ=255)

                           ccGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAg                                                                 >  1:4323588/1‑93 (MQ=255)

                             ggCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTc                                                               <  1:3731091/93‑1 (MQ=255)

                              gCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCa                                                              >  2:1212139/1‑93 (MQ=255)

                                 ggTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGc                                                           >  2:233720/1‑93 (MQ=255)

                                 ggTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGc                                                           >  1:486901/1‑93 (MQ=255)

                                 ggTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGc                                                           <  1:2281273/93‑1 (MQ=255)

                                 ggTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGc                                                           >  2:4970008/1‑93 (MQ=255)

                                 ggTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGc                                                           >  1:3608028/1‑93 (MQ=255)

                                  gTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCt                                                          <  1:2136541/93‑1 (MQ=255)

                                    tataTCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAg                                                        >  2:3507685/1‑93 (MQ=255)

                                     atatCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGc                                                       >  2:4381028/1‑93 (MQ=255)

                                     atatCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGc                                                       >  1:4760448/1‑93 (MQ=255)

                                      tatCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCa                                                      >  1:4216466/1‑93 (MQ=255)

                                      tatCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCa                                                      >  2:5574164/1‑93 (MQ=255)

                                       atCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCAc                                                     >  2:2495113/1‑93 (MQ=255)

                                       atCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCAc                                                     <  2:5601173/93‑1 (MQ=255)

                                          gACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCCCCTTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGa                                                  <  1:5351866/93‑1 (MQ=255)

                                          gACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCTTCAATGGCGACAAACAGAGAACCAGAAGTCACGCTAGCACCGa                                                  >  1:4337104/1‑93 (MQ=255)

                                          gACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGa                                                  <  2:3533539/93‑1 (MQ=255)

                                           aCAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGAt                                                 <  2:4611543/93‑1 (MQ=255)

                                           aCAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGAt                                                 <  2:844549/93‑1 (MQ=255)

                                            cAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATg                                                <  2:3232498/93‑1 (MQ=255)

                                              aaTCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCGCCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGc                                              >  2:3553532/1‑93 (MQ=255)

                                              aaTCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGc                                              >  1:2549361/1‑93 (MQ=255)

                                              aaTCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGc                                              >  1:2962658/1‑93 (MQ=255)

                                               aTCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCa                                             >  2:5953849/1‑93 (MQ=255)

                                                 cGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATa                                           <  1:2576709/93‑1 (MQ=255)

                                                    ttCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGcc                                        >  1:1606676/1‑93 (MQ=255)

                                                       aaTTTCCATCCCACTGAAACCAAAAGCACTATCACCTTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCAt                                     >  2:3028471/1‑93 (MQ=255)

                                                        aTTTCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATa                                    >  1:1086126/1‑93 (MQ=255)

                                                           tCCATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCg                                 >  2:3121010/1‑93 (MQ=255)

                                                            ccATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGa                                >  1:6410523/1‑93 (MQ=255)

                                                             cATCCCACTGAAACCAAAAGCACTATCACCTTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGAt                               >  2:2594560/1‑93 (MQ=255)

                                                             cATCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGAt                               >  2:3172738/1‑93 (MQ=255)

                                                               tCCCACTGAAACCAAAAGCACTATCACCTTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGc                             >  1:4596220/1‑93 (MQ=255)

                                                               tCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGc                             >  2:2211521/1‑93 (MQ=255)

                                                               tCCCACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGc                             >  2:4697059/1‑93 (MQ=255)

                                                                cccACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGcc                            >  2:1113670/1‑93 (MQ=255)

                                                                cccACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGcc                            >  2:507635/1‑93 (MQ=255)

                                                                 ccACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGccc                           >  2:6618129/1‑93 (MQ=255)

                                                                  cACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGCCCa                          >  1:3220700/1‑93 (MQ=255)

                                                                  cACTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGCCCa                          >  1:6307610/1‑93 (MQ=255)

                                                                   aCTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGCCCAt                         >  1:3491525/1‑93 (MQ=255)

                                                                   aCTGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGCCCAt                         >  1:3006893/1‑93 (MQ=255)

                                                                     tGAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATACCCCATACCGATGCCCATGa                       >  2:4331387/1‑93 (MQ=255)

                                                                      gAAACCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGCCCATGac                      >  1:6596345/1‑93 (MQ=255)

                                                                         aCCAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGCCCATGACAcc                   <  1:5461237/93‑1 (MQ=255)

                                                                          ccAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGCCCATGACAccc                  >  2:3783039/1‑93 (MQ=255)

                                                                          ccAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGCCCATGACAccc                  <  1:2627343/93‑1 (MQ=255)

                                                                          ccAAAAGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGCCCAGGACAccc                  <  1:4240304/93‑1 (MQ=255)

                                                                            aaaaGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGCCCATGACACCCCa                <  1:5986813/93‑1 (MQ=255)

                                                                               aGCACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGCCCATGACACCCCAGGt             >  2:5294662/1‑93 (MQ=255)

                                                                                gCACTATCACCTTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGCCCATGACACCCCAGGTa            >  2:5235737/1‑93 (MQ=255)

                                                                                gCACTATCACCATCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGCCCATGACACCCCAGGTa            >  1:1667871/1‑93 (MQ=255)

                                                                                 cACTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGCCCATGACACCCCAGGTAc           >  1:6140382/1‑93 (MQ=255)

                                                                                  aCTATCACCTTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGCCCATGACACCCCAGGTAcc          >  2:276757/1‑93 (MQ=255)

                                                                                   cTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGCCCATGACACCCCAGGTACca         >  2:57114/1‑93 (MQ=255)

                                                                                   cTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGCCCATGACACCCCAGGTACca         >  2:4812122/1‑93 (MQ=255)

                                                                                   cTATCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGCCCATGACACCCCAGGTACca         >  1:3748856/1‑93 (MQ=255)

                                                                                    tATCACCTTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGCCCATGACACCCCAGGTACcac        >  2:6283065/1‑93 (MQ=255)

                                                                                      tCACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGCCCATGACACCCCAGGTACCACaa      >  2:2016512/1‑93 (MQ=255)

                                                                                       cACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGCCCATGACACCCCAGGTACCACAAt     <  2:4086190/93‑1 (MQ=255)

                                                                                       cACCCTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGCCCATGACACCCCAGGTACCACAAt     >  2:1487224/1‑93 (MQ=255)

                                                                                          cTTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGCCCATGACACCCCAGGTACCACAATCCa  <  1:2174169/93‑1 (MQ=255)

CGCCATTATTAAAAATAACGATCGTCACCGGCAGGTTATATCGACAAATCGTTTCAATTTCCATCCCACTGAAACCAAAAGCACTATCACCTTCAATGGCGACAACCGGAGAACCAGAAGTCACGCTAGCACCGATGGCATAGCCCATACCGATGCCCATGACACCCCAGGTACCACAATCCA  >  NC_000913/2488545‑2488727

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