BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	653,656	G→A	R191H (CGC→CAC)	citB →	CitB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	653,656	0	G	A	99.3%	490.9 / ‑5.4	140	R191H (CGC→CAC)	citB	CitB
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (65/74); total (66/74)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

Only 100 of 148 total aligned reads displayed.

CCTTGACGCTAAACGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAAC  >  NC_000913/653566‑653745

ccTTGACGCTAAACGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCcaca                                                                                         >  1:2818704/1‑93 (MQ=255)

ccTTGACGCTAAACGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCcaca                                                                                         >  2:4673506/1‑93 (MQ=255)

    gACGCTAAACGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCAcaccac                                                                                     <  1:2903646/93‑1 (MQ=255)

    gACGCTAAACGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCAcaccac                                                                                     <  2:3701125/93‑1 (MQ=255)

     aCGCTAAACGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCGCACCACt                                                                                    >  1:1553052/1‑93 (MQ=255)

     aCGCTAAACGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACt                                                                                    <  1:3394286/93‑1 (MQ=255)

     aCGCTAAACGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACt                                                                                    >  2:2742928/1‑93 (MQ=255)

      cGCTAAACGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTg                                                                                   >  1:739865/1‑93 (MQ=255)

       gCTAAACGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGc                                                                                  >  1:749957/1‑93 (MQ=255)

        cTAAACGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGcc                                                                                 <  1:4639095/93‑1 (MQ=255)

         tAAACGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCa                                                                                <  2:1993206/93‑1 (MQ=255)

         tAAACGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCa                                                                                >  1:832653/1‑93 (MQ=255)

            aCGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGc                                                                             >  2:163300/1‑93 (MQ=255)

            aCGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGc                                                                             >  2:4115296/1‑93 (MQ=255)

               ctgTGCGAAAACTATTTAAAGATCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACGGACCATCAGCCACACCACTGCCAGGCGtt                                                                          <  2:891953/91‑1 (MQ=255)

               cgGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGtt                                                                          >  1:3186336/1‑93 (MQ=255)

                  tGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATc                                                                       <  2:3474010/93‑1 (MQ=255)

                   gCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCCCTGACCATCAGCCACACCACTGCCAGGCGTTATCt                                                                      <  1:520286/93‑1 (MQ=255)

                   gCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCATGCGTTATCt                                                                      <  2:3141822/93‑1 (MQ=255)

                   gCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCt                                                                      <  1:4723645/93‑1 (MQ=255)

                    cGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCtt                                                                     <  1:2538945/93‑1 (MQ=255)

                    cGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCtt                                                                     <  1:2890344/93‑1 (MQ=255)

                    cGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCtt                                                                     <  1:4536340/93‑1 (MQ=255)

                      aaaaCTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGa                                                                   <  1:3256751/93‑1 (MQ=255)

                       aaaCTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGaa                                                                  >  1:2007209/1‑93 (MQ=255)

                       aaaCTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGaa                                                                  >  1:3877045/1‑93 (MQ=255)

                        aaCTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAat                                                                 <  1:656323/93‑1 (MQ=255)

                        aaCTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAat                                                                 >  1:2833989/1‑93 (MQ=255)

                        aaCTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAat                                                                 <  1:2213504/93‑1 (MQ=255)

                         aCTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAata                                                                >  1:34445/1‑93 (MQ=255)

                         aCTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAata                                                                >  2:1670865/1‑93 (MQ=255)

                          cTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAatat                                                               >  2:2971592/1‑93 (MQ=255)

                           tGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATAtt                                                              >  2:4709835/1‑93 (MQ=255)

                           tGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATAtt                                                              >  2:1735488/1‑93 (MQ=255)

                            gTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTg                                                             >  2:2516592/1‑93 (MQ=255)

                            gTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTg                                                             >  2:515236/1‑93 (MQ=255)

                               tAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTgcgc                                                          >  2:244281/1‑93 (MQ=255)

                                aaaGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGcc                                                         >  1:2076677/1‑93 (MQ=255)

                                aaaGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGcc                                                         <  1:1233543/93‑1 (MQ=255)

                                 aaGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCa                                                        >  2:2097390/1‑93 (MQ=255)

                                   gagCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAgc                                                      >  1:3597943/1‑93 (MQ=255)

                                       cTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAgccgcc                                                  >  1:2249402/1‑93 (MQ=255)

                                        tGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCa                                                 <  2:3686969/93‑1 (MQ=255)

                                        tGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCa                                                 <  1:3166136/93‑1 (MQ=255)

                                          gtgtGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATc                                               <  1:747966/93‑1 (MQ=255)

                                           tgtgCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCt                                              <  1:1085820/93‑1 (MQ=255)

                                           ggtgCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCt                                              <  1:1612385/92‑1 (MQ=255)

                                              gCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGat                                           <  1:1056340/93‑1 (MQ=255)

                                              gCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGat                                           <  1:3218743/93‑1 (MQ=255)

                                               cAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGatt                                          >  1:3311571/1‑93 (MQ=255)

                                                aaCATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGatta                                         <  2:1649141/93‑1 (MQ=255)

                                                aaCATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGatta                                         <  1:3674496/93‑1 (MQ=255)

                                                aaCATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGatta                                         >  1:1996009/1‑93 (MQ=255)

                                                  cATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGattatt                                       <  2:1597586/93‑1 (MQ=255)

                                                     aCGGCAGAAACGGTGGCGCAGGCACGGACCATCAGCCACACCACTGCCAGGCGTTATCATGAATATAGCGCCAGCCGCCAGCTGATTATTGcg                                    >  1:4199302/1‑92 (MQ=255)

                                                        gCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGaa                                 >  2:4238582/1‑93 (MQ=255)

                                                        gCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGaa                                 >  2:878198/1‑93 (MQ=255)

                                                          aGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAAt                               <  1:1228274/93‑1 (MQ=255)

                                                           gAAACGTGGGCGCCGGCACTGACCATCAGCCACACCTCTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAAtt                              <  1:1897004/93‑1 (MQ=255)

                                                           gAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAAtt                              >  1:4352167/1‑93 (MQ=255)

                                                             aaCGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGt                            >  1:3752571/1‑93 (MQ=255)

                                                             aaCGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCAGCCATCTGATTATTGCTGAAATTGt                            <  1:1066179/93‑1 (MQ=255)

                                                              aCGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGCAATTGtt                           >  2:2643381/1‑93 (MQ=255)

                                                               cGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTc                          <  2:4502733/93‑1 (MQ=255)

                                                                gggggCGCCGGCACTGACCATCAGCCACCCCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCa                         <  1:2709421/90‑1 (MQ=255)

                                                                ggTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCa                         <  1:2226962/93‑1 (MQ=255)

                                                                 gTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCAc                        >  2:4308405/1‑93 (MQ=255)

                                                                 gTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCAc                        <  2:4036286/93‑1 (MQ=255)

                                                                 gTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCAc                        <  2:3336758/93‑1 (MQ=255)

                                                                  tGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACg                       <  2:1432484/93‑1 (MQ=255)

                                                                   ggCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACgg                      <  2:2110770/93‑1 (MQ=255)

                                                                   ggCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACgg                      <  2:4707034/93‑1 (MQ=255)

                                                                    gcgcAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGc                     <  2:2263390/93‑1 (MQ=255)

                                                                    gcgcAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGc                     <  2:1889529/93‑1 (MQ=255)

                                                                       cAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCaaa                  <  2:749502/93‑1 (MQ=255)

                                                                        aGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAg                 <  1:1894090/93‑1 (MQ=255)

                                                                          gCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGtt               <  1:4136850/93‑1 (MQ=255)

                                                                          gCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGtt               <  1:1396301/93‑1 (MQ=255)

                                                                          gCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGtt               <  1:1854828/93‑1 (MQ=255)

                                                                           cACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTg              <  2:3280557/93‑1 (MQ=255)

                                                                           cACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTg              <  2:2779068/93‑1 (MQ=255)

                                                                            ccTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTgg             <  2:3485297/92‑1 (MQ=255)

                                                                            aCTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTgg             <  1:4669524/93‑1 (MQ=255)

                                                                            aCTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTgg             >  2:1101050/1‑93 (MQ=255)

                                                                              tGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCa           >  1:30324/1‑93 (MQ=255)

                                                                              tGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCa           >  2:1922744/1‑93 (MQ=255)

                                                                               gACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAACGTTGGCAg          >  1:35386/1‑93 (MQ=255)

                                                                               gACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAg          >  1:1888839/1‑93 (MQ=255)

                                                                                 ccATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAAATGTTCACGGCAAAGTTGGCAGAc        <  2:689509/93‑1 (MQ=255)

                                                                                  cATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGAcc       >  1:2407786/1‑93 (MQ=255)

                                                                                  cATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGAcc       <  1:4089109/93‑1 (MQ=255)

                                                                                   aTCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACca      <  2:1858179/93‑1 (MQ=255)

                                                                                   aTCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACca      >  2:1423224/1‑93 (MQ=255)

                                                                                    tCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACcac     <  2:497001/93‑1 (MQ=255)

                                                                                     cAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACcaca    <  1:2692824/93‑1 (MQ=255)

                                                                                      aGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACaa   <  1:1403867/93‑1 (MQ=255)

                                                                                      aGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACaa   >  2:2367472/1‑93 (MQ=255)

                                                                                       gCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAAc  <  1:3486738/93‑1 (MQ=255)

                                                                                       gCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAAc  >  2:1940532/1‑93 (MQ=255)

                                                                                       gCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAAc  >  1:4431107/1‑93 (MQ=255)

CCTTGACGCTAAACGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAAC  >  NC_000913/653566‑653745

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