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.6%	799.2 / ‑5.8	220	R191H (CGC→CAC)	citB	CitB
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (108/111); total (108/112)
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.85e-01

Only 100 of 231 total aligned reads displayed.

CCTTGACGCTAAACGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAAC  >  NC_000913/653566‑653745

 cTTGACGCTAAACGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCacac                                                                                        >  2:4288125/1‑93 (MQ=255)

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

   tGACGCTAAACGCGGTGCGAAAACTGTTTAAGGAGCCTGGTGTGCAAAATACGGCAGAAACGGTGGCGCAGGCACTTACCATCAGCCAcacca                                                                                      <  1:3403707/93‑1 (MQ=255)

   tGACGCTAAACGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCAcacca                                                                                      >  2:2582322/1‑93 (MQ=255)

    gACGCTAAACGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCAcaccac                                                                                     >  1:2402405/1‑93 (MQ=255)

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

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

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

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

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

            aCGGGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGCCACTGACCATCAGCCACACCACTGCCAGGc                                                                             <  1:6163655/93‑1 (MQ=255)

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

              gcgGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGt                                                                           >  2:5288673/1‑93 (MQ=255)

              gcgGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGt                                                                           <  1:3039028/93‑1 (MQ=255)

              gcgGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGt                                                                           <  1:2933134/93‑1 (MQ=255)

              gcgGTGCGAAAACTGTTTAAAGAGCCTGGAGTGCCACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACAACACTGCCAGGCGt                                                                           >  2:3180906/1‑93 (MQ=255)

                ggTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTa                                                                         >  1:4940847/1‑93 (MQ=255)

                ggTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTa                                                                         <  2:1255269/93‑1 (MQ=255)

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

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

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

                     gAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATATTg                                                                    >  1:4455045/1‑93 (MQ=255)

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

                       aaaCTGTTTAAAGAGCCTGGTGTGCAACATACGGCGGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGaa                                                                  <  1:1614833/93‑1 (MQ=255)

                       aaaCTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGaa                                                                  <  2:3968305/93‑1 (MQ=255)

                       aaaCTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGaa                                                                  <  2:4082729/93‑1 (MQ=255)

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

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

                         aCTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAata                                                                <  1:5777881/93‑1 (MQ=255)

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

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

                            gTTTAAAGAGCCTGGTGTGCAGCATACGGCAGAAACGGTGGCGCAGGCACTGCCCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTg                                                             <  1:3105530/93‑1 (MQ=255)

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

                             gttAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTgc                                                            <  1:674319/92‑1 (MQ=255)

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

                                     gCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAgccg                                                    >  1:1342640/1‑93 (MQ=255)

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

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

                                        tGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCa                                                 >  2:3230450/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                                                  cATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGattatt                                       >  2:287121/1‑93 (MQ=255)

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

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

                                                   aTACGGCAGAAACGGTGGCGCAGGCGCTGACCATCAGCCACACCACTGCAAGGCCTTATCTAGAATATTGCGCCAGCCGCCATCTGATTATTg                                      <  2:5777986/93‑1 (MQ=255)

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

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

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

                                                      cGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGCATATTGCGCCAGCCGCCATCTGATTATTGCTg                                   >  2:1757226/1‑93 (MQ=255)

                                                      cGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTg                                   >  1:3919276/1‑93 (MQ=255)

                                                       ggCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGa                                  <  2:2830123/93‑1 (MQ=255)

                                                       ggCAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGa                                  <  2:112899/93‑1 (MQ=255)

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

                                                         cAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGaaa                                <  2:10109/93‑1 (MQ=255)

                                                         cAGAAACGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGaaa                                <  1:2522272/93‑1 (MQ=255)

                                                            aaaCGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTg                             >  1:4131117/1‑93 (MQ=255)

                                                              aCGGTGGCGCAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGtt                           >  2:4722161/1‑93 (MQ=255)

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

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

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

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

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

                                                                 gTGGCGCAGGCACTGACCATCAGCCACACCACGGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCAc                        <  1:1578205/93‑1 (MQ=255)

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

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

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

                                                                   ggCGCAGGCACTTACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACgg                      <  1:4780226/93‑1 (MQ=255)

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

                                                                       cAGGCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGCTCACGGCaaa                  >  1:2109555/1‑93 (MQ=255)

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

                                                                         ggCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGt                <  1:1992961/93‑1 (MQ=255)

                                                                         ggCACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGt                <  2:2599175/93‑1 (MQ=255)

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

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

                                                                           cACTGACCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTg              >  1:1921357/1‑93 (MQ=255)

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

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

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

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

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

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

                                                                                aCCATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAGGTTGGCAGa         <  2:5327312/93‑1 (MQ=255)

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

                                                                                 ccATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGAc        <  2:1375064/93‑1 (MQ=255)

                                                                                 ccATCAGCCACACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGAc        <  2:2131015/93‑1 (MQ=255)

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

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

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

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

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

CCTTGACGCTAAACGCGGTGCGAAAACTGTTTAAAGAGCCTGGTGTGCAACATACGGCAGAAACGGTGGCGCAGGCACTGACCATCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAAC  >  NC_000913/653566‑653745

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