BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,209,664	(G)_5→6	coding (96/834 nt)	mcrA →	McrA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,209,659	1	.	G	96.8%	400.4 / 5.4	126	coding (91/834 nt)	mcrA	McrA
Reads supporting (aligned to +/- strand): ref base . (1/3); new base G (62/60); total (63/63)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.25e-01

Only 100 of 149 total aligned reads displayed.

TTATGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCG  >  NC_000913/1209567‑1209746

ttatGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGt                                                                                            >  2:1261270/1‑93 (MQ=255)

ttatGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGt                                                                                            >  1:3267875/1‑93 (MQ=255)

ttatGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGt                                                                                            >  2:810696/1‑93 (MQ=255)

 tatGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑g                                                                                         <  1:2955691/93‑1 (MQ=255)

 tatGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑g                                                                                         >  2:275047/1‑93 (MQ=255)

 tatGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑g                                                                                         >  2:4182788/1‑93 (MQ=255)

 tatGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑g                                                                                         <  2:450419/93‑1 (MQ=255)

 tatGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑g                                                                                         >  1:462795/1‑93 (MQ=255)

 tatGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑g                                                                                         >  1:786502/1‑93 (MQ=255)

 tatGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑g                                                                                         <  1:490956/93‑1 (MQ=255)

  atGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑gg                                                                                        <  1:4684739/93‑1 (MQ=255)

   tGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑ggg                                                                                       <  1:1881149/93‑1 (MQ=255)

   tGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑ggg                                                                                       <  1:2525644/93‑1 (MQ=255)

        gTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑ggggggcc                                                                                  >  2:2676106/1‑90 (MQ=255)

        gTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑ggggggcc                                                                                  <  2:1352098/93‑4 (MQ=255)

         tttttGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑ggggggccg                                                                                 >  1:4562167/1‑89 (MQ=255)

         tttttGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGCGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑ggggggccg                                                                                 >  2:3059552/1‑89 (MQ=255)

           tttGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑GGGGG‑CCGGGt                                                                              <  1:1456014/93‑1 (MQ=255)

            ttGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCgggg                                                                              >  2:2248214/1‑92 (MQ=255)

              gATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGa                                                                            >  1:2018481/1‑93 (MQ=255)

                taataaTGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACa                                                                          >  2:841862/1‑93 (MQ=255)

                   taatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGaa                                                                       >  1:2218541/1‑93 (MQ=255)

                   taatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGaa                                                                       >  2:33253/1‑93 (MQ=255)

                    aatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGaaa                                                                      >  2:1889004/1‑93 (MQ=255)

                    aatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGaaa                                                                      >  2:1743171/1‑93 (MQ=255)

                    aatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGACGTG‑GGGGG‑CCGGGCGACAGcaa                                                                      >  2:2595296/1‑93 (MQ=255)

                     atGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAAc                                                                     <  2:3174135/93‑1 (MQ=255)

                     atGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAAc                                                                     >  2:3583258/1‑93 (MQ=255)

                      tGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACa                                                                    >  1:920276/1‑93 (MQ=255)

                      tGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACa                                                                    >  1:663706/1‑93 (MQ=255)

                       ggAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACaa                                                                   >  2:3309840/1‑93 (MQ=255)

                         aaTTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAg                                                                 >  1:1666583/1‑93 (MQ=255)

                         aaTTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAg                                                                 <  1:3719594/93‑1 (MQ=255)

                          aTTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGa                                                                >  1:1659384/1‑93 (MQ=255)

                           ttGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGAt                                                               <  2:1293993/93‑1 (MQ=255)

                           ttGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑GGGGGTCCGGGTGACAGAAACAAAGAt                                                               <  1:4584985/93‑1 (MQ=255)

                            tGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGAtt                                                              <  2:3084528/93‑1 (MQ=255)

                              aaCTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTAc                                                            <  1:1836319/93‑1 (MQ=255)

                              aaCTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTAc                                                            >  1:1621528/1‑93 (MQ=255)

                               aCTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACa                                                           >  2:3769539/1‑93 (MQ=255)

                                cTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTAACAGAAACAAAGATTACaa                                                          <  2:1415898/93‑1 (MQ=255)

                                 tGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGG‑GCGGGTGACAGAAACAAAGATTACaa                                                          <  1:3161789/93‑1 (MQ=255)

                                 tGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGTCAGAAACAAAGATTACAat                                                         <  1:3498265/93‑1 (MQ=255)

                                  gAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAata                                                        >  2:559178/1‑93 (MQ=255)

                                  gAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAata                                                        >  1:3819596/1‑93 (MQ=255)

                                  gAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAata                                                        >  2:1180360/1‑93 (MQ=255)

                                    aaGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATc                                                      <  1:2687880/93‑1 (MQ=255)

                                    aaGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATc                                                      >  2:3546972/1‑93 (MQ=255)

                                       cTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGct                                                   <  2:2284093/93‑1 (MQ=255)

                                       cTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGct                                                   <  2:3067840/93‑1 (MQ=255)

                                        tGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGctc                                                  >  2:1308289/1‑93 (MQ=255)

                                        tGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGctc                                                  <  1:1050361/93‑1 (MQ=255)

                                        tGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGctc                                                  >  1:1379492/1‑93 (MQ=255)

                                         gAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGctct                                                 <  2:4683298/93‑1 (MQ=255)

                                          aGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAGACAAAGATTACAATATCGCTCtt                                                <  1:3274437/93‑1 (MQ=255)

                                           gtgtTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTg                                               >  1:2105403/1‑93 (MQ=255)

                                            tgtTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTgt                                              >  1:1231915/1‑92 (MQ=255)

                                             gtTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGAt                                             <  1:3028548/93‑1 (MQ=255)

                                              ttCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGAtt                                            <  2:4316071/93‑1 (MQ=255)

                                              ttCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGAtt                                            >  1:2465817/1‑93 (MQ=255)

                                               tCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATta                                           >  2:4553970/1‑93 (MQ=255)

                                               tCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATta                                           <  1:3419895/93‑1 (MQ=255)

                                                cGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATtat                                          <  2:1937066/93‑1 (MQ=255)

                                                cGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATtat                                          <  2:512049/93‑1 (MQ=255)

                                                  aTAGGTGAAGAGGATGTTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTatat                                        <  2:441952/93‑1 (MQ=255)

                                                  aTAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTatat                                        >  2:700481/1‑93 (MQ=255)

                                                   tAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATc                                       >  1:1422921/1‑93 (MQ=255)

                                                      gTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCAtt                                    <  1:125008/93‑1 (MQ=255)

                                                       tGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTg                                   <  2:2869358/93‑1 (MQ=255)

                                                        gAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGa                                  >  2:2483888/1‑93 (MQ=255)

                                                         aaGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGaa                                 >  2:4262554/1‑93 (MQ=255)

                                                          agagGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGCTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAAc                                >  2:3884827/1‑93 (MQ=255)

                                                           gagGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACg                               >  1:4064643/1‑93 (MQ=255)

                                                            agGATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACgg                              >  2:780945/1‑93 (MQ=255)

                                                             ggATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggt                             <  2:4719908/93‑1 (MQ=255)

                                                              gATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggt                             <  2:4622843/93‑1 (MQ=255)

                                                              gATGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggtt                            >  2:423982/1‑93 (MQ=255)

                                                                tGGTGTTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttgg                          >  1:4180347/1‑93 (MQ=255)

                                                                  gtgtTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttggtt                        >  1:2416021/1‑93 (MQ=255)

                                                                  gtgtTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttggtt                        <  2:3861222/93‑1 (MQ=255)

                                                                   tgtTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGgttggttg                       >  2:3283802/1‑93 (MQ=255)

                                                                    gtTTATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGa                      <  1:642849/93‑1 (MQ=255)

                                                                      ttATGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGAtt                    <  1:749182/93‑1 (MQ=255)

                                                                         tGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTg                 <  2:355902/93‑1 (MQ=255)

                                                                         tGGTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTg                 >  1:2639011/1‑93 (MQ=255)

                                                                          ggTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTgg                >  1:1998981/1‑93 (MQ=255)

                                                                           gTCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGgt               <  1:4656858/93‑1 (MQ=255)

                                                                            tCTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGgtg              >  2:1232641/1‑93 (MQ=255)

                                                                             cTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGgtgt             <  1:3737838/93‑1 (MQ=255)

                                                                             cTAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGgtgt             >  1:2314329/1‑93 (MQ=255)

                                                                              tAATCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATAGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTa            >  1:2891104/1‑93 (MQ=255)

                                                                               aaTCCTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTAt           >  2:3062760/1‑93 (MQ=255)

                                                                                  ccTTGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATccc        <  2:1625911/93‑1 (MQ=255)

                                                                                    ttGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATcacca      >  2:2731404/1‑89 (MQ=255)

                                                                                    ttGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCaa      <  1:2859231/93‑1 (MQ=255)

                                                                                    ttGAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCaa      <  1:2859232/93‑1 (MQ=255)

                                                                                      gAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGt    <  1:4289056/93‑1 (MQ=255)

                                                                                      gAGTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGt    >  1:851335/1‑93 (MQ=255)

                                                                                       ggTCGTG‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATACCAAGTc   <  2:2425891/92‑1 (MQ=255)

                                                                                       aGTCGT‑‑GGGGGTCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGGATACCAAGTc   <  2:2425889/93‑1 (MQ=255)

TTATGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑‑GGGGG‑CCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCG  >  NC_000913/1209567‑1209746

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