BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,717,416	G→A	M221I (ATG→ATA)	ung →	uracil‑DNA glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,717,416	0	G	A	99.6%	1095.6 / NA	262	M221I (ATG→ATA)	ung	uracil‑DNA glycosylase
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (121/139); minor base C (0/1); total (121/141)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 262 total aligned reads displayed.

TCGCCGCTTTCGGCGCATCGTGGATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATGCCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAGGCTTTATTCTGACGCCACC  >  NC_000913/2717318‑2717513

tCGCCGCTTTCGGCGCATCGTGGATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATAcc                                                                                                 <  1:3950507/101‑1 (MQ=255)

tCGCCGCTTTCGGCGCATCGTGGATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATAcc                                                                                                 <  2:3950994/101‑1 (MQ=255)

   ccgcTTTCGGCGCATCGTGGATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGt                                                                                              >  1:3951216/1‑101 (MQ=255)

    cgcTTTCGGCGCATCGTGGATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTa                                                                                             <  1:3950897/101‑1 (MQ=255)

    cgcTTTCGGCGCATCGTGGATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTa                                                                                             <  1:3950890/101‑1 (MQ=255)

     gcTTTCGGCGCATCGTGGATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTAt                                                                                            >  1:3951221/1‑101 (MQ=255)

     gcTTTCGGCGCATCGTGGATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTAt                                                                                            <  1:3950876/101‑1 (MQ=255)

     gcTTTCGGCGCATCGTGGATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTAt                                                                                            <  1:3950774/101‑1 (MQ=255)

      cTTTCGGCGCATCGTGGATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTAtt                                                                                           <  1:3950747/101‑1 (MQ=255)

      cTTTCGGCGCATCGTGGATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTAtt                                                                                           <  2:3950705/101‑1 (MQ=255)

       tttCGGCGCATCGTGGATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTa                                                                                          <  1:3950661/101‑1 (MQ=255)

        ttCGGCGCATCGTGGATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTAc                                                                                         >  1:3951226/1‑101 (MQ=255)

        ttCGGCGCATCGTGGATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTAc                                                                                         <  1:3950998/101‑1 (MQ=255)

          cGGCGCATCGTGGATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCg                                                                                       >  1:3951229/1‑101 (MQ=255)

           ggCGCATCGTGGATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCgg                                                                                      <  1:3950730/101‑1 (MQ=255)

             cgcATCGTGGATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCa                                                                                    <  1:3950753/101‑1 (MQ=255)

               cATCGTGGATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCaga                                                                                  >  1:3951234/1‑101 (MQ=255)

                aTCGTGGATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCagag                                                                                 >  1:3951237/1‑101 (MQ=255)

                  cGTGGATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCagagag                                                                               <  2:3950686/101‑1 (MQ=255)

                   gTGGATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAgagt                                                                              >  2:3951239/1‑101 (MQ=255)

                    tGGATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCagagag                                                                               >  2:3951241/1‑99 (MQ=255)

                    tGGATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAgagtg                                                                             <  1:3950767/101‑1 (MQ=255)

                     ggATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAgagtgc                                                                            <  1:3950944/101‑2 (MQ=255)

                     ggATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAgagtga                                                                            >  1:3951244/1‑101 (MQ=255)

                      gATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAgagtgag                                                                           >  1:3951249/1‑101 (MQ=255)

                       aTTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATCCCAGTATTACCGGCAGAtagtgagt                                                                          <  2:3950840/101‑1 (MQ=255)

                        ttCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTa                                                                         >  2:3951251/1‑101 (MQ=255)

                        ttCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTa                                                                         <  2:3950827/101‑1 (MQ=255)

                        ttCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTa                                                                         >  1:3951252/1‑101 (MQ=255)

                         tCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCagag                                                                                 >  2:3951259/1‑92 (MQ=255)

                         tCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTaa                                                                        >  2:3951257/1‑101 (MQ=255)

                         tCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTaa                                                                        <  2:3950756/101‑1 (MQ=255)

                          cTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTaaa                                                                       <  1:3950966/101‑1 (MQ=255)

                          cTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTaaa                                                                       <  2:3950879/101‑1 (MQ=255)

                           tttGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAAt                                                                      >  2:3951263/1‑101 (MQ=255)

                           tttGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAAt                                                                      >  1:3951261/1‑101 (MQ=255)

                            ttGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAAtt                                                                     <  2:3950904/101‑1 (MQ=255)

                            ttGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAAtt                                                                     <  2:3950788/101‑1 (MQ=255)

                            ttGGCTGCAACCATTTTGTGCTGGCAAATCAATGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAAtt                                                                     <  1:3950830/101‑1 (MQ=255)

                            ttGGCTGCAACCATTTTGTGCTGGCAAATAAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAAtt                                                                     >  2:3951265/1‑101 (MQ=255)

                             tGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAAttt                                                                    >  1:3951268/1‑101 (MQ=255)

                               ggtGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCa                                                                                    <  2:3950657/81‑1 (MQ=255)

                               gCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGc                                                                  <  2:3950856/101‑1 (MQ=255)

                                cTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCg                                                                 <  2:3950833/101‑1 (MQ=255)

                                 tGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCgg                                                                <  2:3950801/101‑1 (MQ=255)

                                    aaCCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGa                                                             >  2:3951272/1‑101 (MQ=255)

                                       cATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAAt                                                          >  1:3951275/1‑101 (MQ=255)

                                         ttttGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGc                                                        <  1:3950702/101‑1 (MQ=255)

                                         ttttGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGc                                                        <  1:3950798/101‑1 (MQ=255)

                                          tttGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGcc                                                       >  2:3951279/1‑101 (MQ=255)

                                           ttGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATTCCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCg                                                      <  2:3950727/101‑1 (MQ=255)

                                           ttGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAAttg                                                                    >  1:3951280/1‑86 (MQ=255)

                                            tgtgCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCgg                                                     <  1:3950920/101‑1 (MQ=255)

                                            tgtgCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCgg                                                     >  1:3951284/1‑101 (MQ=255)

                                              tgCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGAt                                                   <  1:3950913/101‑1 (MQ=255)

                                              tgCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGAt                                                   <  2:3950761/101‑1 (MQ=255)

                                               gCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATg                                                  >  2:3951286/1‑101 (MQ=255)

                                                cTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATgg                                                 >  1:3951287/1‑101 (MQ=255)

                                                  ggCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCa                                               <  1:3950688/101‑1 (MQ=255)

                                                   gCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCag                                              <  1:3951036/101‑1 (MQ=255)

                                                    cAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAAt                                                          >  2:3951291/1‑88 (MQ=255)

                                                     aaaTCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAAt                                                          >  1:3951294/1‑87 (MQ=255)

                                                     aaaTCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCagag                                            <  2:3951118/101‑1 (MQ=255)

                                                      aaTCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATgg                                                 >  2:3951296/1‑95 (MQ=255)

                                                       aTCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGtt                                          <  1:3950970/101‑1 (MQ=255)

                                                        tCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTg                                         <  1:3951163/101‑1 (MQ=255)

                                                         cAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTGc                                        <  1:3950849/101‑1 (MQ=255)

                                                         cAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTGc                                        >  1:3951298/1‑101 (MQ=255)

                                                          aGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCgggg                                                              >  2:3951302/1‑78 (MQ=255)

                                                              gCTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCag                                              >  1:3951303/1‑90 (MQ=255)

                                                               cTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCg                                  <  2:3950635/101‑1 (MQ=255)

                                                               cTGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCg                                  <  1:3951078/101‑1 (MQ=255)

                                                                tGGAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCgg                                 >  2:3951307/1‑101 (MQ=255)

                                                                 ggAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGcc                                                       >  1:3951309/1‑78 (MQ=255)

                                                                  gAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCt                               <  2:3950862/101‑1 (MQ=255)

                                                                  gAACAACGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCt                               <  2:3950873/101‑1 (MQ=255)

                                                                   aacaacGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCgg                                 <  2:3950825/98‑1 (MQ=255)

                                                                    acaacGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGaaa                                                           >  2:3951313/1‑71 (MQ=255)

                                                                    acaacGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGa                             >  1:3951312/1‑101 (MQ=255)

                                                                       acGTGGCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGAttt                          <  2:3951024/101‑1 (MQ=255)

                                                                           ggCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTg                              >  1:3951318/1‑93 (MQ=255)

                                                                            gCGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAg                     <  1:3950901/101‑1 (MQ=255)

                                                                             cGAGACGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCagag                                            >  2:3951320/1‑77 (MQ=255)

                                                                               agaCGCCGATTGACTGGATACCAGTATTACCGGCAGAgagtg                                                                             >  2:3951327/1‑42 (MQ=255)

                                                                               agaCGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGa                                                             >  1:3951326/1‑58 (MQ=255)

                                                                               agaCGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAGGCt                  >  1:3951323/1‑101 (MQ=255)

                                                                               agaCGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAGGCt                  <  2:3951187/101‑1 (MQ=255)

                                                                               agaCGCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAGGCt                  <  1:3950858/101‑1 (MQ=255)

                                                                                   gCCGATTGACTGGATACCAGTATTACCGGCAGAgagt                                                                              >  2:3951330/1‑37 (MQ=255)

                                                                                   gCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCa                                               <  1:3950884/68‑1 (MQ=255)

                                                                                   gCCGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAGGCTTTAt              <  2:3950777/101‑1 (MQ=255)

                                                                                     cGATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAGGCTTTATTc            >  1:3951333/1‑101 (MQ=255)

                                                                                      gATTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAGGCTTTATTCt           <  1:3950864/101‑1 (MQ=255)

                                                                                       aTTGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAAt                                                          >  1:3951336/1‑53 (MQ=255)

                                                                                        ttGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAGGCTTTATTCTGa         <  1:3950978/101‑1 (MQ=255)

                                                                                         tGACTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAGGCTTTATTCTGAc        <  2:3950926/101‑1 (MQ=255)

                                                                                           aCTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGTGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAGGCTTTATTCTGACGc      >  2:3951339/1‑101 (MQ=255)

                                                                                           aCTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAAt                                                          >  1:3951340/1‑49 (MQ=255)

                                                                                           aCTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAGGCTTTATTCTGACGc      <  2:3950972/101‑1 (MQ=255)

                                                                                            cTGGATACCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAGGCTTTATTCTGACg       <  2:3950934/99‑1 (MQ=255)

TCGCCGCTTTCGGCGCATCGTGGATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATGCCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAGGCTTTATTCTGACGCCACC  >  NC_000913/2717318‑2717513

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