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.2%	879.4 / ‑4.2	250	R191H (CGC→CAC)	citB	CitB
Reads supporting (aligned to +/- strand): ref base G (0/2); new base A (128/120); total (128/122)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-01

Only 100 of 267 total aligned reads displayed.

CCTTGACGCTAAACGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAAC  >  NC_000913/653566‑653745

ccTTGACGCTAAACGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCcaca                                                                                         <  2:6599537/93‑1 (MQ=255)

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

 cTTGACGCTAAACGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCacac                                                                                        <  1:4951030/93‑1 (MQ=255)

  ttGACGCTAAACGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCAcacc                                                                                       >  1:7143368/1‑93 (MQ=255)

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

      cGCTAAACGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTg                                                                                   <  1:1036973/93‑1 (MQ=255)

       gCTAAACGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGGCACACCACTGc                                                                                  >  2:1351712/1‑93 (MQ=255)

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

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

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

          aaaCGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAg                                                                               <  2:1056974/93‑1 (MQ=255)

          aaaCGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAg                                                                               <  2:2115820/93‑1 (MQ=255)

            aCGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGc                                                                             <  2:3837011/93‑1 (MQ=255)

                ggTGCGAAAACTGTTTAAAGAGCCTGGGGTGCAAGATACGGCAGAAACGGTGGCGCAGGCACTGTCCATTAGCCACACCACTGCCAGGCGTTa                                                                         >  2:6074982/1‑93 (MQ=255)

                 gTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGGGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTAt                                                                        >  1:2648010/1‑93 (MQ=255)

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

                  tGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATc                                                                       >  1:978272/1‑93 (MQ=255)

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

                   gCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCt                                                                      >  1:3818021/1‑93 (MQ=255)

                    ggAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGCCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCtt                                                                     <  1:2844962/92‑1 (MQ=255)

                    cGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCtt                                                                     >  2:4438348/1‑93 (MQ=255)

                     gAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTg                                                                    >  1:883995/1‑93 (MQ=255)

                      aaaaCTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGa                                                                   >  1:4261779/1‑93 (MQ=255)

                       aaccTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGaa                                                                  <  2:3994650/90‑1 (MQ=255)

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

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

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

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

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

                         aCTATTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAata                                                                <  2:4349167/93‑1 (MQ=255)

                          cTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAatat                                                               <  2:472433/93‑1 (MQ=255)

                           tGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATAtt                                                              <  2:2022111/93‑1 (MQ=255)

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

                             tttAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTgc                                                            >  2:4123503/1‑93 (MQ=255)

                             tttAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTgc                                                            <  2:1544586/93‑1 (MQ=255)

                             tttAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTgc                                                            <  2:4669353/93‑1 (MQ=255)

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

                               tAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTgcgc                                                          <  2:5735617/93‑1 (MQ=255)

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

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

                                   gagCCTGGTGTGCAACATACGGCAGAAACGGGGGCGCAGGCACTGAACATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAgc                                                      >  2:5545906/1‑93 (MQ=255)

                                    agCCTGGTGTGCAGCATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAgcc                                                     <  1:1121482/93‑1 (MQ=255)

                                    agCCTGGTGTGCAACATACGGCAGAATCGGTGGCGCAGGAACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAgcc                                                     >  2:3281350/1‑93 (MQ=255)

                                    agCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAgcc                                                     <  1:5797261/93‑1 (MQ=255)

                                      ccTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAgccgc                                                   >  1:7708949/1‑93 (MQ=255)

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

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

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

                                         ggTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCAt                                                <  2:6170723/93‑1 (MQ=255)

                                          gtgtGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATc                                               >  2:4874642/1‑93 (MQ=255)

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

                                             tgCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGa                                            >  1:4507618/1‑93 (MQ=255)

                                              gCAACATACGGCAGATACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGat                                           <  1:7543254/93‑1 (MQ=255)

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

                                              gCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGat                                           >  2:6497417/1‑93 (MQ=255)

                                              gCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGat                                           >  1:331866/1‑93 (MQ=255)

                                              gCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGat                                           >  2:6786113/1‑93 (MQ=255)

                                              gCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGat                                           >  2:3538614/1‑93 (MQ=255)

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

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

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

                                                 aCATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGattat                                        >  2:5020156/1‑93 (MQ=255)

                                                   aTACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTg                                      >  1:250494/1‑93 (MQ=255)

                                                    tACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGc                                     >  1:3441372/1‑93 (MQ=255)

                                                     aCGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCt                                    <  2:2566416/93‑1 (MQ=255)

                                                       ggCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGa                                  >  2:3323348/1‑93 (MQ=255)

                                                          aGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAAt                               >  1:4600297/1‑93 (MQ=255)

                                                           gAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAAtt                              <  1:4451255/93‑1 (MQ=255)

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

                                                           gAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATAATGCGCCAGCCGCCATCTGATTATTGCTGAAAtt                              >  2:2238269/1‑93 (MQ=255)

                                                            aaaCGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTg                             <  1:4148620/93‑1 (MQ=255)

                                                            aaaCGGTGGCGAAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTg                             <  2:5417250/93‑1 (MQ=255)

                                                             aaCGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGt                            <  2:4568083/93‑1 (MQ=255)

                                                             aaCGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGt                            <  2:5317103/93‑1 (MQ=255)

                                                                ggTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCa                         >  2:5190088/1‑93 (MQ=255)

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

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

                                                                      gcAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCaa                   >  1:2396807/1‑93 (MQ=255)

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

                                                                        aGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAg                 >  1:2551207/1‑93 (MQ=255)

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

                                                                             cTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGc            >  1:940773/1‑93 (MQ=255)

                                                                             cTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGc            >  2:5949523/1‑93 (MQ=255)

                                                                              tGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCa           <  2:3086423/93‑1 (MQ=255)

                                                                               gACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAg          <  1:783120/93‑1 (MQ=255)

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

                                                                                 ccATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGAc        >  2:7284200/1‑93 (MQ=255)

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

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

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

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

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

                                                                                    tCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACcac     >  1:136104/1‑93 (MQ=255)

                                                                                    tCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACcac     >  2:189806/1‑93 (MQ=255)

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

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

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

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

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

                                                                                       gCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAAc  >  1:7811606/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: `‑`