BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,157,555	T→C	K536E (AAG→GAG)	ygiRQ ←	Fe‑S oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,157,555	0	T	C	90.1%	743.0 / 59.6	253	K536E (AAG→GAG)	ygiRQ	Fe‑S oxidoreductase
Reads supporting (aligned to +/- strand): ref base T (11/14); new base C (109/119); total (120/133)
Fisher's exact test for biased strand distribution p-value = 8.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.94e-01

Only 100 of 264 total aligned reads displayed.

AATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATA  >  NC_000913/3157465‑3157645

aaTAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTGAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTcca                                                                                          >  2:3559951/1‑93 (MQ=255)

aaTAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTcca                                                                                          >  1:1047636/1‑93 (MQ=255)

aaTAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTcca                                                                                          >  1:1703382/1‑93 (MQ=255)

aaTAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTcca                                                                                          <  1:856224/93‑1 (MQ=255)

 aTAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTccat                                                                                         <  2:5633857/93‑1 (MQ=255)

    ggggATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCcatcat                                                                                      >  1:1207531/1‑93 (MQ=255)

     gggATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCatcatc                                                                                     >  2:6483825/1‑93 (MQ=255)

     gggATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCatcatc                                                                                     <  2:5726964/93‑1 (MQ=255)

     gggATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCatcatc                                                                                     >  1:504731/1‑93 (MQ=255)

      ggATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCt                                                                                    <  1:4381245/93‑1 (MQ=255)

      gaaTCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCt                                                                                    <  1:1637622/91‑1 (MQ=255)

        aTCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTc                                                                                  <  1:2719141/93‑1 (MQ=255)

        aTCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTc                                                                                  >  1:4814365/1‑93 (MQ=255)

            gATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATa                                                                              >  2:437517/1‑93 (MQ=255)

              tACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAAc                                                                            >  1:1802306/1‑93 (MQ=255)

               tgtGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACg                                                                           <  2:4512694/91‑1 (MQ=255)

                  gCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGcc                                                                        >  1:1952145/1‑93 (MQ=255)

                  gCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGcc                                                                        >  2:3531116/1‑93 (MQ=255)

                   ctctTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGccc                                                                       >  1:1996713/1‑93 (MQ=255)

                   ctctTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGccc                                                                       <  2:712419/93‑1 (MQ=255)

                    tctTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCct                                                                      <  1:5146757/93‑1 (MQ=255)

                    tctTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCct                                                                      <  1:492745/93‑1 (MQ=255)

                    tctTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCct                                                                      <  2:5214022/93‑1 (MQ=255)

                     ctTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCctt                                                                     <  2:2858864/93‑1 (MQ=255)

                      tttACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCcttc                                                                    >  1:2816746/1‑93 (MQ=255)

                      tttACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCcttc                                                                    <  1:3973914/93‑1 (MQ=255)

                        tACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCcttctt                                                                  <  1:4234694/93‑1 (MQ=255)

                         aCCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCttcttc                                                                 <  2:2878679/93‑1 (MQ=255)

                         aCCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCttcttc                                                                 >  1:1806264/1‑93 (MQ=255)

                          cctgcctgTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCg                                                                <  2:3105863/93‑1 (MQ=255)

                           ctgcctgTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCgg                                                               <  1:4511769/93‑1 (MQ=255)

                           ctgcctgTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCgg                                                               >  1:3579016/1‑93 (MQ=255)

                             gcctgTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTa                                                             <  1:2177365/93‑1 (MQ=255)

                             gcctgTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTa                                                             >  2:4342338/1‑93 (MQ=255)

                              gctgTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCCGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTAt                                                            <  2:134804/92‑1 (MQ=255)

                              cctgTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTAt                                                            >  2:5405063/1‑93 (MQ=255)

                                tgTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGt                                                          >  1:1113855/1‑93 (MQ=255)

                                tgTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGt                                                          <  2:3838840/93‑1 (MQ=255)

                                tgTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGt                                                          >  1:4645475/1‑93 (MQ=255)

                                 gTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGtt                                                         >  2:61893/1‑93 (MQ=255)

                                  ttttGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTc                                                        <  1:6587162/93‑1 (MQ=255)

                                   tttGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGTCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTcc                                                       <  2:4614039/93‑1 (MQ=255)

                                   tttGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTcc                                                       <  1:2583921/93‑1 (MQ=255)

                                     tGAGTAAGTATCGAACCGCTCTTTAAAGCGGACATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCCTCTTCGGTATGATCCgg                                                     >  1:525184/1‑93 (MQ=255)

                                     tGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCgg                                                     >  1:4446551/1‑93 (MQ=255)

                                     tGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCgg                                                     <  1:3741711/93‑1 (MQ=255)

                                     tGAGTAAGTATCAAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCgg                                                     <  1:4269575/93‑1 (MQ=255)

                                      gAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCggg                                                    <  2:17839/93‑1 (MQ=255)

                                       agtaagtaTCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCgggg                                                   <  2:452828/93‑1 (MQ=255)

                                       agtaagtaTCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCgggg                                                   >  1:6450627/1‑93 (MQ=255)

                                        gtaagtaTCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGc                                                  >  2:2998085/1‑93 (MQ=255)

                                         taagtaTCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCa                                                 >  1:5671361/1‑93 (MQ=255)

                                         taagtaTCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCa                                                 >  1:1332412/1‑93 (MQ=255)

                                         taagtaTCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCa                                                 <  2:1670670/93‑1 (MQ=255)

                                             taTCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCt                                             <  1:5535597/93‑1 (MQ=255)

                                              aTCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGATCCGGGGCAATCtt                                            >  1:3035210/1‑93 (MQ=255)

                                               tCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGCTAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTc                                           >  1:3902741/1‑93 (MQ=255)

                                               tCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTc                                           >  2:38772/1‑93 (MQ=255)

                                               tCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTc                                           >  2:4768054/1‑93 (MQ=255)

                                               tCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTc                                           >  2:6088232/1‑93 (MQ=255)

                                                   aCAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGAt                                       >  2:4031993/1‑93 (MQ=255)

                                                     aGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATaa                                     <  2:5288207/93‑1 (MQ=255)

                                                     aGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATaa                                     >  1:1165197/1‑93 (MQ=255)

                                                     aGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATaa                                     >  2:6416787/1‑93 (MQ=255)

                                                      gCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAAc                                    <  1:2400169/93‑1 (MQ=255)

                                                      gCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAAc                                    <  1:5559650/93‑1 (MQ=255)

                                                       ctctTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAAcc                                   <  1:6106645/93‑1 (MQ=255)

                                                       ctctTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAAcc                                   >  2:1477137/1‑93 (MQ=255)

                                                       ctctTTAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAAcc                                   <  1:5934470/93‑1 (MQ=255)

                                                          tttAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAAccgcc                                >  2:3579235/1‑93 (MQ=255)

                                                          tttAAAGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAAccgcc                                >  1:2786504/1‑93 (MQ=255)

                                                           ttAAAGCGGTCATAGCTACCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAAccgccg                               <  1:5810244/93‑1 (MQ=255)

                                                            tAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGa                              <  2:4564341/93‑1 (MQ=255)

                                                             aaaGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGAc                             >  1:4133638/1‑93 (MQ=255)

                                                              aaGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACg                            <  2:596863/93‑1 (MQ=255)

                                                              aaGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACg                            >  1:2125311/1‑93 (MQ=255)

                                                              aaGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACg                            <  2:1916580/93‑1 (MQ=255)

                                                              aaGCGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACg                            <  2:5859833/93‑1 (MQ=255)

                                                               aGCGGTCATAGCTGCCCATGCGCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGt                           <  2:2374269/93‑1 (MQ=255)

                                                                gCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTg                          >  2:5116632/1‑93 (MQ=255)

                                                                 cGGTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGa                         >  1:985768/1‑93 (MQ=255)

                                                                  ggTCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGAt                        <  2:6134545/93‑1 (MQ=255)

                                                                    tCATAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATgg                      <  2:2535825/93‑1 (MQ=255)

                                                                    tCATAGCTGCCCATGCCCGGCTCCATCATCTTCCATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATgg                      >  1:2915091/1‑93 (MQ=255)

                                                                      aTAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGt                    <  1:3136206/93‑1 (MQ=255)

                                                                      aTAGCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGt                    >  2:4417048/1‑93 (MQ=255)

                                                                         gCTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGc                 <  2:2241793/93‑1 (MQ=255)

                                                                          cTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGcc                <  1:410106/93‑1 (MQ=255)

                                                                          cTGCCCATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGcc                >  2:4786911/1‑93 (MQ=255)

                                                                           tGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCa               >  2:566556/1‑93 (MQ=255)

                                                                             cccATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGt             <  2:3257310/93‑1 (MQ=255)

                                                                               cATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTc           <  1:1446902/93‑1 (MQ=255)

                                                                               cATGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGCCGCCAGTTc           >  1:3865681/1‑93 (MQ=255)

                                                                                aTGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCt          >  2:3455652/1‑93 (MQ=255)

                                                                                aTGCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCt          >  2:4307751/1‑93 (MQ=255)

                                                                                  gCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCttt        <  2:2500261/93‑1 (MQ=255)

                                                                                  gCCCGGCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCttt        <  1:6209400/93‑1 (MQ=255)

                                                                                      ggCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGata    <  2:2287719/93‑1 (MQ=255)

                                                                                       gCTCCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGatat   >  1:6044975/1‑93 (MQ=255)

                                                                                        cTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGAtata  <  2:2064370/93‑1 (MQ=255)

AATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATAACCGCCGACGTGATGGGTCGCCAGTTCTTTGATATA  >  NC_000913/3157465‑3157645

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