BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,794,501	A→G	S103P (TCT→CCT)	ygaU ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,794,501	0	A	G	98.7%	476.6 / ‑3.1	149	S103P (TCT→CCT)	ygaU	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (75/72); total (76/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 152 total aligned reads displayed.

TTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATATTCCCCACCG  >  NC_000913/2794411‑2794590

  aGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGTCAGGCTTa                                                                                       <  2:767478/93‑1 (MQ=255)

   gCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTaa                                                                                      >  1:3390315/1‑93 (MQ=255)

   gCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTaa                                                                                      >  1:233180/1‑93 (MQ=255)

   gCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTaa                                                                                      >  2:2275426/1‑93 (MQ=255)

   gCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTaa                                                                                      <  1:3923999/93‑1 (MQ=255)

    cATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAAc                                                                                     <  1:1185419/93‑1 (MQ=255)

     aTCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACg                                                                                    >  2:305713/1‑93 (MQ=255)

     aTCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACg                                                                                    >  1:2502482/1‑93 (MQ=255)

     aTCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACg                                                                                    >  2:3588489/1‑93 (MQ=255)

     aTCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACg                                                                                    >  1:247636/1‑93 (MQ=255)

      tCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACgg                                                                                   >  1:787536/1‑93 (MQ=255)

      tCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACgg                                                                                   <  1:2016436/93‑1 (MQ=255)

      tCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACgg                                                                                   >  2:4108462/1‑93 (MQ=255)

       cGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGt                                                                                  >  1:3451498/1‑93 (MQ=255)

         gTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGtat                                                                                >  1:3329111/1‑93 (MQ=255)

          tttATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGtata                                                                               <  2:472014/93‑1 (MQ=255)

          tttATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGtata                                                                               <  1:3841677/93‑1 (MQ=255)

             aTTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATaaaa                                                                            >  1:3423422/1‑93 (MQ=255)

              ttCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAAc                                                                           >  2:152275/1‑93 (MQ=255)

                cGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAActg                                                                         <  1:4234747/93‑1 (MQ=255)

                 gCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAActgg                                                                        >  1:610931/1‑93 (MQ=255)

                 gCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAActgg                                                                        >  1:3726867/1‑93 (MQ=255)

                  cTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAActggc                                                                       <  2:1821048/93‑1 (MQ=255)

                   ttCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAActggct                                                                      <  2:238506/93‑1 (MQ=255)

                   ttCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAActggct                                                                      >  2:2002486/1‑93 (MQ=255)

                      gAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACtggctggc                                                                   >  1:3839764/1‑93 (MQ=255)

                      gAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACtggctggc                                                                   <  1:733620/93‑1 (MQ=255)

                      gAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCATAACGGTATAAAACttgctggc                                                                   >  1:1743028/1‑93 (MQ=255)

                       aaGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCa                                                                  >  1:2878051/1‑93 (MQ=255)

                        aGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAg                                                                 >  2:2540953/1‑93 (MQ=255)

                         gATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGt                                                                <  2:1150522/93‑1 (MQ=255)

                         gATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGt                                                                >  2:4235821/1‑93 (MQ=255)

                           ttttATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTgg                                                              >  2:4412675/1‑93 (MQ=255)

                           ttttATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGGGTCGCCCGGCTTAACGGTATAAAACTGGCTGTCAGCgg                                                              >  2:4192538/1‑93 (MQ=255)

                            tttATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGc                                                             >  2:3944184/1‑93 (MQ=255)

                            tttATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGc                                                             <  1:4539912/93‑1 (MQ=255)

                             ttATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCt                                                            >  2:1264705/1‑93 (MQ=255)

                             ttATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCt                                                            >  2:494603/1‑93 (MQ=255)

                             ttATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCt                                                            >  1:56629/1‑93 (MQ=255)

                             ttATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCt                                                            <  2:4678309/93‑1 (MQ=255)

                               aTTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTgg                                                          >  2:2409446/1‑93 (MQ=255)

                                ttGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGgt                                                         <  2:387514/93‑1 (MQ=255)

                                ttGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGgt                                                         >  1:4740725/1‑93 (MQ=255)

                                 tGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGgtg                                                        <  2:4231421/93‑1 (MQ=255)

                                 tGTACAGATTAGCGTTACCGTAGACCTGTTCGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGgtg                                                        <  1:1974711/93‑1 (MQ=255)

                                  gTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGgtgt                                                       >  2:3897815/1‑93 (MQ=255)

                                     cAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGc                                                    >  2:695527/1‑93 (MQ=255)

                                      aGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGGGTCGcc                                                   >  2:2725592/1‑93 (MQ=255)

                                       gATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGGGTCGCCg                                                  >  2:1014977/1‑93 (MQ=255)

                                        aTTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGt                                                 <  2:2214731/93‑1 (MQ=255)

                                        aTTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGt                                                 >  1:265911/1‑93 (MQ=255)

                                           aGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGtttt                                              >  1:696172/1‑93 (MQ=255)

                                              gTTACCGTAGACCTGTTTGGAAATGGCGCTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCAc                                           >  2:1628424/1‑93 (MQ=255)

                                              gTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGCGTCGCCGGTTTCAc                                           >  1:3786625/1‑93 (MQ=255)

                                               ttACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCAcc                                          >  1:3252691/1‑93 (MQ=255)

                                               ttACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGGGTGGCCGTTTTCAcc                                          >  1:1159913/1‑93 (MQ=255)

                                                tACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCt                                         <  2:502220/93‑1 (MQ=255)

                                                  ccGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGa                                       <  1:2236285/93‑1 (MQ=255)

                                                  ccGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGGGTCGCCGTTTTCACCCGa                                       >  2:1724137/1‑93 (MQ=255)

                                                   cGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGat                                      >  1:125612/1‑93 (MQ=255)

                                                    gTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGatc                                     >  1:1289381/1‑93 (MQ=255)

                                                    gTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGatc                                     <  2:4090985/93‑1 (MQ=255)

                                                     tAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGatca                                    <  2:889664/93‑1 (MQ=255)

                                                      aGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGatcat                                   >  1:1327537/1‑93 (MQ=255)

                                                       gACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGatcatc                                  >  1:192219/1‑93 (MQ=255)

                                                       gACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGatcatc                                  <  2:2044718/93‑1 (MQ=255)

                                                       gACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCAACTGatcatc                                  >  1:99425/1‑93 (MQ=255)

                                                         ccTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGa                                <  2:957775/93‑1 (MQ=255)

                                                         ccTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGa                                >  1:3373748/1‑93 (MQ=255)

                                                          cTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGac                               >  2:2532754/1‑93 (MQ=255)

                                                          cTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGac                               <  2:4028998/93‑1 (MQ=255)

                                                          cTGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGac                               <  2:703825/93‑1 (MQ=255)

                                                           tGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGaca                              <  1:2341346/93‑1 (MQ=255)

                                                           tGTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGaca                              >  1:1403774/1‑93 (MQ=255)

                                                            gTTTGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGacac                             <  2:483755/93‑1 (MQ=255)

                                                              ttGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTg                           >  1:2162969/1‑93 (MQ=255)

                                                              ttGGAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTg                           <  1:3953418/93‑1 (MQ=255)

                                                                ggAAATGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGc                         <  1:3582463/93‑1 (MQ=255)

                                                                   aaTGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAAt                      <  1:1817652/93‑1 (MQ=255)

                                                                   aaTGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAAt                      <  2:2559609/93‑1 (MQ=255)

                                                                    aTGGCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATa                     <  1:4276263/93‑1 (MQ=255)

                                                                     aggCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCGATAc                    <  1:1240517/92‑1 (MQ=255)

                                                                      ggCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATAcc                   <  1:1908759/93‑1 (MQ=255)

                                                                       gCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCg                  >  2:4552722/1‑93 (MQ=255)

                                                                       gCACTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCg                  <  2:4498075/93‑1 (MQ=255)

                                                                         aCTCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGa                <  1:2552181/93‑1 (MQ=255)

                                                                           tCAGAGTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGaaa              <  2:2773691/93‑1 (MQ=255)

                                                                             atagTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAata            <  2:2512272/91‑1 (MQ=255)

                                                                             agagTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAata            >  2:3294839/1‑93 (MQ=255)

                                                                              gagTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCCATACCGGAAatat           >  2:1111982/1‑93 (MQ=255)

                                                                              gagTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAatat           >  1:1385185/1‑93 (MQ=255)

                                                                               agTGTCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATAtt          >  2:1828090/1‑93 (MQ=255)

                                                                                gtgtCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATATTc         >  1:1857270/1‑93 (MQ=255)

                                                                                gtgtCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATATTc         <  1:4023707/93‑1 (MQ=255)

                                                                                 tgtCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATATTcc        <  2:3229895/93‑1 (MQ=255)

                                                                                  ggcGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATATTccc       <  2:4015896/91‑1 (MQ=255)

                                                                                   tCGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATATTcccc      <  1:1724830/93‑1 (MQ=255)

                                                                                   tCGCCAGGCCTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATATTcccc      <  1:3161230/93‑1 (MQ=255)

                                                                                    cGCCAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATATTccccc     >  2:4186994/1‑92 (MQ=255)

                                                                                       cAGGCTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATATTCCCCAccc  >  1:2572926/1‑92 (MQ=255)

TTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTTTTCACCTGATCATCGACACTGGCAATACCGGAAATATTCCCCACCG  >  NC_000913/2794411‑2794590

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: `‑`