BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	709,880	T→C	intergenic (‑11/+278)	fur ← / ← fldA	Fur/flavodoxin 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	709,880	0	T	C	99.5%	661.2 / ‑5.1	202	intergenic (‑11/+278)	fur/fldA	Fur/flavodoxin 1
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (119/82); total (119/83)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 207 total aligned reads displayed.

TGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAA  >  NC_000913/709790‑709969

tGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGccc                                                                                         >  1:1221702/1‑93 (MQ=255)

 gaagaaCTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCt                                                                                        >  1:510802/1‑93 (MQ=255)

 gaagaaCTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCt                                                                                        >  1:1195437/1‑93 (MQ=255)

  aagaaCTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTg                                                                                       >  2:4739700/1‑93 (MQ=255)

    gaaCTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGtt                                                                                     >  2:1035403/1‑93 (MQ=255)

    gaaCTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGtt                                                                                     >  1:1048063/1‑93 (MQ=255)

      aCTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTAc                                                                                   >  2:2017728/1‑93 (MQ=255)

       cTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACt                                                                                  >  1:4590791/1‑93 (MQ=255)

         tCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTaa                                                                                <  1:4466582/93‑1 (MQ=255)

          ccAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAg                                                                               >  1:1134570/1‑93 (MQ=255)

           cAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGc                                                                              >  2:2465563/1‑93 (MQ=255)

            aGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCg                                                                             >  1:1698795/1‑93 (MQ=255)

             ggATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCgg                                                                            >  2:3518094/1‑93 (MQ=255)

               aTTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGtt                                                                          >  1:2820861/1‑93 (MQ=255)

                tttttAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTCAGCAGTTc                                                                         >  2:593859/1‑93 (MQ=255)

                 ttttAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCa                                                                        >  2:1703790/1‑93 (MQ=255)

                    tAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACtt                                                                     <  1:1761092/93‑1 (MQ=255)

                     aaaCGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTc                                                                    >  1:3966565/1‑93 (MQ=255)

                       aCGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCAt                                                                  >  1:662477/1‑93 (MQ=255)

                       aCGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGGATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCAt                                                                  >  2:1577472/1‑93 (MQ=255)

                          aGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTag                                                               >  1:3852036/1‑93 (MQ=255)

                          aGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTag                                                               >  1:3883005/1‑93 (MQ=255)

                           ggAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTaga                                                              <  1:3548347/93‑1 (MQ=255)

                           ggAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTaga                                                              <  1:92738/93‑1 (MQ=255)

                              aGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTagaaga                                                           >  2:647884/1‑93 (MQ=255)

                               gCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAgaagaa                                                          <  1:4190611/93‑1 (MQ=255)

                                cGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGaagaag                                                         >  2:2355749/1‑93 (MQ=255)

                                 gTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGt                                                        >  2:3277245/1‑93 (MQ=255)

                                 gTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGt                                                        >  1:2030755/1‑93 (MQ=255)

                                  ttACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTg                                                       >  2:2684281/1‑93 (MQ=255)

                                  ttACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTg                                                       >  2:708100/1‑93 (MQ=255)

                                   tACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGa                                                      >  2:1891215/1‑93 (MQ=255)

                                     cTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACa                                                    <  1:2499662/93‑1 (MQ=255)

                                      tttCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAg                                                   >  2:3167863/1‑93 (MQ=255)

                                       ttCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGa                                                  <  2:2799128/93‑1 (MQ=255)

                                       ttCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGa                                                  <  2:3610906/93‑1 (MQ=255)

                                       tgcAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGa                                                  <  2:2105238/91‑1 (MQ=255)

                                        tCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACAAAGCGGTTCACTTCATTAGAAGAAGTGACAGaa                                                 <  2:2309602/93‑1 (MQ=255)

                                         cAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAAt                                                <  2:584800/93‑1 (MQ=255)

                                         cAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAAt                                                <  1:134615/93‑1 (MQ=255)

                                         cAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAAt                                                >  2:1635323/1‑93 (MQ=255)

                                         cAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAAt                                                >  2:3844854/1‑93 (MQ=255)

                                          aGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAAtt                                               <  2:3757048/93‑1 (MQ=255)

                                          aGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAAtt                                               <  2:1178704/93‑1 (MQ=255)

                                           ggCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAAttt                                              >  1:3439537/1‑93 (MQ=255)

                                           ggCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAAttt                                              <  2:2001187/93‑1 (MQ=255)

                                            gCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTg                                             >  2:2973883/1‑93 (MQ=255)

                                             ccAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGc                                            <  2:509887/93‑1 (MQ=255)

                                             ccAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGc                                            >  1:3676696/1‑93 (MQ=255)

                                              cAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGct                                           >  1:4706043/1‑93 (MQ=255)

                                                gCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGctct                                         <  2:1040019/93‑1 (MQ=255)

                                                gCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGctct                                         >  2:113446/1‑93 (MQ=255)

                                                gCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGctct                                         <  2:1362672/93‑1 (MQ=255)

                                                 ctttctttAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCtt                                        >  1:4372978/1‑93 (MQ=255)

                                                 ctttctttAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCtt                                        <  2:619424/93‑1 (MQ=255)

                                                 ctttctttAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCtt                                        <  2:1277957/93‑1 (MQ=255)

                                                  tttctttAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTg                                       >  1:24686/1‑93 (MQ=255)

                                                    tctttAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTgag                                     <  2:887068/93‑1 (MQ=255)

                                                    tctttAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTgag                                     <  1:4298046/93‑1 (MQ=255)

                                                    tctttAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTgag                                     <  1:420602/93‑1 (MQ=255)

                                                     ctttAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTgaga                                    >  2:3411098/1‑93 (MQ=255)

                                                      tttAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGAt                                   <  1:3184182/93‑1 (MQ=255)

                                                         aGGGCGGTTTTTTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTGCACTTCATTAGAAGAAGTGACAGAAGTTGCTCGTGATATAAt                                >  2:4197868/1‑93 (MQ=255)

                                                         aGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGANAGAATTTGCTCTTGAGATAAt                                >  1:1261993/1‑93 (MQ=255)

                                                           ggCGGTATTGTTATCAGTCATGCGGAATTTGCCCTGTTACTAAGCGGTTGTCTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGc                              <  1:1325716/93‑1 (MQ=255)

                                                           ggCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGc                              >  1:4142362/1‑93 (MQ=255)

                                                           ggCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGc                              >  2:3999546/1‑93 (MQ=255)

                                                              ggTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTa                           >  2:3043381/1‑93 (MQ=255)

                                                              ggTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTa                           >  2:3775878/1‑93 (MQ=255)

                                                              ggTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTa                           >  2:2508742/1‑93 (MQ=255)

                                                              ggTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTa                           >  1:4255582/1‑93 (MQ=255)

                                                              ggTATGGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTa                           <  1:2618770/93‑1 (MQ=255)

                                                                tATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATc                         >  2:3263331/1‑93 (MQ=255)

                                                                 aTTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGATCACTTCATTAGAAGAAGCGACAGAATTTGCTCTTGAGATAATGCGTATCa                        >  2:3096094/1‑93 (MQ=255)

                                                                  ttGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCAt                       <  2:249614/93‑1 (MQ=255)

                                                                     ttATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATtat                    <  2:1644263/93‑1 (MQ=255)

                                                                     ttATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATtat                    <  1:3323281/93‑1 (MQ=255)

                                                                      tATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATtata                   <  2:2656630/93‑1 (MQ=255)

                                                                       aTCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAg                  >  1:3868689/1‑93 (MQ=255)

                                                                        tCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGa                 <  2:4361655/93‑1 (MQ=255)

                                                                         cAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGaa                >  2:4052238/1‑93 (MQ=255)

                                                                          aGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAAt               <  1:4226248/93‑1 (MQ=255)

                                                                           gTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAAtt              >  2:2200764/1‑93 (MQ=255)

                                                                           gTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAAtt              <  1:1854856/93‑1 (MQ=255)

                                                                           gTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAAtt              >  1:2887119/1‑93 (MQ=255)

                                                                            tCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTg             <  1:3535586/93‑1 (MQ=255)

                                                                            tCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTg             <  1:293266/93‑1 (MQ=255)

                                                                              aTGCGGAATCTGCCCTGTTGCTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGcc           >  2:3535220/1‑93 (MQ=255)

                                                                                   gAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGcc      <  2:1537210/93‑1 (MQ=255)

                                                                                   gAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGcc      <  1:3607154/93‑1 (MQ=255)

                                                                                   gAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGcc      <  1:4495833/93‑1 (MQ=255)

                                                                                    aaTCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCt     <  1:2795242/93‑1 (MQ=255)

                                                                                    aaTCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCt     <  1:3044303/93‑1 (MQ=255)

                                                                                    aaTCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCt     >  2:1895315/1‑93 (MQ=255)

                                                                                     aTCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTa    >  2:762858/1‑93 (MQ=255)

                                                                                      tCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTaa   >  2:4627271/1‑93 (MQ=255)

                                                                                      tCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTaa   >  1:2778250/1‑93 (MQ=255)

                                                                                       cTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTaaa  <  1:2276245/93‑1 (MQ=255)

                                                                                       cTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTaaa  <  2:3908109/93‑1 (MQ=255)

                                                                                       cTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTaaa  <  1:25588/93‑1 (MQ=255)

TGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAA  >  NC_000913/709790‑709969

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: `‑`