BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,008,990	C→T	R109C (CGT→TGT)	yigM →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,008,990	0	C	T	98.5%	706.1 / ‑2.4	199	R109C (CGT→TGT)	yigM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (2/1); new base T (94/102); total (96/103)
Fisher's exact test for biased strand distribution p-value = 6.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.84e-01

Only 100 of 199 total aligned reads displayed.

ATCTCTACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGG  >  NC_000913/4008901‑4009079

aTCTCTACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTc                                                                                        <  2:804801/93‑1 (MQ=255)

 tctcTACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCt                                                                                       <  1:874832/93‑1 (MQ=255)

 tctcTACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCt                                                                                       <  2:3544740/93‑1 (MQ=255)

  ctctACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTg                                                                                      <  2:4553468/93‑1 (MQ=255)

  ctctACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTg                                                                                      <  1:3972523/93‑1 (MQ=255)

   tctACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTgc                                                                                     <  1:3133539/93‑1 (MQ=255)

     tACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTgcgc                                                                                   >  2:1295739/1‑93 (MQ=255)

      aCCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCt                                                                                  >  1:4435545/1‑93 (MQ=255)

       ccTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTg                                                                                 >  1:1869628/1‑93 (MQ=255)

        cTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTgg                                                                                <  1:4055459/93‑1 (MQ=255)

        cTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTgg                                                                                >  1:1522382/1‑93 (MQ=255)

         tGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTggg                                                                               >  1:562196/1‑93 (MQ=255)

         tGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTggg                                                                               <  2:1771403/93‑1 (MQ=255)

         tGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTggg                                                                               <  1:3131652/93‑1 (MQ=255)

         tGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTggg                                                                               <  1:3151636/93‑1 (MQ=255)

         tGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTggg                                                                               <  2:2763230/93‑1 (MQ=255)

         tGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTggg                                                                               <  1:1846734/93‑1 (MQ=255)

         tGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTggg                                                                               >  1:1828754/1‑93 (MQ=255)

          gACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTgggg                                                                              >  2:3311517/1‑93 (MQ=255)

           gcGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGc                                                                             <  2:1282461/92‑1 (MQ=255)

             ggTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTa                                                                           <  2:464067/93‑1 (MQ=255)

              gTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTACGCGCTGTCTGCGCTGGGGCTAt                                                                          >  1:4569334/1‑93 (MQ=255)

              gTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTAt                                                                          <  2:442284/93‑1 (MQ=255)

              gTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTAt                                                                          >  2:156845/1‑93 (MQ=255)

                ttCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGc                                                                        <  2:2537518/93‑1 (MQ=255)

                   cGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCtt                                                                     <  1:3729727/93‑1 (MQ=255)

                    gAGCTGCTGCTGTTCACCGTGCTGGCGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCttt                                                                    >  1:2579019/1‑93 (MQ=255)

                     aGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTa                                                                   <  2:1407986/93‑1 (MQ=255)

                     aGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTa                                                                   <  2:2558505/93‑1 (MQ=255)

                     aGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTa                                                                   >  1:1426557/1‑93 (MQ=255)

                      gctgctgctGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAg                                                                  >  1:2427321/1‑93 (MQ=255)

                       ctgctgctgTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAgc                                                                 <  2:316951/93‑1 (MQ=255)

                           tgctgTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCt                                                             >  2:1953166/1‑93 (MQ=255)

                             ctgTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTg                                                           >  2:4413574/1‑93 (MQ=255)

                              tgTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGc                                                          <  2:2406233/93‑1 (MQ=255)

                               gTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATTACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCt                                                         >  2:4147483/1‑93 (MQ=255)

                                ttCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTg                                                        >  1:788432/1‑93 (MQ=255)

                                 tCACCGTGCTGCCGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTgg                                                       <  1:3082567/93‑1 (MQ=255)

                                  cACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGc                                                      >  1:351398/1‑93 (MQ=255)

                                  cACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGc                                                      >  1:3250334/1‑93 (MQ=255)

                                   aCCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCg                                                     <  1:3593169/93‑1 (MQ=255)

                                   aCCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCg                                                     >  1:2798662/1‑93 (MQ=255)

                                   aCCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCg                                                     >  1:2687247/1‑93 (MQ=255)

                                   aCCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCg                                                     >  2:2624049/1‑93 (MQ=255)

                                    ccGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCgg                                                    >  1:3462545/1‑93 (MQ=255)

                                     cGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCggg                                                   >  1:1195068/1‑92 (MQ=255)

                                     cGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGt                                                   <  2:3897582/93‑1 (MQ=255)

                                      gTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTg                                                  >  2:1075674/1‑93 (MQ=255)

                                      gTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTg                                                  >  1:3818161/1‑93 (MQ=255)

                                       tGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGa                                                 <  1:2154067/93‑1 (MQ=255)

                                       tGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGCGa                                                 >  1:487103/1‑93 (MQ=255)

                                         cTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGAtt                                               >  2:3477124/1‑93 (MQ=255)

                                         cTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGAtt                                               <  1:1097305/93‑1 (MQ=255)

                                          tGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTg                                              >  1:3347268/1‑93 (MQ=255)

                                           gACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTgg                                             >  2:3775517/1‑93 (MQ=255)

                                           gACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTgg                                             <  1:2269930/93‑1 (MQ=255)

                                            aCGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGt                                            <  1:4632778/93‑1 (MQ=255)

                                              gccgcTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGc                                          >  1:3804524/1‑93 (MQ=255)

                                              gccgcTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGc                                          >  1:3033062/1‑93 (MQ=255)

                                                 gcTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCgg                                       >  2:4350689/1‑93 (MQ=255)

                                                 gcTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCgg                                       >  1:1057487/1‑93 (MQ=255)

                                                  ctctACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCggg                                      >  2:3835066/1‑93 (MQ=255)

                                                  ctctACATCACGCTGATTTATGACATCATGAGTAAGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCggg                                      >  1:4033310/1‑93 (MQ=255)

                                                  ctctACATCACGCTGATTTATGACATCATGAGTAAGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCggg                                      <  2:2504007/93‑1 (MQ=255)

                                                   tctACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGa                                     <  2:4494398/93‑1 (MQ=255)

                                                   tctACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGa                                     <  2:4080415/93‑1 (MQ=255)

                                                    ctACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGat                                    <  2:1263843/93‑1 (MQ=255)

                                                     tACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGatt                                   >  2:2554805/1‑93 (MQ=255)

                                                      aCATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGatta                                  <  2:427634/93‑1 (MQ=255)

                                                        aTCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGattatt                                >  1:1888358/1‑93 (MQ=255)

                                                        aTCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGattatt                                <  2:3218386/93‑1 (MQ=255)

                                                          cACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCg                              >  2:2389857/1‑93 (MQ=255)

                                                           aCGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGc                             >  1:2834201/1‑93 (MQ=255)

                                                            cGCTGATTTATTACATCATGAGTAAGCGCTGTCTGCGCTGGAGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCt                            <  1:2039583/93‑1 (MQ=255)

                                                             gCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTa                           <  2:3713047/93‑1 (MQ=255)

                                                              cTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTAt                          <  2:676551/93‑1 (MQ=255)

                                                              cTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTAt                          >  2:4324827/1‑93 (MQ=255)

                                                               tGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATg                         >  1:3389527/1‑93 (MQ=255)

                                                               tGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATg                         >  2:1670960/1‑93 (MQ=255)

                                                                gATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGa                        >  1:2503991/1‑93 (MQ=255)

                                                                gATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGa                        >  2:2356959/1‑93 (MQ=255)

                                                                gATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGa                        <  1:1474578/93‑1 (MQ=255)

                                                                  tttATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATc                      <  1:1271467/93‑1 (MQ=255)

                                                                   ttATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCa                     <  2:1845427/93‑1 (MQ=255)

                                                                     aTGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAgg                   >  2:2304205/1‑93 (MQ=255)

                                                                      tGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGt                  <  1:3215635/93‑1 (MQ=255)

                                                                       gACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTc                 <  1:1987920/93‑1 (MQ=255)

                                                                       gACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTc                 <  1:689959/93‑1 (MQ=255)

                                                                        aCATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCa                >  2:182439/1‑93 (MQ=255)

                                                                        aCATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCa                <  2:1031612/93‑1 (MQ=255)

                                                                         catcatGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCAc               >  1:1680141/1‑93 (MQ=255)

                                                                         catcatGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCAc               <  1:4542103/93‑1 (MQ=255)

                                                                          atcatGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCAcc              >  2:221829/1‑93 (MQ=255)

                                                                           tcatGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCg             >  2:1492382/1‑93 (MQ=255)

                                                                                aGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCAt        >  2:1579411/1‑93 (MQ=255)

                                                                                 gTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCAtt       <  1:1284710/93‑1 (MQ=255)

                                                                                  tAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCAttt      <  2:1676976/93‑1 (MQ=255)

                                                                                   aaGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTc     >  1:2866196/1‑93 (MQ=255)

                                                                                    aGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCAtttat    >  2:2506609/1‑91 (MQ=255)

                                                                                    aGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCt    <  1:1958310/93‑1 (MQ=255)

ATCTCTACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTGG  >  NC_000913/4008901‑4009079

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