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	98.8%	261.5 / ‑3.5	84	coding (91/834 nt)	mcrA	McrA
Reads supporting (aligned to +/- strand): ref base . (0/1); new base G (50/33); total (50/34)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 102 total aligned reads displayed.

TTATGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑GGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAA  >  NC_000913/1209567‑1209742

ttatGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGt                                                                                      <  2:372168/93‑1 (MQ=255)

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

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

   tGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑ggg                                                                                  >  1:1128870/1‑93 (MQ=255)

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

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

    gCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑gggg                                                                                 >  1:4459807/1‑93 (MQ=255)

    gCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑gggg                                                                                 >  2:1613519/1‑93 (MQ=255)

    gCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑gggg                                                                                 <  2:2347109/93‑1 (MQ=255)

    gCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑gggg                                                                                 <  2:49268/93‑1 (MQ=255)

     cATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑ggggg                                                                                >  1:1981975/1‑93 (MQ=255)

     cATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑ggggg                                                                                >  1:1033670/1‑93 (MQ=255)

     cATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑ggggg                                                                                <  1:4221388/93‑1 (MQ=255)

      aTGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑GGGGGc                                                                               <  1:3551405/93‑1 (MQ=255)

       tGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑ggggggc                                                                              <  1:3772301/93‑3 (MQ=255)

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

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

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

           tttGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCggg                                                                           <  1:3052724/93‑1 (MQ=255)

           tttGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCggg                                                                           <  2:1713588/93‑1 (MQ=255)

                taataaTGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACa                                                                      >  1:1095542/1‑93 (MQ=255)

                taataaTGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACa                                                                      >  2:3714107/1‑93 (MQ=255)

                 aataatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAg                                                                     >  1:1695361/1‑93 (MQ=255)

                  ataatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGa                                                                    >  2:1693457/1‑93 (MQ=255)

                  ataatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGa                                                                    >  1:2986330/1‑93 (MQ=255)

                  ataatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGa                                                                    >  2:3230999/1‑93 (MQ=255)

                  ataatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGa                                                                    >  1:3160793/1‑93 (MQ=255)

                   taatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGaa                                                                   >  1:4318361/1‑93 (MQ=255)

                   taatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGaa                                                                   >  2:1139816/1‑93 (MQ=255)

                   taatGGAATTGAACTAAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGaa                                                                   <  1:2947094/93‑1 (MQ=255)

                    aatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGaaa                                                                  >  1:4516369/1‑93 (MQ=255)

                    aatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGaaa                                                                  <  2:736196/93‑1 (MQ=255)

                     atGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAAc                                                                 >  1:938829/1‑93 (MQ=255)

                      tGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACa                                                                >  1:2024686/1‑93 (MQ=255)

                       ggAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACaa                                                               >  1:3515551/1‑93 (MQ=255)

                       ggAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACaa                                                               >  2:3895646/1‑93 (MQ=255)

                       ggAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACaa                                                               <  1:1643886/93‑1 (MQ=255)

                        gAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAATaaa                                                              <  2:2747400/93‑1 (MQ=255)

                        gAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACaaa                                                              >  1:3586082/1‑93 (MQ=255)

                        gAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACaaa                                                              <  2:3576611/93‑1 (MQ=255)

                        gAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACaaa                                                              <  2:4083431/93‑1 (MQ=255)

                          aTTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGa                                                            <  1:2734124/93‑1 (MQ=255)

                          aTTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGa                                                            >  1:826041/1‑93 (MQ=255)

                           ttGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGAt                                                           <  1:268073/93‑1 (MQ=255)

                           ttGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGAt                                                           >  1:4595415/1‑93 (MQ=255)

                            tGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGAtt                                                          >  1:2453390/1‑93 (MQ=255)

                             gAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTa                                                         <  1:4058180/93‑1 (MQ=255)

                              aaCTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTAc                                                        <  1:1965348/93‑1 (MQ=255)

                              aaCTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTAc                                                        <  2:2593410/93‑1 (MQ=255)

                              aaCTGAAAGCTGAGTGTTCGATAGGGGAAGAGGATGGAGTTTATGGCCTAATCCTTTAATCGTGGGGGGCCGGGTGACAGAAACAAAGATTAc                                                        >  2:3285549/1‑93 (MQ=255)

                                cTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACaa                                                      >  1:1607921/1‑93 (MQ=255)

                                cTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACaa                                                      >  2:1033969/1‑93 (MQ=255)

                                cTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACaa                                                      >  1:866035/1‑93 (MQ=255)

                                 tGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAat                                                     <  2:256580/93‑1 (MQ=255)

                                 tGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAat                                                     >  1:4159508/1‑93 (MQ=255)

                                  gAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAata                                                    >  2:4329266/1‑93 (MQ=255)

                                  gAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAata                                                    >  2:459429/1‑93 (MQ=255)

                                    aaGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATc                                                  <  2:2416813/93‑1 (MQ=255)

                                     aGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCg                                                 >  2:3175502/1‑93 (MQ=255)

                                      gCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGc                                                >  1:4341384/1‑93 (MQ=255)

                                        tGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGctc                                              >  1:1290585/1‑93 (MQ=255)

                                        tGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGctc                                              >  1:2964767/1‑93 (MQ=255)

                                        tGAGTGTTAGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCGGGGGGATCGAAACAAAGATTACAGTATCGctt                                              >  1:916992/1‑92 (MQ=255)

                                         gAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGctct                                             >  2:2125804/1‑93 (MQ=255)

                                         gAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGctct                                             <  2:2197095/93‑1 (MQ=255)

                                          aGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCtt                                            <  1:1825332/93‑1 (MQ=255)

                                           gtgtTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTg                                           <  2:4027840/93‑1 (MQ=255)

                                            tgtTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGa                                          >  1:4478229/1‑93 (MQ=255)

                                             gtTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGAt                                         <  2:2653939/93‑1 (MQ=255)

                                              ttCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGAtt                                        >  1:3860746/1‑93 (MQ=255)

                                               tCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATta                                       <  1:4645026/93‑1 (MQ=255)

                                               tCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATta                                       <  2:3514471/93‑1 (MQ=255)

                                               tCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATta                                       >  1:2692781/1‑93 (MQ=255)

                                                cGATAGGTGAAGAGTATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATtat                                      <  1:1319064/93‑1 (MQ=255)

                                                 gATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATtatt                                     >  2:463354/1‑92 (MQ=255)

                                                   tAGGTGAAGAGGATGGTGTTTATGGTCTACTCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATc                                   <  2:4604488/93‑1 (MQ=255)

                                                   tAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATc                                   >  1:2058272/1‑93 (MQ=255)

                                                       tGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTg                               <  1:401196/93‑1 (MQ=255)

                                                       tGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTg                               <  2:349525/93‑1 (MQ=255)

                                                        gAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGa                              <  1:2143339/93‑1 (MQ=255)

                                                         aaGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGaa                             >  1:4340125/1‑93 (MQ=255)

                                                         aaGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGaa                             <  2:4035515/93‑1 (MQ=255)

                                                             ggATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggt                         >  1:1449395/1‑93 (MQ=255)

                                                              gATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggtt                        <  1:2743552/93‑1 (MQ=255)

                                                              gATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggtt                        <  2:3836850/93‑1 (MQ=255)

                                                               aTGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttg                       >  1:3397392/1‑93 (MQ=255)

                                                                  gtgtTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttggtt                    <  1:172549/93‑1 (MQ=255)

                                                                  gtgtTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttggtt                    >  2:459911/1‑93 (MQ=255)

                                                                  gtgtTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTAAACggttggtt                    >  2:459912/1‑93 (MQ=255)

                                                                    gtTTATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGa                  >  2:3007296/1‑93 (MQ=255)

                                                                     tttATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGAt                 >  2:135183/1‑93 (MQ=255)

                                                                       tATGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTc               >  1:164864/1‑93 (MQ=255)

                                                                        aTGGTCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCt              >  2:4066909/1‑93 (MQ=255)

                                                                            tCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGgtg          >  2:343032/1‑93 (MQ=255)

                                                                            tCTAATCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGgtg          <  2:1684041/93‑1 (MQ=255)

                                                                               aaTCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTAt       >  2:4077778/1‑93 (MQ=255)

                                                                                 tCCTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATcc     >  1:3614764/1‑93 (MQ=255)

                                                                                   cTTGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCa   >  2:732648/1‑93 (MQ=255)

                                                                                    ttGAGTCGTGGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCaa  >  1:1397814/1‑93 (MQ=255)

                                                                                    ttGAGTCGTGGGGGGCCGGGGGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCaa  >  2:3822163/1‑93 (MQ=255)

TTATGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGT‑GGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAA  >  NC_000913/1209567‑1209742

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